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SECTION 1 


SUMMARY 


This report describes the digital computer program developed to study 
the vibration response of a coupled rotor/bifilar/airframe system. 

The theoretical development of the rotor/airframe system equations of 
motion is provided in Reference 1 while the fuselage and bifilar absorber 
equations of motion are discussed in Appendix D. 

The modular block approach used in the make-up of this computer program 
is described in Section 2. Section 3 provides descriptions of the input 
data needed to run the rotor and bifilar absorber analyses. Sample output 
formats are presented and discussed in Section 4. The results for four 
test cases, which use the major logic paths of the computer program, are 
presented in Section 5. In Section 6, the overall program structure is 
discussed in detail, including the segmentation procedure (overlay) needed 
to run the program on the NASA CDC computer system, the routine flow dia- 
grams and a list of the COMMON blocks. Finally, in Section 7* the Fortran 
subroutines are described in detail. 

The work was conducted under NASA contract NASI- 15612, "Bifilar Analysis 
Study" . 



SECTION 2 


PROGRAM DESCRIPTION 


The bi filar absorber analysis was developed to provide the engineer 
with an analytical tool capable of rapid parametric vibration evaluation of 
the entire helicopter. A block diagram of the vibration analysis program 
is shown in Figure 1, where it can be seen that a modular approach has been 
adopted to form the main program. Modularization is achieved when each 
component block outputs the mass, damping and stiffness matrices or force 
vectors with the pertinent degrees-of- freedom (d.o.f. ) including the 6 
d.o.f. of the attachment point, i.e. 6 d.o.f. at the hub for rotor or 
bifilar absorber attachment, or 6 d.o.f. in the cabin or cockpit for fixed- 
system absorber attachment. Where two blocks merge together, the degrees- 
of-freedom of the attachment point are eliminated and replaced by fuselage 
modal degrees-of-freedom. 

2.1 Analytical Model 

Derivation of the coupled analysis is shown in Appendix D. 

A description of the analytical model which consists of primarily 
rotor, fuselage, rotating- system absorbers and fixed-system absorbers is 
given below. 


I 


| 


! 


I 


I 

Rotor 

1 

The rotor system is represented by a modal approach which utilizes j 

the equations of motion from Reference 1. The rotor blade degrees-of- 
freedoms which can be incorporated in the analysis are: up to four blade 

elastic modes (coupled flatwise/edgewise), up to 2 blade torsional blade 
modes (first mode represents a rigid blade while the second one is an 
elastic mode) and rigid blade flapping and lead- lag motions - a total of 8 
blade modes which correspond to 24 d.o.f. (each mode has one symmetric and 
two cyclic components). The rotor/airframe coupling terms are incorporated 
in the analysis using 5 airframe modes corresponding to uncoupled fuselage 
longitudinal (x), lateral (y), vertical (z), roll (e x ) and pitch (e ) 
motions - yaw motion (e z ) is not included. y 

I 

I 

The major assumptions made in the development of the rotor system 
model are listed below: f 

I 

1. Dynamic and aerodynamic effects assume small perturbations about j 

steady initial values of the system generalized coordinates. 

2. Aerodynamic forces are developed using strip theory. 

3. Number of rotor blades must be greater than two due to the polar 

symmetry assumption made in the rotating system generalized coordi- 
nate transformations. j 



4. The following effects are not included in the rotor analysis: 

I | 

a. Forward flight aerodynamics 

b. Rotor speed d.o.f. 

j c. Variable inflow over the rotor disc 

) d. Unsteady aerodynamics 

e. Airframe yaw motion 

i ■■ ■ 4 

Fuselage 

The fuselage dynamic model is a set of linear modal equations which 
are provided in Appendix D. The computer program accepts inputs of 
system modal properties of up to 16 airframe modes. 

Rotating-System Absorbers 

l 

The bifilar analysis includes linear and non-linear inplane rotor- 
head absorbers and linear vertical absorbers. The forced response analysis 
can use up to 5 types of linear absorbers (inplane plus vertical). A 
maximum of 12 non-linear inplane pendulums can be employed in the time 
history analysis. Viscous damping of the absorbers is assumed. 

Fixed-System Absorbers 

| 

Fixed system absorbers are modeled in the analysis as a simple uni- 
directional spring-mass-damper system. The absorbers attachment point I 

must be at a defined modal vector point. Provisions for up to 5 absorbers 
are provided for in the analysis. 


2.2 Program Execution 

The rotor/bifilar coupled program starts by calculating the dynamic 
and aerodynamic rotor/ airframe matrices (assuming that rotor coupling has 
been requested) and couples them with the bifilar analysis fixed system 
dynamic matrices. Then, it expands the matrices to include the contribu- 
tions from the fixed system absorbers and the linear inplane and vertical 
bifilar pendulums. At this point, a decision is made on the type of 
solution to be calculated: forced response or time-history. If the 

forced response solution is requested, then the generalized forces are 
calculated followed by the evaluation of the forced response solution. If 
the time-history solution is required, then the program proceeds to cal- 
culate the dynamic matrices of the non-linear inplane bifilars, adds them 
to the matrices from the linear analysis, solves for the acceleration 
vector and integrates it to obtain the velocity and displacement vectors. 
The final results are harmonically analyzed (up to 10 harmonics can be 
obtained) and printed out. 



Multiple cases, in which any number of input variables are changed, 
can be easily run (a detailed discussion can be found in Section 5 - TEST 
CASES RESULTS). If the rotor blade characteristics are not changed, then 
the rotor matrices are calculated only once, stored and used at a later 
time as needed. 

Computer running time is highly dependent on the total system degrees- ! 
of-freedom used and the type of solution requested. The time-history 
solution requires considerably greater computer time than the forced 
response solution. Typical computer running times for the time history 
solution range from one to seven minutes while the forced response execution- 
usually requires from a few seconds to one-half a minute (see Section 5.5 
for more details). 

When a nonsystem error occurs during the execution of a case, the 
program attempts to continue using the best data available; if the error is 
too fundamental for the case to be meaningful, then a partial data print- 
out, followed by an error message, is given before job termination. 


i 

| 


I 



SECTION 3 



1 

I 

: 

I INPUT DESCRIPTION 

j j 

I 1 

{ 3.1 Rotor Aeroelastic Analysis Input 1 

3.1.1 Computer Listing of Input for Rotor Aeroelastic Analysis 


Symbol 

Location 

Input Item 

Units 

RHO 

1 

Air mass density Lb sec**2/ft**4 

VS 

2 

Speed of sound 

Ft/sec 

TL 

3 

Tip loss factor 

Nd 

VIP 

4 

Rotor axial velocity 

Knots 

OMEGA I 

5 

Rotor rotational speed 

RPM 

RIP 

6 

Rotor radius 

Ft 

EIP 

7 

Blade offset 

Ft 

BLADES 

8 

No. of blades - must be greater than 2 

Nd 

KBETA 

9 

Blade flapping hinge spring constant 

Lb- in/ rad 

KGAMMA 

10 

Blade lag hinge spring constant 

Lb- in/rad 

GAMOI 

11 

Blade prelag angle - lag positive 

Deg 

BETOI 

12 

Blade precone angle - up positive 

Deg 

THETAO 

13 

Blade collective pitch angle at 75% radius 

Deg 

EB 

14 

Blade Young's modulus 

Lb/in**2 

YPH2 

15 

Distance along blade axis from center 
of rotation to push rod 


PHL 

16 

Distance from blade elastic axis to 
push rod - positive toward leading edge 

In 

ZETGAM 

17 

Fraction of critical lag damping 

Nd 

ZETPIT 

18 

Fraction of critical blade pitch 
damping - based on rotor speed 

Nd 

OMGI 

19 

Reference rotor speed for defining 
percent critical lag damping in ground 
resonance studies 

RPM 

XNEMOD 

20 

Number of blade bending modes - up to 4 

Nd 

— 

21 

Open 


ZETBLD 

22-25 

Fraction of critical damping of blade 
bending modes - based on mode frequencies 

Nd 

— 

26-105 

Open 


ALPHAl 

106 

Blade pitch/lag coupling 

Nd 


j 




| DUM1(1) 

j 

| DUM1 ( 2 ) 
DUM1(3) 

| DUM1(4) 

ROTEST 

SYSDEF 


ROTDEF 

ARTIC 

PRINT 


Control Switches 

107 Printout of 30 X 30 rotor dynamic (3) and aerodynamic (2) 
matrices 

0 - no ; 1 - yes 

108 Printout of KXK (compressed) rotor matrices (dyn. + aero. ) 
0 - no ; 1 - yes 

109 Punch out of KXK (compressed) rotor matrices (dyn. + 
aero.) 

0 - no ; 1 - yes 

110 Use of rotor matrices in bifilar anlaysis 

0. do not use 

1. use new rotor matrices j 

-1. use previous rotor matrices j 

111 1. for main rotor j 

2. for tail rotor 

112 Open j 

113 System definition - ABCDEF6H. j 

A - blade bending 

B - blade rigid body pitching j 

C - blade rigid body flapping j 

D - blade rigid body lagging 
E - fixed system modes 
F - blade elastic pitching 

6 - set to 1 * 

H - set to 1 

Element = 0 to include I 

= 1 to exclude j 

1 

114 Rotor definition - X. 

X = 1. - blade hinged flatwise and edgewise j 

= 2. - blade cantilevered flatwise and edgewise j 

= 3. - blade hinged flatwise, cantilevered edgewise 
= 4. - blade cantilevered flatwise, hinged edgewise 
= 5. - gimbaled rotor 

115 Blade pitch input control - XY. j 

X - 1 - pitch bearing follows blade out-of-plane 
root slope 

= 0 - pitch bearing remains in plane of hub or j 

preconed position j 

Y = 1 - pitch bearing follows blade inplane root slope 
= 0 - pitch bearing remains in vertical plane or 
prelagged position 

116-118 Open 

119 Main printout control - X. 

X = 3. - A, basic calculations + dyn. and aero. 

integrals | 

=4. - B, A + blade frequency input i 

= 5. - C, basic calculations only 
= 6. - D, B + blade frequency output 
- 7 or greater - same as for X = 3. 



DUM2(1) 

120 

Propeller moment option in dynamic stiffness matrix 

0. include propeller moment (default value) 

1. exclude propeller moment 

— - — 

121- 

124 

Open 

LAGKII 

125 

Blade lag damper option 

0. include lag damper 

1. exclude lag damper ] 

1 

Lag damper physical characteristics 

f 

LI 

126 

| 

LI, inches 

L2 

127 

L2, inches 

L3 

128 

L3, inches 

L4 

129 

L4, inches 

L5 

130 

L5, inches 

L6 

131 

L6, inches 

CLD 

132 

Damper damping coefficient, lb-sec/in. j 

KLD 

133 

Damper stiffness coefficient, lb/in. j 

ZETELP 

134 

| 

Blade elastic pitch modal damping, nd 

CP 

135 

Effective radius to cantilevered point in torsion 
(defaults to value of blade offset in location 7) , in. j 



136- 

199 

Open 

BN 

200 

Number of elements in blade segment chart - up to 20. 
These are defined over the length of the blade only - 
the first segment is at the root of the blade-normal ly 
the last segment should be no more than 3 percent radius 

RR 

201- 

249 

Segment lengths, inches 

■k-kic 

i 

The following blade properties are input from the center 
of rotation to the blade tip at specified radial positions 

k-kk 

NCHI 

250 

Number of elements in blade chord chart 

CHORD 

251- 

349 

Radius - in : chord - in 

NATWI 

350 

Number of elements in blade aerodynamic twist chart 

ATWIST 

351- 

449 

Radius - in : twist - deg 

(twist down is positive and is zero at 75 percent 
radius) 


.7 



NSTWI 

450 

STWIST 

451- 

549 

NCGI 

550 

CG 

551- 

649 

NACI , 

650 

AC 

651- 

749 

NEAI 

750 

EA 

751- 

849 


NGJFI 

850 

RGJF 

851- 

949 

NMBI 

950 

RMB 

951- 

1049 

NIEI 

1050 

RIE 

1051- 

1149 

NIFI 

1150 

RIF 

1151- 

1249 

NFM 

1250 

RFM&FMI 

1251- 

1349 

NEM 

1350 

REM&EMI 

1351- 

1449 

NTM 

1450 

RTM&TMI 

1451- 

1549 

NGJI 

1550 

RGJ 

1551- 


1649 


8 


Number of elements in blade structural twist chart 
Radius - in : twist - deg 

(twist down is positive. The twist at 75 percent 
radius is defined by the blade geometry consistent 
with the definition of aerodynamic twist) 

Number of elements in blade center of gravity chart 
Radius - in : CG - in 

(CG is positioned relative to and positive ahead of 
elastic axis) j 

Number of elements in blade aerodynamic center chart 
Radius - in : AG - in 

(AC is positioned relative to and positive ahead of 
elastic axis) | 

Number of elements in blade elastic axis chart 
Radius - in : EA - in I 

(EA is positioned relative to and positive aft of semi-j 
chord) 

I 

. 

The following blade properties are defined for a 
specified radial segment 

Number of elements in fl exbeam torsional stiffness 
(cross-beam rotor design) chart 
GJ - lb-in**2 : segment length - in 

Number of elements in blade weight chart 
Weight - lb/ in : segment length - in 

Number of elements in blade edgewise second moment of 
area chart 

IYY - In**4 : segment length - in 

Number of elements in blade flatwise second moment of 
area chart 

IXX - In**4 : segment length - in 

1 

Number of elements in blade flatwise mass moment of j 
inertia chart (about CG) i 

IF - lb-in-sec**2/in : segment length - in j 

Number of elements in blade edgewise mass moment of 
inertia chart (about CG) 

IE - lb-in-sec**2/in : segment length - in 

Number of elements in blade torsional mass moment of 
inertia chart (about CG) ! 

IT - lb-in-sec**2/in : segment length - in 

1 

Number of elements in blade torsional stiffness - GJ 
chart 

GJ - lb-in**2 : segment length - in 



XBR Parameters 


I YA 
1 YB 
i FRR 


1650 Distance from center of rotation to outer snubber, in. 

1651 Distance from center of rotation to inner snubber, in. 

1652 Distance from center of rotation to flexbeam root, in. 

1653- Open 

1669 


Tail Rotor Control System Parameters 


PBM 1670 Weight at blade push-rod, lb. 

— 1671- Open 

1674 

Cl 1675 Damping associated with weight above, Ib-sec/in. 

1676- Open 

1679 

PBK 1680 Stiffness associated with weight above, 1 b/ in. 

— 1681- Open 

1766 

Main Rotor Control System Parameters 

PBMM 1767 Weight at blade pushrod, lb. 

— 1768 Open 

C1M 1769 Damping associated with weight above, Ib-sec/in. 

1770- Open 

1772 

PBKM 1773 Stiffness associated with the weight above, Ib/in. 

— 1774- Open 

1778 

RSB 1779 Distance from center of rotation to pushrod connection 

on swashplate, in. 

RS 1780 Distance from center of rotation to servo actuator 

connections on swashplate, in. 

— 1781- Open 

1849 

Blade Section Aerodynamic Data 

DRGDAT 1850- Drag data 

2749 

LIFDAT 2750- Lift data 

3649 

3650- Pitching moment data 
4548 


PMDAT 



3.1.2 Description of Input for Rotor Aeroelastic Analysis 

Additional information on the input parameters is provided in this 
section. 

Each input quantity is given in the following format: 

Location No., Quantity, Units. 

Important details and comments. 

1. Air Mass Density, lb sec 2 /ft^. 

If set to zero, all aerodynamic calculations are omitted. 

2. Speed of Sound, ft/sec. 

Used to compute local blade Mach number for lift, drag, and 
pitching moment calculations. 

3. Tip Loss Factor, Nondimensional . 

Provides lift and pitching moment loss in the blade tip region. 
Drag is not affected. 

Value should be equal to or greater than one minus the non- 
dimensional length of the blade tip segment. A value of one 
constitutes no loss. 

4. Rotor Axial Velocity, Knots. 

Represents climb or sideslip velocity for a main or tail rotor 
respectively. For propellers this is the forward speed of the 
aircraft. Positive velocities are in the same direction as the 
thrust. 

5. Rotor Rotational Speed, rpm. 

When calculating blade bending frequencies at low rpm's, computer 
time is greatly increased. It is suggested that rpm not be less 
than 50 cpm. 

6. Rotor Radius, ft. 

Measured from center of rotation. 

7. Blade Offset, ft. 

Distance from center of rotation to flap and lag hinge. These 
hinges are assumed to be coincident. 


10 



8. Number of Blades, Nondimensional . 

| 

Any number greater than 2 since the analysis assumes that the 
rotors have polar symmetry. The analysis will execute for N = 2 

but the results are incorrect. 1 

| 

9. Blade Flapping Hinge Spring Constant, lb in/rad . f 

If the rotor definition, location 114, stipulates a hinged 
root boundary condition, the flapping spring provides root 
flapping restraint in the calculation of the blade elastic modes. 

It is also used in the rigid-body flapping equation. 

10. Blade Lag Hinge Spring Constant, lb in/ rad. 

Same comments as above with "lagging" substituted for "flapping". 

11. Blade Prelag Angle, Degrees. 

Lag positive. 

12. Blade Precone Angle, Degrees. 

Up positive, 

13. Blade Collective Pitch, Degrees. 

Aerodynamic blade angle at 75% radius. Leading edge up positive. 

If affects blade thrust and blade bending frequencies and mode 
shapes. 

14. Blade Young's Modulus, lb/in 2 . 

Used in calculation of blade elastic modes and steady elastic 
deflections. Appears explicitly in steady elastic deflection 
calculations and therefore must have a value if the blades are | 
flexible. 

15. Distance Along Blade Axis From Center of Rotation to Blade Pushrod, 

in, ] 

Used in calculation of pitch/flap coupling for rigid-body flapping 
and for pitch/flap and pitch/lag coupling for elastic blades. 

16. Distance From Blade Elastic Axis to Pushrod, in. j 

Used in calculation of pitch/flap coupling for rigid body flapping, 
pitch/lag coupling for rigid body lagging, and pitch/flap and 
pitch/lag coupling for elastic blades. Positive toward leading 
edge. 





17. 

: 

\ 

\ 

i 

Fraction of Critical Lag Damping, Nondimensional . 

| 

Used only in the rigid-body lag equation. Based on 
uncoupled lag frequency. This value is ignored if the 
lag damper option is exercised (loc 125 set to 0). 

00 

• 

Fraction of Critical Blade Pitch Damping, Nondimensional. 

Used only in the blade pitch equation. Based on rotor 
speed. If, for example, at a given rotor speed the blade j 

torsional frequency is 5 per rev and we wish to incorporate 
10% critical damping, we would input 0.5. 

19. 

Reference Rotor Speed, RPM. 

Holds the rigid-body lag and pitch damping coefficients, 
C, constant at the value corresponding to that at the 
chosen reference rotor speed. If actual rotor speed 
variations are made, the percentage of critical damping 
will vary. If this is not desired, a zero should be 
input. 

20. 

Number of Blade Bending Modes, Nondimensional. j 

1 

l 

Up to four elastic modes can be used. If rigid-body 
modes are also being used, the program automatically 
recognizes this and will provide the correct elastic j 

modes. For example, if two bending modes are requested ] 

and rigid-body flapping is being used, the program finds 
the first three blade modes, eliminates that mode which 
corresponds to rigid-body flapping, and provides the 
remaining two modes. A zero locks out elastic modes. j 

22.-25. 

| 

Fraction of Critical Damping of Blade Bending Modes, 
Nondimensional. 

1 

Used only in the blade bending equations. Based on modal 
frequencies. 

106. 

| 

Blade Pitch-Lag Coupling, Nondimensional. j 

Defined as degrees pitch-up per degress lead for rigid- 
body lagging motion, or degrees pitch-up per degree blade 
tip inplane lead angle measured at blade root for elastic ] 

modes. Lead/pitch-up positive. j 

} 

1 

i 

! 

j 

! 12 J 

:L_j I 1 




107. 


108. 

109. 


no. 


Ill 


113. 


114 . 


Printout option for 30 X 30 rotor matrices. The dynamic 
stiffness, damping and mass matrices and the aerodynamic 
damping and stiffness matrices are printed out if switch 
is set to 1. 


Printout option for compressed rotor matrices. Same as 
above except that only those degrees- of- freedom requested 
are printed out if switch is set to 1. 

Punchout option for compressed rotor matrices. This 
option was used before rotor was coupled to bifilar 
anlaysis internally in the program. It is excercised if 
set to 1. 


: 


Coupling switch governing use of rotor matrices in the 
bifilar analysis. 

0. - do not use rotor matrices in the bifilar analysis 

1. = calculate and use rotor matrices in the bifilar analysis 
-1. = use previously calculated rotor matrices in the bifilar 

analysis 

Control Switch for Main or Tail Rotor. 


1 = Main Rotor 

2 = Tail Rotor 


System Definition. 


Exercises primary control in the program. Overrides any 
contradictory controls. There is one exception: when 

location 114 equals 5, i.e. , a gimbaled rotor. In this 
case location 114 has executive control whereby digits 1, 

3 and 4 of the system definition are ignored. Up to 8 
blade degrees- of- freedom (4 bending modes plus 2 torsional 
modes plus rigid body flapping and inplane motions) can 
be used in the rotor analysis. 

Rotor Definition. 


1 

I 


When equal to 1, 3, or 4, motion at the hinges is restrained 
by the springs in locations 9 and 10 as appropriate. 

When equal to 5, gimbaled rotor, the program automatically 
uses rigid-body flapping, and 4 elastic modes with boundary 
conditions suitably selected to correctly define the 
first 5 gimbaled rotor modes. See Reference 1 for a 
complete explanation of this. 



Rotor Definition - X. 

X = 1 - Blade Hinged Flatwise & Edgewise 
= 2 - Blade Cantilevered Flatwise & Edgewise 
= 3 - Blade Hinged Flatwise, Cantilevered Edgewise 
= 4 - Blade Cantilevered Flatwise, Hinged Edgewise 
= 5 - G infra led Rotor 

Blade Pitch Input Control. 

This simply determines whether inboard or outboard blade 
feathering is being employed. The feathering bearing is 
always at the root of the blade. 

Blade Pitch Input Control - XY. 

X = 1 - Pitch Bearing Follows Blade Out of Plane 
Root Slope 

= 0 - Pitch Bearing Remains In Plane of Hub or 
Preconed Position 

Y = 1 - Pitch Bearing Follows Blade Inplane Root 
SI ope 

= 0 - Pitch Bearing Remains In Vertical Plane or 
Prelagged Position 

Main Printout Control. 

Provides a graduated printout capability for debugging, 
etc. 

It is generally good practice to use option 5 for the 
first case of any run in order to establish that all of 
the input is correct. 

Main Printout Control - X, 

X = 3. - A, Basic Calculations + Dynamic & Aerodynami 
Integrals 

= 4. - B, A + Blade Frequency Input 
= 5. - C, Basic Calculations Only 
= 6. - D, B + Blade Frequency Output 
>_ 7. - A (Same as X = 3. ) 




Propeller Moment Option - X, 

The propeller moment contribution to the blade rigid 
torsion degree-of- freedom can be accounted for or neglected 
depending on the rotor blade design. Removal of the 
propeller moment is necessary when the rotor design 
employs blade counterweight devices. The propeller 
moment is given by: 

Mp = ipQ w here 

Ip - \ U E - I F ) dx 

The inertia, Ip, is the difference between the blade 
edgewise and flatwise mass moments of inertia integrated 
along the blade span. Rotor speed is n and blade torsional 
deflection is given by e. The moment affects only the 
dynamic stiffness matrix. 

X = 0. - Include Propeller Moment 
= 1. - Do Not Include Propeller Moment 

Blade Lag Damper Option - X. 

X = 0. - Include Lag Damper 
= 1. - Do Not use Lag Damper 

This option allows the use of a blade lag damper which 
includes all blade motions kinematics and damping and 
stiffness constants. 

A typical arrangement is shown in the sketch below. 







126.-134. Refer to the lag damper sketch above. 

135. Effective radius to cantilevered point in torsion, in. 

All elastic deformation of the blade occurs outboard of 
this point. For example, in a conventional blade this 
might be the radius to the pitch horn; in a crossbeam 
design, it may be the radius to the outboard snubber (if 
the torque tube is assumed rigid) or it may be the radius 
to the push rod connection for a torsional ly flexible 
torque tube. The user must use his own judgement in 
deciding where the cantilever point is situated. If the 
input cantilever point radius is less than the offset, it 
is set equal to the offset. 

200. Number of Elements in Blade Segment Chart. 

Up to 20 may be used. The choice of elements is important, 
particularly in relation to cross-beam type rotors where 
it is generally necessary to have smaller elements inboard 
where the spar is highly twisted. Thus, we may have, for 
example, 10 elements describing the inboard 25% of the 
blade, and 10 elements describing the remaining 75%. 

201. -220. Segment Lengths, in. 

These are defined over the length of the blade only. That 
is, their sum should equal the blade radius minus the 
offset. The tip segment should generally be of the order 
of (1-tip loss factor) times the blade radius. 

General Comments for Inputs 250 to 849 . 

All of these quantities are input as values corresponding 
to a radial position. The first value is always at the 
center of rotation and the last at the tip. Linear 
interpolation is used between adjacent values. Each 
chart can have up to 49 pairs of radius/value coordinates. 

250.-349. Blade Chord, in. 

350.-449. Blade Aerodynamic Twist, deg. 

Defined positive nose down and must be zero at 75% radius. 
This distribution is used in the calculation of aerodynamic 
effects . 





-549. Blade Structural Twist, deg. 


Defined positive nose down and will have a value at 75% 
radius consistent with the geometric properties of the 
blade. This distribution is used in the calculation of 
the blade elastic modes, steady elastic deflections, and 
blade center of gravity offset effects. 


550.-649. Blade Center of Gravity, in. 

Positioned relative to, and defined positive ahead of, 
elastic axis. 

650.-749. Blade Aerodynamic Center, in. 

Positioned relative to, and defined positive ahead of, 
elastic axis. 

750.-849. Blade Elastic Axis, In. 


Positioned relative to, and defined positive aft of, 
blade semi-chord (for a normal blade, input = - chord/4). 

General Comments for Inputs 850 to 1649 . 

All of these quantities are input as values over a segment 
length. If an offset exists, the first pair should be 
"zero, offset". The total of the segment lengths should 
equal the rotor radius. The weight and stiffness charts 
can have up to 49 value/segment length pairs. The mass 
moment of inertia charts can have up to 24 value/segment 
length pairs. 


850.-949. Flexbeam Torsional Stiffness Distribution. 

2 

Stiffness: lb-in , Segment Length: in. 

If an effective offset is input, stiffness may be set to 
zero over this distance. 


| 





950.-1049. 

1050.-1149. 

1150.-1249. 

1250.-1299. 

1350.-1399. 

1450.-1499. 

1550.-1649. 

1650. 

1651. 


Blade Weight Distribution. 

Weight: Ib/in, Segment Length: in. 

Blade Flatwise Area Moment Distribution. 

4 

Area Mom: in , Segment Length: in. 

Blade Edgewise Area Moment Distribution: 

4 

Area Mom: in , Segment Length: in. 

Blade Edgewise Mass Moment of Inertia Distribution. 

2 

Edg. Mass Mom: lb in sec /in. Segment Length: in. 

Blade Flatwise Mass Moment of Inertia Distribution. 

2 

Flat. Mass Mom: lb in sec /in. Segment Length: in. 

Blade Torsional Mass Moment of Inertia Distribution. 

2 

Tors. Mass Mom: lb in sec /in. Segment Length: in. 

Blade Torsional Stiffness Distribution. 

2 

Stiffness: lb- in , Segment Length: in. 

Zeros may be input up to the cantilever point radius 
(loc. 135). Do not include flexbeam stiffness in 
this chart. 

Distance From Center of Rotation to Outer Snubber, in. 
Distance From Center of Rotation to Inner Snubber, in. 




I 


1 

1 

j 


I 

S 


Locations 1650 and 1651 have a value other than zero 
only when cross-beam type rotors are being treated. 
Nonzero values instruct the program to apply the 
collective pitch at the outer snubber location. The 
spar, or structural, blade angle inboard of the outer 
snubber location is then made equal to the sum of the 
collective pitch (which is reduced linearly from the 
value at the outer snubber to zero at the blade root) 
and the input structural twist. The aerodynamic blade 
angle inboard of the outer snubber is made equal to 
the value at the outer snubber. 

The torque tube is assumed to be rigid in calculations 
of pitch/flap coupling. 



1652. 


Distance From Center of Rotation to FI exbeam Root, in. 


This is the radius at which torsional deformations of 
the flexbeam may be assumed to be zero. This radius 
is also used in defining the twist distribution for 
blade bending mode computations. 

Tail Rotor Control System Parameters 

1670. Weight at Blade Pushrod, lb. 

1675. Damping Associated With Weight Above, Ib-sec/in. 

1680. Series Stiffness of One Arm of Pitch Beam, Pushrod, 

and Pitch Horn for Pitch Beam Assumed Cantilevered 
at its Center, lb/ in. 

Main Rotor Control System Parameters 

1767. Weight at Blade Pushrod, lb. 

1769. Damping Associated With Weight Above, Ib-sec/in. 

1773. Stiffness Associated With Pushrod, Swashplate Connecti 

Pitch Horn, etc..., lb/in. 

1779. Distance From Center of Rotation to Pushrod Connection 
on Swash Plate, in. 

1780. Distance From Center of Rotation to Servo Actuator 
Connection on Swash Plate, in. 

Blade Section Aerodynamic Data 


Drag, lift, and pitching moment data are each input in sets of CD, CL 
or CM pairs corresponding to up to 12 Mach numbers. 

Any or all of the data charts can be omitted. 

In each chart there are up to 35 pairs of angle-of-attack and CD, CL, 
or CM data corresponding to each Mach number. 

A total of 75 locations are allocated to each Mach number. 

Mach numbers and angles-of-attack are input in ascending order. 

If charts are employed use no less than 2 Mach numbers and 5 pairs 
of data per chart. 



The input angle-of-attack range should exceed the expected range for 
the condition being analyzed. 

Unsymmetric airfoil data can be used. 

Only one set of airfoil data can be used. 

Example of input format: 


1850 

1851 

1852-1924 

1925 

1926 

1927-1999 

Etc. 


Number of pairs for first Mach number 

First Mach number 

Angle-of-attack: CD pairs 

Number of pairs for second Mach number 

Second Mach number 

Angle-of-attack: CD pairs 


1850.-2749. Drag Data. 

2750.-3649. Lift Data. 

3650.-4549. Pitching Moment Data. 


! 

! 



3.2 Bifilar Analysis Input 
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3.2.1 Computer Listing of Input for Bifilar Analysis 


x 

x 

x 

X 
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xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
COMPUTER LISTING OF INPUT TO THE BIFILAR ANALYSIS 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


SYMBOL 

LOCATION 

DESCRIPTION 



UNITS 

NINBF 

1 

KINDS OF INPLANE BIFILARS 

= < 

5 


NAP 

2 

NUMBER OF AIRCRAFT STATION POINTS 
TO EVALUATE RESPONSE 

=< 

4 


WF 

3 

FORCING FREQUENCY 



CYCLES/REV 

NIMP 

4 

DIMENSION OF IMPEDANCE MATRIX 
(LOC 4 = LOC 9 X 2 ) 





5 

OPEN 




XNARM 

6 

TOTAL NUMBER OF INPLANE BIFILARS 

=< 

10 


OMEGAR 

7 

ROTOR SPEED 



RPM 

NFABS 

8 

NUMBER OF FIXED SYSTEM ABSORBERS 

=< 

5 


NF 

9 

NUMBER OF FIXED SYSTEM MODES 

= < 

16 


IWRITE 

10 

=1. GET EXTENSIVE PRINTOUT 





11 

OPEN 




NVBIF 

12 

KINDS OF VERTICAL BIFILARS 

=< 

5 


XNARM 

13 

TOTAL NUMBER OF VERTICAL BIFILARS 

=< 

10 


RHM 

14 

=1. CALCULATE ROTOR HEAD MOTION ONLY 


ICHECK 

15 

=1. WRITE M,C,K MATRICES AND IMPEDANCE 



16 

OPEN 




ROTPNT 

17 

=0. NO PRINTOUT OF ROTOR MATRICES 
=1. GET PRINTOUT OF ROTOR MATRICES 



NLBFSW 

18 

=0. NO NON-LINEAR INPLANE BIFILARS 
=1. USE NON-LINEAR INPLANE BIFILARS 




19 

OPEN 





FIXED SYSTEM ABSORBER DATA 


FABSM 

FABSWN 

FABSD 

20-29 

30-39 

40-49 

FIXED SYSTEM ABSORBER MASS ARRAY 
FIXED SYSTEM ABSORBER FREQUENCY ARRAY 
FIXED SYSTEM ABSORBER DAMPING ARRAY 

BUGS 

CPM 

FIXED SYSTEM MODAL DATA 

XMG 

XWN 

XDP 

50-69 

70-89 

90-109 

GENERALIZED MODAL MASS ARRAY 
GENERALIZED MODAL FREQUENCY ARRAY 
GENERALIZED MODAL DAMPING ARRAY 

BUGS 

CPM 

FORCES AND MOMENTS DATA 

FHC 

FHS 

FTC 

110-119 

120-129 

130-139 

COSINE COMPONENT OF HUB FORCES & MOMENTS 
SINE COMPONENT OF HUB FORCES & MOMENTS 
COSINE COMPONENT OF TAIL ROTOR LOADS 

LB a IN-LB 
?? 

11 



FTS 140-149 SINE COMPONENT OF TAIL ROTOR LOADS " 
FFC 150-159 COSINE COMPONENT OF HORIZONTAL TAIL LOADS " 
FFS 160-169 SINE COMPONENT OF HORIZONTAL TAIL LOADS " 
FEC 170-179 COSINE COMPONENT OF LOADS OF ANOTHER POINT " 
FES 180-189 SINE COMPONENT OF LOADS OF ANOTHER POINT " 


INPLANE BIFILAR PENDULUMS DATA 


BIFM 190-199 INPLANE BIFILAR MASS ARRAY BUGS 

BARM 200-209 BIFILAR ARM FROM CENTER OF ROTATION IN 

BIFWN 210-219 INPLANE BIFILAR TUNING FREQUENCY ARRAY CYCLES/REV 

BIFDAM 220-229 INPLANE BIFILAR DAMPING ARRAY 


VERTICAL BIFILAR PENDULUMS DATA 


BIFM 230-239 VERTICAL BIFILAR MASS ARRAY BUGS 

BIFARM 240-249 BIFILAR ARM FROM CENTER OF ROTATION IN 

BIFWN 250-259 VERTICAL BIFILAR TUNING FREQUENCY ARRAY CYCLES/REV 

BIFDAM 260-269 VERTICAL BIFILAR DAMPING ARRAY 


270-449 OPEN 

XPHI 450-549 TRANSFER MATRIX TO MODAL COORDINATES 

FOR MAIN ROTOR IMPEDANCE (6 X NF) 


MODE SHAPES OF 2 AIRCRAFT STATIONS WHERE LOADS ARE APPLIED 


XPF 550-649 FIRST AIRCRAFT STATION (6 X NF) 

XPE 650-749 SECOND AIRCRAFT STATION (6 X NF) 


MODE SHAPES OF FIXED SYSTEM ABSORBERS & TAIL ROTOR HUB 


XPHABS 750-849 FIXED SYSTEM ABSORBER COUPLING MATRIX (NFABS X NF) 

(20 SPACES RESERVED FOR EACH ABSORBER) 

XPT 850-949 TAIL ROTOR HUB TRANSFER MATRIX 

______ 950-999 OPEN 


MODE SHAPES OF AIRCRAFT STATION POINTS 


XPHAP 1000-1399 COUPLING MATRIX OF AIRCRAFT STATION POINTS 

(6 X NF X NAP) 

1000-1099 FOR FIRST AIRCRAFT STATION 

1100-1199 FOR SECOND AIRCRAFT STATION 

1200-1299 FOR THIRD AIRCRAFT STATION 

1300-1399 FOR FOURTH AIRCRAFT STATION 

!400-1489 OPEN 


PRINTOUT OPTIONS 


FSPNT 

FSRPNT 

FSAPNT 

IBPPNT 

VBPPNT 

IBQPNT 

VBQPNT 

I GAMMA 

ICUSE 


=0. PRINTOUT IS NOT WANTED 
=1. PRINTOUT IS WANTED 

1490 A = FIXED SYSTEM MATRICES 

1491 B = A + ROTOR MATRICES 

1492 C = B + FIXED SYSTEM ABSORBERS 

1493 D = C + INPLANE BIFILAR PENDULUMS 

1494 E = D + VERTICAL BIFILAR PENDULUMS 

1495 INPLANE BIFILAR MATRICES ONLY (9 X 9) 

1496 VERTICAL BIFILAR MATRICES ONLY (9 X 9) 

1497 PRINTOUT OF "GAMMAS" (DEGREES-OF-FREEDOM) 

1498 =0. DO NOT USE INITIAL CONDITIONS OPTION 



=1. USE INITIAL CONDITIONS FROM PREVIOUS CASE 
(SEE LOCATIONS 1720-1939) 

U99 OPEN 

XXXXXXXXXXX END OF INPUT FOR LINEAR BIFILAR ANALYSIS XXXXXXXXXXXXXXXXX 
xxxxxxxx INPUT FOR NON-LINEAR INPLANE BIFILARS ANALYSIS xxxxxxxxxxxxx 


INPLANE BIFILAR PENDULUMS CHARACTERISTICS 


BM 1500-1519 BIFILAR WEIGHTS LBS 

DP 1520-1539 BIFILAR DAMPING 

WP 1540-1559 BIFILAR NATURAL FREQUENCIES CYCLES/REV 

BR 1560-1579 ARM DISTANCES FROM CENTER OF ROTATION IN 

PSIR 1580-1599 RELATIVE AZIMUTH LOCATIONS DEGREES 


HUB FORCES X MOMENTS HARMONICS INPUT 


FORMAT IS " COSCKXPSI) + SIN<K*PSl) " 

WHERE K=1 THROUGH A MAXIMUM VALUE OF 10 

XX XXX XXXXXXXXXXXXX XXXXXXXXX XXXXXXXXXXXXX XXXXXXXXXXXXXXXXXX 


FXX 

1600-1619 

LONGITUDINAL HUB FORCE 

(X) 

LBS 

FYY 

1620-1639 

LATERAL 

HUB FORCE 

(Y) 

LBS 

FZZ 

1640-1659 

VERTICAL 

HUB FORCE 

(Z) 

LBS 

XXM 

1660-1679 

ROLL 

HUB MOMENT 

(THETAX) 

IN-LB 

YYN 

1680-1699 

PITCH 

HUB MOMENT 

(THETAY) 

IN-LB 

ZZM 

1700-1719 

YAW 

HUB MOMENT 

(THETAZ) 

IN-LB 

GAMMA 

1720-1739 

INITIAL 

BIFILAR ANGULAR 

DISPLACEMENT 

RADIANS 

GAMMAD 

1740-1759 

INITIAL 

BIFILAR ANGULAR 

VELOCITY 

RAD/SEC 

PS IP 

1760 

FORCE RAMP FACTOR( 3. IS 

RECOMMENDED) 


DPSI 

1761 

DELTA AZIMUTH ANGLE 


DEGREES 

FPST 

1762 

MAXIMUM 

AZIMUTH ANGLE 


DEGREES 


NBIF 1763 NUMBER OF BIFILARS ( MAX=12) 

NHFH 1764 NO. OF INPUT HUB FORCE HARMONICS (MAX=10) 

XNH 1765 NO. OF OUTPUT STATE VARIABLES HARMONICS (MAX=1 0 ) 

PSI 1766 INTEGRATION VARIABLE (NOT AN INPUT ITEM) 

NFP 1767 NUMBER OF POINTS ON AIRCRAFT WHERE FORCE 

IS APPLIED BESIDES ROTOR HEAD (=< 2) 

(MODE SHAPES ARE INPUT IN LOC. 550-749) 


HUB & STATE VARIABLES INITIAL VALUES 


HUBD 1768-1773 HUB INITIAL DISPLACEMENTS 

HUBV 1774-1779 HUB INITIAL VELOCITIES 

STDP 1780-1859 STATE VARIABLES INITIAL DISPLACEMENTS 

TVL 1860-1939 STATE VARIABLES INITIAL VELOCITIES 


FORCES X MOMENTS HARMONICS INPUT OF 2 ADDITIONAL AIRCRAFT POINTS 


FORMAT IS " COS(K*PSI) + SIN(K*PSI) " 
WHERE K=1 THROUGH A MAXIMUM VALUE OF 10 
(MODE SHAPES ARE INPUT IN LOC. 550-749) 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


P1FX 

1940-1959 

POINT 1 

LONGITUDNAL 

(X) 

LBS 

PI FY 

1960-1979 

w 

LATERAL 

(Y) 

LBS 

P1FZ 

1980-1999 

it 

VERTICAL 

(Z) 

LBS 

P1MX 

2000-2019 

ii 

ROLL 

(THETAX) 

IN-LB I 

P1MY 

2020-2039 

ii 

PITCH 

(THETAY) 

IN-LB 


24 



P1M2 

2040-2059 

If 

YAW 

CTHETAZ) 

IN-LB 

P2FX 

2060-2079 

POINT 2 

LONGITUDINAL 

(X) 

LBS 

P2FY 

2080-2099 

II 

LATERAL 

CY) 

LBS 

P2FZ 

2100-2119 

If 

VERTICAL 

(Z) 

LBS 

P2MX 

2120-2139 

If 

ROLL 

(THETAX) 

IN-LB 

P2MY 

21 AO-2159 

If 

PITCH 

(THETAY) 

IN-LB 

P2MZ 

2160-2179 

2180-2199 

If 

OPEN 

YAW 

CTHETAZ) 

IN-LB 


SUCCESSIVE CASES CONTROL SWITCH 


CODE 2200 =1. COMPLETES LAST CASE 

=0. GOES BACK TO ROTOR ANALYSIS FOR NEXT CASE 


END OF INPUT LOCATIONS xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 





3.2.2 Description of Input for Bi filar Analysis 

| 

Additional information on the input parameters is provided in this 
| section. 

Each input quantity is given in the following format: 

Location No., Quantity, Units 
Important Details and Comments 


Up to 5 kinds of inplane bifilars can be used in the 
analysis if no vertical bifilars are present. The combined 
number of kinds of inplane and vertical bifilars cannot be 
greater than 5. This corresponds to a maximum number of 
15 degrees - of- freedom (including collective and 2 cyclic 
motions). Ex. loc 1=2 if one set of 3P and one set of 5P 
inplane bifilars are employed. 

2. Aircraft Station Points, Nondim. 

Due to computer storage requirements, a maximum of 4 points 
can be used. j 

3. Forcing Frequency, Cycles/Rev. j 

| 

Frequency at which the bifilars will respond divided by 
the rotor speed. j 

1 

4. Impedance Matrix Dimension, Nondim. f 

I 

This location is 2 times the number of fixed system modes 
loaded in location 9. 

6. Number of Inplane Bifilars, Nondim. 

The actual number of inplane bifilars present in each kind | 
as defined in loc. 1. It should be noted that, if more 
than one kind of inplane bifilars is used, then the number 
of bifilars in each kind has to be the same. The maximum 
number allowed is 10. 

7. Rotor Speed, RPM. 

Used to non-dimensional ize frequency units. 


Linear Analysis Input (Loc 1-1498) 


Kinds of Inplane Bifilars, Nondimensional . 




Number of Fixed System Absorbers, Nondim. 

The total number of fixed system absorbers is limited to 5. 
Number of Fixed System Modes, Nondim. 

| 

The maximum number of fixed system modes is 16 which allows 
for 6 rigid body modes (if desired) and 10 flexible modes. 

Rigid airframe modes should not be used in the non-linear 
bifilar analysis. 

Program will not execute if number of modes is zero. 

Extensive Printout Option, Nondim. j 

Set this location to 1. to obtain additional printout for 
debugging use. j 

Kinds of Vertical Bifilars, Nondim. 

Up to 5 kinds of vertical bifilars can be used if no 

inplane bifilars are present. Otherwise, the total of j 

the two kinds is 5. j 

j 

Number of Vertical Bifilars, Nondim. 

Actual number of vertical bifilars present in each kind 
as defined in loc. 12. Same number of bifilars is 
assumed for each kind. Maximum number is 10. 

) 

Rotor Head Motion Switch, Nondim. 

If this control is set to 1. then only the rotor head motion 
will be calculated. ] 

Additional Printout Option, Nondim. 

If set to 1., stiffness, damping and mass matrices are 
printed out for debugging purposes. 

Rotor Matrices Printout Option, Nondim. 

If set to 1., the input matrices to the bifilar anlaysis 
from the rotor aeroelastic program are printed out. Order 
is stiffness, damping and mass. The stiffness and damping 
matrices include aerodynamics if air density is non-zero. j 

Linear/Non-Li near Analysis Option, Nondim. 

If set to zero, then linearized inplane bifilar equations 

GH 27 




of motion are used and the forced response of all components 
is calculated. Sample runs using 16 and 28 degrees-of- 
freedom required 3 seconds and 25 seconds respectively. 


If set to 1., then the full non-linear inplane bifilar 
equations of motion are used and a time history solution 
is calculated. Computer time for the non-linear option 
is highly dependent on the number of degrees-of- freedom 
used. A sample run with 17 d.o.f. took 1 minute and 8 
seconds requiring 13 rotor revolutions for convergence of 
the bifilar motions. A sample run with 29 d.o.f. took 7 
minutes and 7 seconds of computer time and required 16 
rotor revolutions for convergence. 

2 

20.-29. Fixed System Absorber Masses, Lb-sec /in. or Bugs 

Absorber weight in pounds divided by 386.4. Although 
10 locations are provided, only 5 can be used (see loc. 


30.-39. Fixed System Absorber Frequencies, Cpm j 

40.-49. Fixed System Absorber Damping, Nondim. I 

Damping associated with the absorber i.e., .01 correspond- 
ing to 1 % critical. j 

' 

2 

50.-69. Fixed System Generalized Masses, Lb-sec /in or Bugs 

Maximum number of masses is actually 16 (see loc. 9) 
although 20 computer locations have been allowed. 

: 

70.-89. Fixed System Generalized Frequencies, Cpm 

90.-109. Fixed System Generalized Damping, Nondim. 

Forces and Moments Data -General Comments, Lb & In-lb 

110.-189. The Forces and Moments are Defined in the Fixed System. 

The order of the forces and moments input data is: x, y, 

z, e^, e , e_, which correspond to longitudinal, lateral, 
vertical - ; roll, pitch and yaw motions respectively. The 
coordinate system used is a right-handed (system with 
the x-axis defined as positive aft and y-axis as positive 
out of the right wing). 

| 

Thus only six locations are needed to define an input force 
or moment although 10 computer locations have been allocated. I 



190. 


200 . 

210 . 


220 . 

230. 

240. 

250 . 


Examples of input of hub forces (loc. 110-129): 

a) For a 4P inplane bifilar, to load a pure 5P hub force, 
then F x (cosine) = -F y (sine) (loc. 110 = - loc. 121) 
and Fy (cosine) = Fx \sine) (loc. Ill = loc. 120). 

b) To load a pure 3P hub force, then F x (cosine) = F y 
(sine) (loc. 110 = loc. 121) and F v (cosine) = -Fx 
(sine) (loc. Ill = - loc. 120). y 

i 

Additional forces/moments inputs can be specified for the 
tail rotor (loc. 130-149), the horizontal tail (150-169) j 

and one other arbitrary point (loc. 170-189). j 

2 

199. Inplane Bifilar Masses, Lb-sec /in. or Bugs 

Bifilar weights divided by 386.4 are loaded for each kind I 
of bifilar used (see loc. 1). The limit on the number 
of bifilar kinds which can be used is 5 although 10 computer 
locations are available. j 

I 

209. Inplane Bifilar Arm From Center of Rotation, In. 

219. Inplane Bifilar Tuning Frequency, Cycles/Rev J 

The linear inplane bifilar tuning frequency is defined by 

}' 

F 2 = R/r where R = bifilar arm (loc. 200-209) 

r = bifilar distance from hinge to j 

bifilar center of mass ! 

Example: Given R = 18.22 in, r = 2.02444, then F = 3.0 ! 

per rev j 

229. Inplane Bifilar Damping, Nondim. j 

2 

239. Vertical Bifilar Masses, Lb-sec /in. or Bugs 

Same comment as for inplane bifilar (see loc. 12). 

| 

249. Vertical Bifilar Arm From Center of Rotation, In. 

| 

259. Vertical Bifilar Tuning Frequency, Cycles/Rev 

The linear vertical bifilar tuning frequency is defined by 
r 2 = ( R+ r ) /r where R = bifilar arm (loc. 240-249) 

r = bifilar distance from hinge to j 

bifilar center of mass 



Example: Given R = 18.50 in, r = 1.23333 in, then F = 4.0 

per rev 

260.-269. Vertical Bifilar Damping, Nondim. 

450.-549. Transfer Matrix to Modal Coordinates for Main Rotor 
Impedance, In/ in & Rad/ in. 

The main rotor hub transfer matrix has the dimensions 
6 X NF (where NF is the number of fixed system modes 
defined in loc. 9). Since the maximum value of NF is 
16, then 96 computer locations are necessary to define the 
largest transfer matrix. 

The first mode is loaded into locations 450-455, the 
second mode into locations 456-461, etc. .. until the 
last or 16th mode into locations 540-545. 

For each mode, the order of input is x, y, z, e x , e y , e z . 
The units of the linear motions are in/in while rotations 
are in rad/in. The user must be careful to load the 
mode shapes corresponding to the fixed system modes 
masses, frequencies and damping values from locations 
50-109. 

550.-649. Mode Shapes for First Aircraft Station Where Loads Are 
Applied, In/in. & Rad/in 

650.-749. Mode Shapes for Second Aircraft Station Where Loads Are 
Applied, In/in & Rad/in 

Same comments as for the main rotor hub transfer matrix 
above apply for these inputs. 

750.-849. Fixed System Absorber Coupling Matrix, In/in & Rad/in 

The fixed absorber coupling matrix has the dimension 
NFABS X NF (where NFABS is the number of fixed system 
absorbers from location 8). The modal response in one 
direction (vertical, lateral or longitudinal) for each 
fixed system mode (as defined in location 9) is loaded 
in locations 750-769 for the first absorber, in locations 
770-789 for the second, and so on. 

850.-949. Tail Rotor Hub Transfer Matrix, In/in & Rad/in 

Same comments as above for main rotor hub transfer 
matrix (loc. 450-549). 



1000.-1399 



Mode Shapes of Aircraft Station Points, In/in & Rad/in 

The response of 4 aircraft stations can be analyzed 
according to the mode shapes input in locations 
1000-1099 for the first station, 1100-1199 for the 
second, 1200-1299 for the third, and 1300-1399 for 
the fourth station. 

The mode shapes are loaded in the same manner as 
described above for the main rotor hub (loc. 450- 
549). 

1490.-1494. Printout Options, Nondim. 

If any printout option switch is set to 0., then the 
printout of the corresponding matrices is suppressed. 

If printout of the matrices is desired, then the proper 
switch should be loaded as 1. 

The build-up of the matrices is as follows: 

(loc. 1490) 1. Fixed system matrices (16 X 16 max) 

(loc. 1491) 2. Add rotor matrices (24 X 24 max -* 

40 X 40 max total) 

(loc. 1492) 3. Add fixed system absorbers (5 X 16 

max -* 45 X 45 max. total) 

(loc. 1493) 4. Add linear inplane bifilars 

(loc. 1494) 5. Add vertical bifilars (combined 

with inplane bifilars, the matrices 
are 15 X 15 max •* 60 X 60 max. total) 

Thus, the maximum number of degrees-of- freedom which 
can be handled at the present time with the linear 
rotor/bifilar coupled program is 60. 

Location 1493 also controls the printout of the "Final 
Combined Mass Matrix" for the non-linear analysis option. 

1495.-1496. Bifilar Pendulums Printout Switches, Nondim. 

The individual contributions to the coupled system 
matrices from the inplane and vertical bifilars can 
be obtained through the use of the printout options 
as defined in locations 1495 and 1496 respectively. 

1497. The results of the forced response analysis for all 

the system degrees-of- freedom can be printed out 
through this switch. 
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1498. Initial Conditions Option, Nondim. 

This option applies only when the non-linear analysis 
is requested (loc. 18 = 1. ). If set to 1.0, then 
the initial conditions to be used should be loaded in 
locations 1720-1939. Input to these locations is 
printed out at the end of a non-linear analysis run. 

Non-Linear Analysis Input (Loc. 1500-2179) 

1500.-1519. Inplane Bifilar Pendulum Weights, Lbs. 

The weight of each bifilar pendulum in pounds is 
loaded according to the number of pendulums from loc. 
1763. The maximum number of pendulums allowed is 
12 . 


1520.-1539. Inplane Bifilar Pendulum Damping, Nondim. 

1540.-1559. Inplane Bifilar Pendulum Frequencies, Cycles/ Rev 


1560.-1579. 

1580.-1599. 


1600.-1719. 


Same comments as provided for the linear inplane 
bi filers (loc. 210-219) apply here. 

Inplane Bifilar Pendulum Arms From Center of Rotation, 
In. 

Relative Azimuth Locations of Inplane Bifilar 
Pendulums, Deg 

Example: If four inplane pendulums are analyzed, 

then locations 1580-1583 are respectively 
0., 90., 180., 270. 

Hub Forces and Moments Input, Lbs and in. -lb 

For the non-linear analysis, the fixed system hub 
forces and moments are input in harmonics format. 

Up to 10 harmonics can be used. Cosine and sine 
components of each harmonic are input alternately. 

Ex. Loc 1600 - Longitudinal motion - cosine 
component of first harmonic 
Loc 1601 - Longitudinal motion - sine component 
of first harmonic 

Loc 1602 - Longitudinal motion - cosine 
component of second harmonic 

Etc. . . 


Location 1764 is used in conjunction with the input 
loads. 



1720,- 

1740.' 

1760. 

1761. 

1762. 

1763. 


1739. Initial Bifilar Angular Displacements, Radians 

These locations are used if the initial conditions 
switch (loc. 1498) is 1. These values are printed 
out at the end of a converged time history solution 
for each inplane bifilar. 

1759. Initial Bifilar Angular Velocity, Rad/Sec 

Same comments as above for the initial displacements. 
Force Ramp Factor, Nondim. 

The hub loads are imposed on the rotor/bifilar/fixed 
system coupled system as a ramp input dependent on 
the maximum azimuth angle (loc. 1762) for the time 
history solution and the ramp factor. 

Example: If loc. 1762 = 4320 degrees (or 12 

revolutions) and loc. 1760 = 3.0, then the hub 
loads (loc 1600-1719) are applied linearly in the 
azimuth interval from zero to 4320/3 (which equals 
1440 degrees or 4 revolutions). 

Azimuthal Increment For Time History Solution, Degrees 

A value of 2 degrees is recommended. However, if the 
time history does not converge, lower values can be 
tried to eliminate what may possibly be a numerical 
instability. It should be noted that the computer 
execution time is directly proportional to the size 
of this input quantity. 

Maximum Azimuth Angle, Degrees 

It is recommended that values corresponding to 10 to 
20 rotor revolutions (3600 and 7200 degrees respect- 
ively) be used in this location. As for loc. 1761 
above, the computer execution time is also dependent 
on this input. If the time history does not meet 
the convergence criteria, then it will terminate 
when the integration azimuthal angle equals the 
input value in loc. 1762. 

Number of Inplane Bifilar Pendulums, Nondim. 

A maximum number of 12 bifilars can be used. 



1764. 


1765. 


1767. 


1768. 


1774- 


1780. 


Number of Input Hub Force Harmonics, Nondim. 

A maximum number of 10 harmonics, corresponding to 
the input loads in loc. 1600-1719, can be used. 

Number of Output State Variables Harmonics, Nondim. 

This location governs the number of harmonics analyzed 
and printed out for the following variables: 

1. Bi filar pendulum response, degrees 

2. Hub response (x, y, z, e x , By, e z ), g's 

3. Aircraft stations response (x, y, z), g's 

Up to 10 harmonics may be requested. 

Number of Aircraft Additional Points Where Loads Are 
Applied, Nondim. 


This number is equal or less than 2. If non- zero, 

then load appropriate mode shapes in locations 550-749. j 


-1773. Hub Initial Displacements, In or Rad 

At the completion of the time history solution, the 
hub displacements are printed out to be used for 
successive cases if desired. The output yields, 
in order, the longitudinal, lateral, vertical, roll, 
pitch and yaw motions. 

1779. Hub Initial Velocities, In/sec or Rad/sec 

Same comments as above for the hub initial displace- 
ments. 

-1859. State Variables Initial Displacements, In or Rad 

The initial displacements of all the degrees-of- 
freedom (except the non-linear inplane bifilar 
pendulums) are printed out for use in successive 
cases. The order of printout is: 


;■) 

i 

| 


| 


| 

I 


1. Fixed system 

2. Rotor (if required) 
3Fixed system absorbers 

4. Linear inplane bifilars 

5. Linear vertical bifilars 



State Variables Initial Velocities, In/sec or rad/sec. 

Same comments as made above for the initial displace- 
ments. 

Aircraft Additional Points Forces and Moments Input, 
Lbs and in-lb. 

Refer to comments on hub forces (locations 1600-1719). 

2200. Successive Cases Control Switch, Nondim. 

If successive runs are to be made, then location 
2200 is set to zero. The program then goes back 
to the rotor aeroelastic analysis and starts execution 
of the next case. The last bifilar analysis case 
must have a 1. in location 2200 for proper termination 
of the computer run. 


The maximum number of degrees- of- freedom which can handled by the 
linear and non-linear analyses are respectively 60 and 72. A breakdown of 
the individual component maximum d.o.f. is presented in the chart below. 


Component 

Description 

Maximum No 

. of D.O.F. 

Input 

Linear 

Non-Linear 

Locations ) 

1 . 

Linear Bifilars 
(inplane + vertical) 

15 

15 

1 & 12 

2. 

Fixed System Absorbers 

5 

5 

8 

3. 

Fixed System Modes 

16 

16 

9 

4. 

Rotor Blade 

24 

24 

- 

5. 

Non-linear Inplane 
Bifilars 

— « 

12 

1763 

Total ► 

60 

~7T 



3.3 Input Data Format 


All data is input via cards or card-like images with the following 
(loader) format. Column 2 represents the number of data values on the 
card, columns 3-6 give the location numbers of the first data values, 
successive data values are loaded into successive locations, and columns 
7-66 contain the data values in fields of 12, the default format being 
5E12.4. A minus sign in column 1 indicates the end of a case. Subsequent 
cases are loaded immediately after this card in the same format. See 
subroutine L0ADIT for a more detailed description of data cards. 


1860.-1939. 


1940.-2179. 




SECTION 4 


OUTPUT DESCRIPTIONS 


if 



Rotor Aeroelastic Analysis Output 


If the rotor coupling option is exercised (location 110 is not 
zero), then a listing of the input rotor data is first printed out. A 
j sample page of this output can be seen in Figure 2. The input cards are 
listed out as read by the computer. For successive cases, only the new 
input items will be printed out to reflect the changes made between runs. 
This information is provided for all printout options. 

The output formats which follow are obtained for all printout options, 
as provided in location 119. Additional printout of blade bending freq- 
uency calculations and dynamic and aerodynamic integrals evaluations can 
be obtained if the printout switch is not equal to 5. These formats are 
not shown here for sake of brevity - they are mainly used for debugging 
purposes. 


Input and calculated dynamic and aerodynamic characteristics of the 
rotor blade and fixed system are presented in Figures 3a through 3h. 
Additional explanations are provided in Table 1 below. Some differences 
may appear in the output formats from the input data due to corrections 
applied by the computer program to eliminate possible inconsistencies in 
the input control options. 

The abbreviations used in the output formats are listed and discussed 
in the table below. It should be noted that some input quantities have 
been preset within the computer program since some program capabilities 
were not needed for the coupling of the rotor and bifilar systems. 


Output 
Symbol 

PHIXPH 
PHIZPH 

.... . Lag Damper Quantities ..... . 

PHELD Edgewise Bending Mode at Lag Damper 

PHEPLD Edgewise Bending Mode Slope at Lag 

Damper 

PHFLD Flatwise Bending Mode at Lag Damper 

PHFPLD Flatwise Bending Mode Slope at Lag 

Damper 

QEOLD Edgewise Steady Deflection at Lag 
Damper 


TABLE 1. Rotor Analysis Output Description 

Quantity Description Input Present Figure 

Location Value Number 

Blade Bending Mode Coupling Factors 3a 

Blade Bending Mode Coupling Factors 




Output 

Quantity Description Input 

Present 

Figure 

Symbol 

Location 

Value 

Number 

QEOPLD 

Edgewise Steady Slope at Lag Damper 


- 

3a j 

QFOLD 

Flatwise Steady Deflection at Lag 

- 

- 

i 

j 


Damper 




QFOPLD 

Flatwise Steady Slope at Lag Damper 

- 

- 


PHLD 

Torsional Mode Shape at Lag Damper 

- 

- 


THTLD 

Blade Twist at Lag Damper 

- 


i; 

PHOS 

Torsional Mode Shape at Outboard 
Snubber 

' 


l 

| 

XNAMOD 

Number of Fixed System Modes 

_ 

5. 

3b 

V F 

Forward Flight Speed 

- 

0. 







ROTEST 

Rotor Definition 

Ill 

— 


FTEST 

Flight Definition 

- 

1. 

j 

SYSDEF 

System Definition 

113 

- 

f 

\ 

ROTDEF 

Rotor Definition 

114 

- 

•| 

ARTIC 

Blade Pitch Input Control 

115 

- 

: j 

PHASE 

Phasing Matrix Printout Control 

- 

0. 

i 

VECT 

Eigenvector Printout Control 

- 

0. 

] 

TRMASC 

Tail Rotor Main Mass Control 

- 

1. 


SUMASC 

Tail Rotor Subsidiary Mass Control 

- 

111. 

j 

TSERVC 

Tail Rotor Servo Control 

- 

1. 

■j 

MRMASC 

Main Rotor Mass Control 

- 

1. 

\ 

MSERVC 

Main Rotor Servo Control 

- 

111. 

'I 

CIR 

Circulatory Unsteady Aerodynamics 
Control 

- 

1. 

i 

CIRN 

Noncirculatory Unsteady Aerodynamics 
Control 

“ 

1. 

i 

LAGKII 

Lag Damper Control 

125 

- 

i 

ZETBLD 

Fraction of Critical Damping of Blade 
Bending Modes 

22-25 








ZETG 

Fraction of Critical Damping of 
Fixed System Modes 

- 

0. 

1 

MG 

Generalized Mass of Fixed System 
Modes 

— 

0. 


OMF 

Frequency of Fixed System Modes 

- 

0. 

i 

PHY 

Lateral Fixed System Mode Shape 
(Second Mode Only Equals 1.0) 

*" 

0. & 1. 


PHX 

Longitudinal Fixed System Mode Shape 
(First Mode Only Equals 1.0) 

— 

0. & 1. 

i 

PHZ 

Vertical Fixed System Mode Shape 
(Third Mode Only Equals 1.0) 

— 

0. & 1. 

; 



Output 

Quantity Description 

Input 

Present 

Symbol 

Location 

Value 

PHTY 

Pitch Fixed System Mode Shape 


0. & 1. 

PHTX 

(Fifth Mode Only Equals 1.0) 
Roll Fixed System Mode Shape 


0. & 1. 

R 

(Fourth Mode Only Equals 1.0) 
Radius of the Blade Element Mid- 

201-220 


AC 

Points From Center of Rotation 
Blade Aerodynamic Center 

650-749 


CG 

Blade Center of Gravity 

550-649 

- 

EA 

Blade Elastic Axis 

750-849 

- 

QEO 

Blade Edgewise Steady Deflection 


- 

QFO 

Blade Flatwise Steady Deflection 


- 

QEOP 

Derivative of QEO with respect 


- 

QFOP 

to R 

Derivative of QFO with respect 



DT 

to R 

Elemental Thrust 



DH 

Elemental Drag 


- 

DM 

Elemental Pitching Moment 


- 

D 

Derivative with Respect to 


- 

dTutT 

D 

Tangential Velocity 
Derivative with Respect to 



dTupT 

D 

Vertical Velocity 
Derivative with Respect to 



W) 

CL 

Angle-of-Attack 
Coefficient of Lift 

2750-3649 


CD 

Coefficient of Drag 

1850-2749 

- 

CM 

Coefficient of Pitching Moment 

3650-4548 

- 

D 

Derivative with Respect to 


- 

DA 

D 

Angle-of-Attack 

Derivative with Respect to Mach 



DM 

0 STRUC- 

Number 

Structural Twist 

450-549 


TURAL 
0 AERO- 

Aerodynamic Twist 

350-449 


DYNAMIC 

Up 

Hover Inflow Velocity 



Ut 

Tangential Velocity 


- 

U 

Total Velocity 


- 

0 

Inflow Angle 


- 

ALPHA 

Angle-of-Attack 


- 

PHE(I) 

Edgewise Part of i™ Blade Bend- 

- 

- 

PHF(I) 

ing Mode 

Flatwise Part of i^* 1 Blade Bend- 



PHEP(I) 

ing Mode 

Edgewise Slope Part of i# 1 Blade 



PHFP(I) 

Bending Mode 

Flatwise Slope Part of i th Blade 
Bending Mode 

- 

- 


Figure 

Number 


3b 


3c 


3d 


3e 


i 

| 


3f 





The rotor blade radial distributions of edgewise, flatwise and 
] torsional mass moments of inertia, mass, and edgewise and flatwise area 
| moments of inertia are shown in Figure 3g. Care must be exercised in 
the input of these quantities to make sure that the sum of the blade 
| segments equals the blade radius (location 6). 

It either rigid blade flapping or inplane motion is used in the 
rotor analysis, then the rotor blade flapping mass, first and second 
moments of inertia and the blade lag frequency are printed out as can 
be seen in Figure 3h. In this figure are also shown calculations of the 
blade bending mode generalized masses (defined as the blade mass times 
the sum of the squares of the flatwise and edgewise components) and of 
other blade parameters. 

The total number of degrees-of- freedom used in the rotor analysis 
is indicated in Figure 3h. The individual degrees-of- freedoms are 
identified by integers according to the schedule given below. 


Number 

De g ree- o f- F reedom 

Moti on 

1-4 

Blade Bending Modes (up to 4) 

Symmetri c 

5-6 

Blade Torsional Modes (up to 2) 


| 

ji 

7 

Blade Rigid Body Flapping 


| 

8 

Blade Rigid Body Lead- Lag 

\ 

f 

5 

9-16 

Blade Bending Modes 

Cycl i c 

17-20 

Blade Torsional Modes 


:j 

21-22 

Blade Rigid Body Flapping 


j 

23-24 

Blade Rigid Body Lead- Lag 

1 

V ] 

25-29 

Fixed System Modes (Fixed at 5) 


■j 


Thus, for the example in Figure 3h, it is seen that two blade 
bending modes, rigid body flapping and lead- lag and five fixed system 
modes are employed for a total of 17 degrees- of- freedom. 

The maximum number of rotor/fixed system degrees-of- freedom is 29 
(8 blade modes times 3 plus 5 fixed system modes). 

A sample output matrix is presented in Figure 4. For this case, 
location 108 was set to 1. to yield the printout of the compressed 
(17 X 17) rotor/fixed system matrices. The complete 29 X 29 matrices can 
be displayed if location 107 is set to 1. The order of the elements in 
each column follows the schedule shown above. For example, the fourth 
element in the second row corresponds to the lead- lag symmetric stiff- 
ness contribution to the second blade bending mode symmetric equation 
of motion. 

The compressed matrices include rotor aerodynamic contributions. 

They are coupled directly to the bifilar analysis. If no rotor coupling 
is desired, then the compressed matrices are stored for future use. 



4.2 Bi filar Analysis Output 

Typical output formats presenting results from the bifilar analysis 
are shown in Figures 5 through 8. The output parameters from Figures 5a 
j through 5k are common for both linear and non-linear bifilar analyses. 

! The final linear analysis results are presented in Figure 6 while the 
non-linear results can be seen in Figures 7 and 8. 

Figure 5a shows the number of degrees- of- freedom utilized in a 
given computer run, the number of aircraft stations where the response 
is calculated and the printout options requested. For the example 
shown, it is seen that nine fixed system modes, one fixed system absorber 
and one kind each of inplane and vertical bifilars are requested. In 
addition, rotor coupling is employed and the response of four aircraft 
stations is to be analyzed. All printout switches have been activated 
(set at 1.0) to show examples of the output generated. 

The bifilar analysis starts off with the fixed system degrees-of- 
freedom. Then, it expands the fixed system stiffness, damping and mass 
matrices to include in sequence the contributions from the rotor, fixed 
system absorbers, linear inplane bifilars and finally the linear vertical 
bifilars. Then, either a forced response or a time-history solution is 
calculated depending on the input to location 18. The build-up of the 
degrees-of- freedom is shown in Figures 5b through 5k for the stiffness 
matrix only. The damping and mass matrices are handled in an identical 
manner. 

The basic fixed system stiffness matrix is presented in Figure 5b. 

It is seen that for this example the stiffness matrix is a square diagonal 
matrix of order 9. This printout is governed by location 1490. 

The rotor/fixed system stiffness matrix to be coupled with the 
bifilar analysis is shown in Figure 5c and d and is of order 18. This 
matrix is basically the same as that from Figure 4 except that now the 
fixed system degrees-of- freedom appear first and include an extra equation 
corresponding to the yaw degree-of- freedom (which is not present in the 
rotor aeroelastic analysis). Consequently, the total number of degrees- 
of-freedom is 18. This printout option is controlled by location 17. 

Every time a new system component is added, the printout shown in 
Figure 5e is automatically obtained. It shows what the present number 
of degrees-of-freedom (d.o.f.) is (for this example, 9 fixed system 
d.o.f.), the number to be added (12 rotor d.o.f.) and the final system 
d.o.f. (a total of 21 d.o.f.). 

The combined fixed system/rotor stiffness matrix is displayed in 
Figure 5f. Only the first nine equations are shown here for brevity. 

This output is obtained if location 1491 is set to 1.0, 



Next, the fixed system/ rotor coupled system is expanded to include 
j the contribution of the fixed system absorber. In this example, only 
one absorber is used and thus the number of d.o.f. becomes 22, as can be 
1 seen in Figure 5g. Location 1492 controls this printout option. 

If linear inplane bifilar pendulums are used (see location 1), then 
the fixed system/bifilar coupled matrices (mass, damping and stiffness) 
can be obtained for each kind of bifilar (see loc. 1495). The output 

matrices in all cases are square matrices of order 9. The first 6 

d.o.f. correspond to the hub motions (x, y, z, e x , 6y and e z respect- 
ively) and the next three to the inplane bifilar symmetric and two 

cyclic modes. In Figure 5h, the coupled stiffness (QQK) matrix is 

presented. 

The system matrices are now expanded to include the contribution of 
the inplane bifilar pendulums. The stiffness matrix of the fixed system/ 
rotor/fixed system absorber/inplane bifilar coupled system can be seen in 
Figure 5i. Since only one kind of inplane bifilar is used in this example, 
the final number of d.o.f. is increased by 3 for a total of 25. This 
output is generated if location 1493 is set to 1. 

i 

If linear vertical bifilar pendulums are employed in the analysis 
(see loc. 12), the corresponding fixed system/vertical bifilar matrices 
can be displayed (see loc. 1496) as seen in Figure 5j. The format is the 
same as that for the inplane bifilars. 

Now the vertical bifilar d.o.f. are added to the system. The final 
stiffness matrix for this example is shown in Figure 5k. The total d.o.f. 
to be analyzed is 28. This output format is controlled through location 
1494. 1 

The results of the linear bifilar analysis are presented in Figure 
6. The cosine and sine components of the generalized forces appearing 
near the top of Figure 6a are obtained by multiplication of the hub 
forces/moments vectors (see loc 110-119 for the cosine component and loc 
120-129 for the sine component) by the transpose of the rotor hub transfer 
matrix (located in loc. 450-549). The units are lbs and inch-lbs. 

The "GAMMAS" printed out in Figure 6a are the generalized coordinates 
of the rotor/ fixed system absorber/bifilar coupled system and are obtained 
from the forced response solution. First, the cosine component for all 
d.o.f. is printed out and then the sine component. For this example, the 
total number of d.o.f. of the rotor, the fixed system absorber and the 
inplane and vertical bifilars equals 19. Thus, 38 values of "GAMMAS" are 
printed out; the printout switch is location 1497. The units are in inch 
for the fixed system absorber and non-dimensional for the rotor and bifilars.) 




The "GAMMAS" are sorted out according to the system components 
present and printed out accordingly, as shown in Figure 6a and 6b. The 
calculated amplitudes are in inch for the fixed system absorber and in 
degrees for both inplane and vertical bifilar pendulums, all phase 
angles are shown in degrees. The method used to calculate the bifilar 
amplitudes is shown below. 

OUTPUT FORMAT 



Cosine 

Sine 

Amp! i tude 

Phase 

Symmetric Equation 
Cyclic Equation - sine 
Cyclic Equation - cosine 

A °C 

A SC. 

ACC 

A 0S 

A SS 

A CS 

> 

a N-1 

An+1 

phi n 

PHIn-i 

phin+i 


2 2 

A N = A 0C + A os * 57 - 3 °/ NB 


A N-1 = A l 2 + A 2 2 * 57 * 3 °/ NB where A i ~ A cc + A SS & A 2 = A CS " A SC 
A N+1 = + a 4 2 * $7... 30/MB where Ag = A cc - A ss & A 4 - A c$ + A $c 

and NB = number of bifilar pendulums 

The input frequencies are listed out in Figure 6b; the units are 
in Hz. 

The forcing frequency in Hz is obtained by multiplying the forcing 
frequency in cycles/rev (loc 3) by the rotor speed (loc 7) and dividing 
by 60. 

The conversion factor to g's is obtained by dividing the square 
of the forcing frequency in rad/ sec by the factor 386.40. 

The fixed system generalized coordinates are also shown in Figure 
6b. The units are inch. They are utilized to calculate the dynamic 
response of the aircraft stations and of the rotor head which are shown 
in Figures 6b through 6d. As indicated in the printout, the aircraft 
and hub response is in g's. 

Sample results from the time history solution when non-linear 
inplane bifilar absorbers are employed (loc 18 is 1.) are printed in 
Figures 7 and 8. If the total number of degrees-of- freedom exceeds 72, 
then a message will be printed out to that effect and the non-linear 
analysis proceeds to the next case to be analyzed. 





In Figure 7a, the top line lists out the different components 
j degrees-of-freedom requested in a specific computer run. For the 
J example shown, the total number of d.o.f. is 29 which is obtained as 
| f ol 1 ows : 

j System Component 


D.O.F. 


1 . 

Fixed system modes (NF) 

9 ] 

2. 

Rotor d.o.f. (KROTOR) 

12 

3. 

Fixed system absorber (NFABS) 

1 1 

4. 

Kinds of (linear) inplane bifilars (0X3) 

0 1 

5. 

Kinds of vertical bifilar (1X3) 

3 i 

6. 

Number of non-linear inplane bifilars 

4 j 



29 ; 

The 

time history solution proceeds at first to 

collect all the 


acceleration terms on the left hand side of the equations of motion. 

The right hand side contains the stiffness and damping terms and the 
forcing functions. Successive integrations of the accelerations yield 
the velocity and displacement vectors at the next time increment. Then, 
the system accelerations are computed again and the procedure continues 
until a converged time history is obtained or the maximum azimuthal 
position (specified in location 1762) is reached. 

The output formats presented in Figures 7a through 7e are for zero 
azimuth. The initial right hand side (r.h.s.) terms are displayed in 
Figure 7a. They are all zero to start since no initial conditions of 
the state variables displacements and velocities have been located in 
the input locations 1780 through 1939. This printout appears for azimuth 
positions up to 3 degrees. 

j 

The next printout shown in Figure 7a is that of the left-hand-side 
mass matrix whose order is the sum of the six hub d.o.f. and the number 
of non-linear inplane bifilars (loc. 1763) which for this example is 
four. The bifilar force vector is also listed out in this figure. 

Small values appear in this vector due to coupling terms between the 
fixed system and the bifilar pendulums. These outputs are presented for 
azimuth positions up to 5 degrees. 

The rotor head mode shapes appear in Figure 7b. This matrix is an 
image of the input in locations 450 through 549. It is shwon only for 
azimuth positions up to 2 degrees. 

The next printout in Figure 7b is that of the "Expanded Bifilar 
Mass Matrix" whose order is the sum of the number of fixed system modes 
(loc 9) and the non-linear bifilars (loc 1763). Similarly, the "Expanded 
Bifilar Force Vector" is printed out in Figure 7c. Both output formats 
are generated for azimuth ang 1 es up to 5 degrees* . * ***..*. J 




The contributions of the remaining system components are now added 
] to the fixed system/non-linear inplane bifilar system. The resulting 
: mass matrix and force vector can be seen in Figures 7d and 7e respectively. 

The final matrix order is 29, as previously stated in Figure 7a. These 
printouts are obtained for azimuth angles up to 15 degrees. 

The solution vector of the state variables and the non-linear bifilar 
displacements and velocities at zero azimuth is shown in Figure 7e. It 
is printed out for azimuth angles up to 30 degrees. This vector is now 
| used to calculate the right- hand- side terms at the next azimuth position 
(as defined in location 1761) which are displayed at the bottom of 
Figure 7e. 

The time history solution proceeds around the rotor azimuth until 
either the convergence criterion set on the bifilar motions is satisfied 
or the maximum azimuth value (loaded in location 1762) is reached. At 
the end of each rotor revolution, bifilar displacements and hub motions 
are printed out as exhibited in Figure 7f. For the example shown, 16 
rotor revolutions were necessary before the convergence criterion was 
met, i.e. the angular displacements of the first two bifilars (G1 and 
G2) for two successive revolutions must be within .002 radian (or .1146 
degree). 

j} 

: \ 

The rotor/bifilar analysis then proceeds to calculate the harmonic 

response of the non-linear bifilar pendulums displacements, the hub six 
degrees-of-freedom (the order is x, y, z, e x , 6y and e z ) and the aircraft 
station(s) linear motions (x, y and z), as can oe seen in Figures 8a and 
8b. The pendulum output is in degrees while the hub and A/C station(s) 
output is in g's. All phase angles are in degrees. The four rows of 
output describing the hub and A/C response are respectively the cosine 
and sine components, the total amplitude and phase angle. For the 
example shown, the bifilar motions are about 9.8 degrees each; the 
longitudinal (x) hub response is .13864 g's. | 

| 

The final output format of the non-linear analysis results consists 
of the initial values of bifilar, hub and state variables displacements 
and velocities. These results can be loaded into locations 1720 through 
1939 and will be used as initial conditions for the subsequent computer 

run provided location 1498 is set to 1.0. 

. | 

Additional information can be printed out if the control switches 
in locations 10 and 15 are activated. This is only needed if debugging 
of the bifilar analysis calculations is desired. 



SECTION 5 

TEST CASES RESULTS 


! 5.1 Test Cases Description 

The rotor/bifilar coupled analysis has been executed for the follow- 
j ing four test cases: 

Case 1. Includes rotor and uses linear bifilar analysis. 

Case 2. Includes rotor and uses non-linear bifilar analysis. 

Case 3. Doesn't include rotor and uses linear bifilar analysis. 

Case 4. Doesn't include rotor and uses non-linear bifilar analysis. 

These cases test the major program logic paths. 

For each case, the component degrees- of- freedom utilized are listed 
in the chart below. 


COMPONENT D.O.F. 


Test 

Case 

Rotor 

Blade 

Fixed 

System 

Fixed 

Absorber 

Li near 
Inplane 

Bifilar 

Vertical 

Non-Linear 

Bifilars 

Total 

D.O.F 

1 

12 

9 

1 

3 

3 

0 

28 

2 

12 

9 

1 

0 

3 

4 

29 

3 

0 

9 

1 

3 

3 

0 

16 

4 

0 

9 

1 

0 

3 

4 

17 


The rotor degrees-of- freedom include two blade bending modes and 
rigid body flapping and lead-lag motion. Only one kind of inplane and 
vertical linear bifilars is used in the test cases. However, the analysis 
has been checked out for cases utilizing several kinds of bifilar pendulums. 
The inplane bifilars are tuned to 4 per rev. 

5.2 Job Control Language 

The Job Control Language (JCL) needed to execute the coupled program 
on the IBM 370/169 computer system is presented in Appendix A. This 
JCL must be modified for use on the NASA CDC computer system. A brief 
description of the JCL setup is discussed below. 



Description 


Describes job name and characteristics (class, time, etc.)* I 

Executes program module E90BCFIN. 

Locates program module in ET473.SEBBY.L0AD. 

Reads input data from ET473.BIFILAR.DATA (NASARUN) using Unit 5. 
Provides paper output using Unit 6. 

Provides punched cards output using Unit 7. j 

Stores calculated data internally in Unit 8. 

Stores input data for first case in Unit 11 to be used for 
successive cases. f 

Ends JCL setup. j 

S 

Test Cases Input Data j 

The input data needed to run the four test cases is listed in Appendix j 

i 


The input format which must be followed to run multiple cases is 
described in Table 2 below. i 

I 

TABLE 2. Program Multiple Cases Setup j 

Data 
Block 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


The format above is shown for two cases only for brevity. Data blocks 
7 through 10 are repeated for each case. 

If the first case does not use rotor data, then the data blocks 
numbered 1 through 3 above are replaced by a single card as follows: 


Input Data Description Case 

Number 

Rotor blade data (loc. 1-4549) 1 

Last rotor blade data card 
(minus sign in column 1) 

Title for rotor blade data 
Title for bifilar data 
Bifilar data (loc. 1-2199) 

Last bifilar data card 

(loc 2200 = 0. - not last case) ” 

Rotor blade data 2 (last) 

Last rotor blade data card 
(minus sign in column 1) 

Bifilar data 

Last bifilar data card 

(loc 2200 = 1. - last case) V 


Card 

Numb( 

1 

2 

3 

4 

5 

6 

7 

8 


5.3 


-1 110 0 . 



5.4 Test Cases Output 

The results of the rotor/bifilar coupled program for the four test 
cases are presented in Appendix C. Only the important results are shown 
in the Appendix to minimize the size of this report (the actual run 
consisted of 165 pages of output with only the most important printout 
switches being activated). 

Some important results from the bifilar analysis are summarized in 
the Table 3 below. 


TABLE 3. Bifilar Analysis Test Cases Results 


Test 

Inplane Bifilar 

Rotor Hub Response 

Case 

Response 

Longitudinal (X) 

Lateral (Y) | 

Vertical (Z) | 

Number 

Ampi „ 

Phase 

Ampl . 

Phase 

Ampl . 

Phase 

Ampl. 

Phase 

L 


(deg) 

(deg) 

(g's) 

(deg) 

(g's) 

(deg) 

(g's) 

(deg) 

i 

9.29 

-77 

. 134 

113 

.094 

178 

.006 

132 

r 

c 


1. 9.76 

96 

. 139 

-84 

.100 

-20 

.008 

-53 



2. 9.77 

-174 









3. 9.80 

-84 







\ 

! 

4. 9.80 

6 







3 

9.52 

-86 

.184 

113 

.081 

137 

.010 

46 

l\ 


9.57 

88 

.193 

-79 

.086 

-56 

.010 

-134 



9.58 

177 









9.63 

-93 



■ 




N 


9.62 

-2 



.1 





For all test cases, the input force is a 4 per rev fixed system force 
with lateral sine and longitudinal cosine components of 500 pounds. For 
the non-linear bifilar analysis results, the response of each of the four j 
inplane bifilars is shown in the table above (cases 2 and 4). | 

I 

5.5 Test Cases Computer Time f 

The computer total running time for the four test cases was 8 
minutes and 43 seconds. The break-down in computer time per case is 
shown in the following table. 


Test Case 
Number 

Number 
of D.O.F. 

Bifilar 

Analysis 

Computer Time 

Minutes 

Seconds 

1 

28 

Linear 

0 

25 

2 

29 

Non-linear 

7 

07 

3 

16 

Linear 

0 

03 

4 

17 

Non-linear 

1 

08 


i 

— Total ■ :i 

r 

- 43'—== 












From the table above, a comparison between cases 1 and 2 and between 
cases 3 and 4 reveals that the time history analysis requires considerable 
j greater computer time for a complete converged solution than the linear 
| analysis. Also, the computer running time increases tremendously as the 
i number of system degrees-of- freedom is increased when comparing cases 1 
and 3 and 2 and 4. 




SECTION 6 ' ' | 

OVERALL PROGRAM STRUCTURE 

The rotor/bifilar program is basically made up of two parts: the 

rotor aeroelastic analysis and the fixed system/fixed absorber/bifilar 
pendulums analysis. The bifilar anlaysis can be executed with and without 
coupling with the rotor analysis while the opposite is not possible. In 
addition, the bifilar portion of the program can use either a forced 
response solution for linear inplane bifilar pendulums or a time history 
solution for non-linear inplane bifilars. The main purpose of the rotor 
aeroelastic analysis is to provide the rotor blade stiffness, damping 
and mass matrices for coupling with the bifilar analysis. j 

6.1 Segmentation Structure 

Due to the large size of the coupled program, it was necessary to 
implement a segmentation structure to permit operation within a 64K (decimal )( 
for CDC computer use. A basic breakdown of the 10 control segments needed 
is shown in Table 4 below and in the schematic of Figure 9. j 

1 

TABLE 4. Segmentation Structure Description 


Segment 

Leading Fortran 

Number of Segment 

Segment 

Number 

Routine 

Routines 

Description 

1 

SHAKIT 

4 

Controls overall program logic j 
and specifically the rotor 
analysis program flow. 

2 

PRELIM 

34 

Reads rotor input, initializes 
data and performs many basic j 

rotor cal cul ati ons . 

3 

DYNMAT 

12 

Calculates rotor dynamic matrices 

4 

AERMAT 

17 

Calculates rotor aerodynamic j 

matrices. 

5 

EIGER 

2 

Combines dynamic and aero- 
dynamic rotor matrices, com- 
presses and links them to the 
bifilar analysis. 

6 

MAINSV 

3 

Controls bifilar analysis 
program flow. 

7 

SYSCTL 

6 

Calculates contributions from 
fixed system modes, fixed 
absorber, inplane and vertical 
linear bifilar and couples rotor 
matrices. 

8 

HUBIMP 

1 

Computes rotor hub impedance 
and transfer matrices. 

9 

CMPUTE 

4 

Controls generalized forces 
calculations and solves for 
the forced response for the 
. linear bifilar analysis option. : 
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! 10 NLBIF 9 Performs time-history 

solution for the non-linear 

| bifilar analysis option. 


Total = 92 

A total of 92 Fortran routines have been developed for this program: 

69 of them comprise the rotor analysis portion while 23 make up the bifilar 
analysis portion. 

From the schematic presented in Figure 9, it seen that the rotor 
aeroelastic analysis is performed in segments 1 through 5 while the bifilar 
analysis is handled by segments 6 through 10. There is no lateral transfer 
of data between any two segments; data can only be transferred in a vertical j 
sense to segment 1 for the rotor portion and to segments 1 and 6 for the 
bifilar portion. 


Table 5 below lists in alphabetical order the Fortran sub-routines 
and the common blocks needed for each segment. 

TABLE 5. Segmentation Structure Routines and COMMON Blocks 

Segment Number Fortran Subroutine Name COMMON Block Name 

1 SHAKIT 

INTEG 
LOAD IT 
QTFG 
V 


2 PRELIM MATEO 

BLIN4 MATF 

ELI MATR 

ELO MIND 

EXTEND MISC 

E159X MODES 

FILL MSHAPE 

FOLL 0RTH0G 

FREQUN OVUN 

FULL PFMULT 

GMPRDD PICK 

Y MATE I 

3 DYNMAT DMATEX 

BLELPD DMDMAT 

D I SCON DMMMAT 

Y DISINT DMSMAT 

4 AERMAT AERIN5 

AERINT AERIN6 

AERIN1 AERIN7 

AERIN2 AERIN8 

.. AERIN3 AERLST 

v AERIN4 AEROI 



DYNINP 

LAGDAM 


E0F6 

NIMIC 


INDAT 

PRNTSW 


INEIG 

INEIGN 

TMDS 

PINT 

CONT 


POUT 

DAT 


PRODM 

EMATI 


PROUT 

EMATO 


REMOVE 

FMAT 


ROOTX 

FREQ 


SECAER 

PHTNO 


SIMLIN 

PMAT 


SKIPLN 

PRAM1 


SORTAB 

TORFIN 


STDEFL 

WORKA 


DYNINT 

DYNOUT 


DYNIN1 

DYNIN2 

DYNLST 

NAMIC 


AEROII 

AEROI 


AMATEX 

AER02 


AMDMAT 

AER03 


AMSMAT 

AER04 


BLELPA 

AER05 

CDCAER 




Segment Number 

Fortran Subroutine 

Name COMMON Block 

Name 

5 

EIGER 




COMPRSS 



6 

MAINSV 

NDOF 

TOTMAT 



INCOND 

NLDAT1 

XFRDAT 

\ 

! 

INPUTV 

PRSWTH 


7 

SYSCTL 





ADDOFR 





FIXABS 





FIXSYS 





LINBIF 



\ 

f 

LVBIF 



8 

HUBIMP 



9 

CMPUTE 





FORCER 





GENFOR 



\ 

/ 

OUTPUT 



1( 

) 

NLBIF HARMON 

HARM 




BIFEXP INTEQ 

NLDAT2 




BIFILR OUT 





COMB IN RHS 



\ 

! 

CONVER 




All the Fortran sub-routines listed above except four are computer 
independent. The computer dependent routines are: SHAKIT (segment 1), 

LOADIT (segment 1), PRELIM (segment 2) and MAINSV (segment 6). This is 
due to CDC computer requirements for identification of the main routine, 
file read error and end of file transfers and word size definition for 
alpha-numeric read statements. The program coding allows the programmer 
to convert easily to the IBM 370/168 computer system by commenting out the 
appropriate block(s) of statements. 

6.2 Flow Diagrams 

Computer logic flow diagrams are presented in Figures 10 through 16 
for the 10 segments. A summary of the flow charts is provided in the 
table below. 

Figure Number Flow Chart Description Segment Number(s) 


10 

Main Program Flow Chart 

1 & 5 

11 

"PRELIM" Flow Chart 

2 

12' 

"MODES" Flow Chart 

2 

13 

"DYNMAT" Flow chart 

3 

14 

"AERMAT" Flow Chart 

4 

15 

Bifilar Analysis Flow Chart 

6 & 8 

16a 

"SYSCTL" Flow Chart 

7 

16b 

"CMPUTE" Flow Chart 

9 

16c 

"NLBIF" Flow Chart 

10 




It should be noted that "IMSL" routines are needed for the calculations 
performed in "HUBIMP" and "NLBIF" . Table 6 below lists the 10 "IMSL" 

| routines used in the bifilar analysis portion of the coupled program: 

| TABLE 6. IMSL Routines 

j 

No. IMSL Routine Name 


1 LEQT2F 

2 LINV2F 

3 LUDATF 

4 LUELMF 

5 LUREFF 

6 UERTST 

7 UGETIO 

8 VXADD 

9 VXMUL 

10 VXSTO 

These routines must be supplied by the government for successful! 
operation of the rotor/bifilar coupled program. 

6.3 "COMMON" Blocks 

The "COMMON" blocks used in the rotor/bifilar analysis are presented 
in Figure 17 as they appear in each Fortran sub-routine. Both routines 
and "COMMON" blocks are listed alphabetically for easy reference. 


I 


l 

| 
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SECTION 7 

I 

SUBROUTINE DESCRIPTIONS 


The Fortran subroutines needed to execute the rotor/bifilar coupled 
program are described in detail in this section. For each routine, the 
following information is provided, where applicable: 

1 . Name 

2. Purpose 

3. Method 

4. Usage 

5. Subroutines Called 

6. Error Returns 

7. Restrictions 


In addition, an alphabetical listing of the routines and the correspond- 
ing page numbers are shown in Table 7 below for easy reference. 




TABLE 

7. Program Subroutines Listing 




No. 

Name Page 

No. 

Name 

Page No. 

Name Page 

No. 

Name 

Page 

1 

ADDOFR 

24 

CMPUTE 

47 

FORCER 

70 

NLBIF 

131 

2 

AERINT 

25 

COMBIN 

48 

FREQUN 

71 

ORTHOG 

132 

3 

AERIN1 

26 

CONVER 

49 

FULL 

72 

OUT 

133 1 

4 

AERIN2 

27 

DISCON 

50 

GENFOR 

73 

OUTPUT 

134 i 

5 

AERIN3 

28 

DISINT 

51 

GMPRDD 

74 

OVUN 

135 | 

6 

AERIN4 

29 

DMATEX 

52 

HARMON 

75 

PFMULT 

136 

7 

AERIN5 

30 

DMDMAT 

53 

HUB IMP 

76 

PICK 

138 

8 

AERIN6 

31 

DMMMAT 

54 

INCOND 

77 

PINT 

139 

9 

AERIN7 

32 

DMSMAT 

55 

INPUTV 

78 

POUT 

140 

10 

AERIN8 

33 

DYNINT 

56 

INTEG 

79 

PRELIM 

141 | 

11 

AERLST 

34 

DYNIN1 

57 

INTEQ 

80 

PRO DM 

144 

12 

AERMAT 

35 

DYNIN2 

58 

LINBIF 

81 

PROUT 

145 1 

13 

AERO I 

36 

DYNLST 

59 

LOADIT 

82 

QTFG 

146 

14 

AEROII 

37 

DYNMAT 

60 

LVBIF 

83 

REMOVE 

147 

15 

AMATEX 

38 

EIGER 

61 

MAINSV 

84 

RHS 

148 • 

16 

AMDMAT 

39 

ELI 

62 

MATE I 

85 

ROOTX 

149 

17 

AMSMAT 

40 

ELO 

63 

MATEO 

86 

SECAER 

150 

18 

BIFEXP 

41 

EXTEND 

64 

MATE 

87 

SHAKIT 

151 

19 

BIFILR 

42 

E159X 

65 

MATR 

88 

SIMLIN 

152 ! 

20 

BLELPA 

43 

FILL 

66 

MIND 

89 

SKIPLIN 

153 

21 

BLELPD 

44 

FIXABS 

67 

MI SC 

90 

SORTAB 

154 

22 

BLIN4 

45 

FIXSYS 

68 

MODES 

91 

STDEFL 

155 

23 

CMPRSS 

46 

FOLL 

69 

MSHAPE 

92 

SYSCTL 

157 


U3" ] 



NAME: 


ADDOFR 


I 

j PURPOSE: 
METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 

54 j 


To add degrees- of- freedom to the mass, damping and f 

stiffness system matrices in the bifilar analysis. 

| 

| 

The new matrices to be added of order NL X NL are 
split up in four parts: 

1. An upper diagonal matrix of order (NL-NA)X(NL-NA) 

2. A lower diagonal matrix of order NA X NA 

3. An upper off diagonal matrix of order (NL-NA) X NA 

4. A lower off diagonal matrix of order NA X (NL-NA) 

The lower diagonal matrices are added to the original 
matrices of order NP X NP, which now become of order 
(NP + NA) X (NP + NA"), The upper off diagonal matrix 
is pre-multi plied by the fixed system transfer matrix 
XPH (locations 450-549) while the lower off diagonal 
matrix is post-multiplied by XPH. The upper diagonal 
matrix is pre-multi plied and post-multiplied by XPH. j 


CALL ADDOFR (NL, NA, NP) 

NL - Order of matrices to be added 

NA = Number of degrees-of- freedom to be added 

NP = Present order of matrices (after ADDOFR, 
order of matrices is NP + NA) 

None 

None 

None 




NAME: 

j PURPOSE: 
METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


AERINT 


To set initial values of all blade aerodynamic integrals j 
to zero and to control the calculation of the aerodynamic! 
integrals needed to construct the aerodynamic damping 
and stiffness matrices. 


Aerodynamic integrals independent of the number of 
bending modes are calculated by calls of AERIN1 and 
AERIN7 for lift, AERIN3 and AERIN8 for drag, and 
AERIN4 for pitching moment. Integrals involving bend- 
ing mode dependent functions are calculated in AERIN2, 
AERIN5 and AERIN6 for lift, drag and pitching moment 
respectively. 


CALL AERINT 

I 

AERIN1, AERIN2, AERIN3, AERIN4, AERIN5, AERIN6, AERIN7, I 
AERIN8 1 


None 


RESTRICTIONS: 


Due to computer storage restrictions, the aerodynamic 
integrals had to be calculated in 8 separate routines, 
i.e. AERIN1 through AERIN8. 



AERIN1 


NAME: 


PURPOSE: 


METHOD: 


'] USAGE: 


To calculate the aerodynamic integrals which contain 
thrust derivatives and are independent of the number 
of blade bending modes. 


These integrals are formed in AEROI from the product 
of 5 radial functions, and are referred to as un- 
subscripted ATK and AT1K integrals. 

Since all the integrals involve thrust derivatives, 
the upper limit of integration is the blade radius 
multiplied by the tip loss factor. 


CALL AERIN1 


SUBROUTINES AEROI I, AEROI 

CALLED: 


ERROR RETURNS: None 


j RESTRICTIONS: None 


: 
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NAME: 

1 PURPOSE: 

1 

I METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


AERIN2 


To calculate the aerodynamic integrals which contain 
bending mode dependent functions and thrust derivatives. 


I 


| 


These integrals are formed in AEROI from the product 
of 5 radial functions. They are referred to as ATI I 
and AT2I for doubly subscripted integrals and ATJ, 

AT1J, AT2J, AT3J for singly subscripted integrals. 

Since all the integrals involve thrust derivatives, the j 
upper limit of integration is the blade radius multiplied] 
by the tip loss factor. 


CALL AERIN2 
AEROI I, AEROI 


None 


None 


\ 

! 

I 

i 

! 

| 

I 

j 

j 


| 


I 


l 

j 

l 


1 

i 


1 



NAME: 


PURPOSE: 


| METHOD: 


j USAGE: 

.1 

| 

f 

SUBROUTINES 
I CALLED: 


ERROR RETURNS: 


RESTRICTIONS: 


AERIN3 


To calculate the aerodynamic integrals which contain 
drag derivatives and are independent of the number of 
blade bending modes. 


These integrals are formed in AEROI from the product 
of 5 radial functions, and are referred to as un- 
subscripted ADK and AD1K integrals. 

Since all the integrals involve drag derivatives, 
they are computed over the whole blade. 


CALL AERIN3 
AEROI I, AEROI 


None 


None 




j NAME: 

| 

j PURPOSE: 

j 

j 

i 

j METHOD: 

\ 

j 

j USAGE: 

SUBROUTINES 

CALLED: 

■I 

ERROR RETURNS: 
RESTRICTIONS: 




AERIN4 


To calculate the aerodynamic integrals which contain 
pitching moment derivatives and are independent of the 
number of blade bending modes. 


These integrals are formed in AEROI from the product 
of 5 radial functions, and are referred to as un- 
subscripted AM1K, AM2K, AM3K integrals. 

Since all the integrals involve pitching moment deriv- 
atives, the upper limit of integration is the blade 
radius multiplied by the tip loss factor. 


CALL AERIN4 


AEROI I, AEROI 


None 


None 






AERIN5 


NAME: 


i PURPOSE: To calculate the aerodynamic integrals which contain 

bending mode dependent functions and drag derivatives. 

I 

{ 

| METHOD: These integrals are formed in AEROI from the product 

of 5 radial functions. They are referred to as ADI I 
and AD2I for doubly subscripted integrals and ADJ, AD1J 
AD2J and AD3J for singly subscripted integrals. 

Since all the integrals involve drag derivatives, they 
are computed over the whole blade. 


I USAGE: 


CALL AERIN5 


i SUBROUTINES AEROI I, AEROI 

CALLED: 


! ERROR RETURNS: 

None 

j RESTRICTIONS: 

None 


\ 


f. 


i 




60 






1 NAME: AERIN6 


PURPOSE: 


METHOD: 


| 

| 

\ 


! 

USAGE: 


I SUBROUTINES 
CALLED: 


To calculate the aerodynamic integrals which contain 
one bending mode dependent functions and pitching 
moment derivatives. 


These integrals are formed in AERO I from the product 
of 5 radial functions, and are referred to as singly ! 
subscripted AM1J and AM2J integrals. Since all the I 

integrals involve pitching moment derivatives, the upper | 
limit of integration is the blade radius multiplied by 
the tip loss factor. I 


CALL AERIN6 
AEROII , AERO I 


! 


ERROR RETURNS: None 


I RESTRICTIONS: None 


| 

! 

j 


i 

j 
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NAME: 

PURPOSE: 


AERIN7 


To calculate the aerodynamic integrals which contain 
thrust derivatives and are independent of the number 
of blade bending modes. 


METHOD: These integrals are formed in AEROI from the product 

of 5 radial functions, and are referred to as unsubscripted 
AT2K and AT3K integrals. 

Since all the integrals involve thrust derivatives, the 
upper limit of integration is the blade radius multiplied ! 
by the tip loss factor. j 


USAGE: CALL AERIN7 


! SUBROUTINES 
CALLED: 


ERROR RETURNS: 


RESTRICTIONS: 


AEROI I, AEROI 

None 

None 


\ 


] 


\ 

i 


I 

■I 

i 

j 

i 

s 

1 

j 

! 

| 


i 

| 





NAME: 

PURPOSE:' 

| METHOD: 

I 

i 

| 

USAGE: 

I SUBROUTINES 
CALLED: 

ERROR RETURNS: 

■{ 

\ 

1 

: RESTRICTIONS: 

'j 

■I 

i 

i 

| 


L 




AERIN8 


To calculate the aerodynamic integrals which contain 
drag derivatives and are independent of the number of 
bending modes. 


These integrals are formed in AEROI from the product of 
5 radial functions, and are referred to as unsubscripted 
AD2K and AD3K integrals. Since all the integrals involve 
drag derivatives, they are calculated over the whole 
blade. I 


CALL AERIN8 
AEROI I, AEROI 


1 

: 


None 


None 


:< 




j 

j 


j 

j 
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I NAME: AERLST 

\ 


PURPOSE: 

To print the aerodynamic integrals. 


METHOD: 

If the print option is set at 3 or 4, then the 
are printed. Otherwise control is returned to 
with no integrals printed. 

i ntegral 
AERMAT 


Only those integrals calculated are printed. 


USAGE: 

CALL AERLST 


SUBROUTINES 

CALLED: 

None 


ERROR RETURNS: 

None 


RESTRICTIONS: 

None 



'] 

1 


i 


: 




j NAME: 
j PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 



AERMAT 


To control the calculation of the aerodynamic damping 
and stiffness matrices. 


"5 


j 


5 


The aerodynamic integrals needed to calculate the 
matrices are evaluated in AERINT, stored in COMMON 
blocks AER01 through AER05, and printed out by AERLST. 
AERMAT then makes calls to AMDMAT to calculate the 
damping matrix, and to AMSMAT to calculate the stiff- 
ness matrix. Then, it extends the matrices to include 
blade torsional bending mode terms in AMATEX and calls 
BLELPA to calculate elastic torsional contributions. 


CALL AERMAT 

AERINT, AERLST, AMDMAT, AMSMAT, AMATEX, BLELPA 
None 


None | 

I 






AEROI 





NAME: 


PURPOSE: 

To calculate the integral of the function y = F1(R)* 
F2(R)*F3(R)*F4(R)*F5(R) in the interval (rl, r2). 

METHOD: 

The 5 radial functions are multiplied together at each 
radial station specified. It should be noted that the 
second argument in the AEROI calling sequence is 
designated complex and corresponds to an aerodynamic 
derivative function. The integral is evaluated in two 
parts, real and imaginary, by calls to INTEG. If there 
are no unsteady aerodynamic corrections, then the 
imaginary part of the integral is set to zero. 

USAGE: 

CALL AEROI (FI, F2, F3, F4, F5, ANS) 


FI, F3, F4, F5 = Real functions defined at the radii 

in R from AEROI I. One or more may be 
equal to 1.0 at every R. 

.■ 

i 

1 

i 

F2 = Complex function defined at the radii 

in R. If no unsteady aerodynamic 
correction has been applied, the 
imaginary part will be zero. 

i 

I 

i 

i 

; 

ANS = The complex integral of the product 

of FI through F5. If no unsteady 
aerodynamic corrections have been 
applied, then the imaginary part will 
be zero. 

i SUBROUTINES 
! CALLED: 

i 

INTEG 

ERROR RETURNS: 

1 

i 

None 

j 

S RESTRICTIONS: 

None 







NAME: 


AEROII 


PURPOSE: 


j 

: METHOD: 

! 

i 


i 

i USAGE: 


; 

j 

■s 

j 

j 

s 

| 

! SUBROUTINES 
! CALLED: 

j 

! ERROR RETURNS: 

I 

! RESTRICTIONS: 


To set up calculations of the integral of the function 
y = a(R)*b(R)*c(R)*d(R)*e(R) over the interval (rl,r2). 


I 

This routine sets the radial values where the functions 
are defined, the limits of integration, the number of 
subdivisions to be used in the trapezoidal rule inte- 
gration method and the switches for unsteady aerodynamics J 

j 

Subsequent calls to AEROI supply the 5 functions whose j 
product is to be integrated. 


CALL AEROII (R, N, RL, RU, K, CIR, CIRN) 

R = Array of radii at which the functions are defined. 

N = The number of points in R and in each function 

array. 

RL = The lower limit of integration. 

RU = The upper limit of integration. 

K = The number of subdivisions to be used in perform- 

ing the integration by the trapezoidal rule. 

CIR = 1 No circulatory unsteady aerodynamics used. 

CIRN = 1 No nonci rculato ry unsteady aerodynamics used. 

None 

None 

1. N must be equal or less than 25. 

2. CIR and CIRN have been set to 1 for no unsteady 
aerodynamics. 






| NAME: 

| 

! PURPOSE: 

j 

| METHOD: 

] 

1 

.j 

! 

| USAGE: 

j 

! SUBROUTINES 
i CALLED: 

3 ERROR RETURNS: 

I RESTRICTIONS: 


AMATEX 


To expand the aerodynamic matrices to include the blade 
elastic torsional mode. 


The original aerodynamic damping (AMD) and stiffness 
(AMS) matrices are increased by 3 rows and columns 
and redefined as AMDN and AMSN to accomodate one 
collective and two cyclic modes associated with the 
blade elastic torsional mode. 


CALL AMATEX 


None 


None 


The matrices are limited to 30 x 30. 





NAME: 


AMDMAT 


j PURPOSE: 

| METHOD: 

I 

i 

i 

i 

| USAGE: 

j 

1 

! SUBROUTINES 
| CALLED: 

j ERROR RETURNS: 

RESTRICTIONS: 


To calculate the aerodynamic damping matrix. 


The matrix is calculated from the expressions given 
in Reference 1 using the aerodynamic integrals cal- 
culated in AERINT. 


CALL AMDMAT 
None 


None 


None 


I 

X 






NAME: 


AMSMAT 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To calculate the aerodynamic stiffness matrix. 


The elements are calculated from the expressions given 
in Reference 1 using the aerodynamic integrals cal- 
culated in AERINT. 


CALL AMSMT 
None 


None 


RESTRICTIONS: 


None 



j NAME: 

j 

PURPOSE: 

j 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


BIFEXP 


To transfer the non-linear bifilar mass matrix and the 
bifilar and rotor hub force vector to fixed system 
coordinates. 


The bifilar mass (square matrix of order 6+NBIF) and 
the force vector (of order 6+NBIF) are transferred to 
the fixed system coordinates by proper multiplications 
with the fixed system mode shapes vector (loc. 450-549). 
The final expanded mass matrix and force vector have 
the dimensions of NF + NBIF. The method used is 
similar to that discussed in routine ADDOFR. The 
expanded mass and force vector are passed through 
labelled C0MM0N/NLDAT2. 

They are printed out for azimuth positions up to 5 
degrees. The fixed system mode shapes vector is 
printed out for azimuths up to 2 degrees. 


CALL BIFEXP (NRHS, NF) 

NRHS = Total number of non-linear bifilars plus 6 
(maximum is 18). 

NF = Number of fixed system modes (maximum is 16). 
None 

None 


Expanded matrix and force vector cannot have dimensions 
greater than 28. 



NAME: 


j PURPOSE: 


1 METHOD: 


USAGE: 


SUBROUTINES 

CALLED: 


ERROR RETURNS: 



BIFILR 


To calculate the non-linear bifilar mass matrix and 
the bifilar and rotor hub force vector. 


For the time history solution, the bifilar acceleration 
terms are left on the left- hand- side of the equations 
of motion while the r.h.s. contains the bifilar 
damping and stiffness terms and the input hub forces. 
The bifilar mass matrix, S (dimensioned 6+NBIF), is 
calculated and stored in labelled COMMON/INEIGN. 

Then, the bifilar damping and stiffness terms are 
used to develop the r.h.s. vector, T (dimensioned 
6+NBIF). Next, the rotor hub forces, loaded in 
locations 1600 through 1719, are evaluated and added 
to the bifilar contributions in vector T. 


The hub forces are added to the system gradually 
according to the ramp factor (location 1760) and the 
azimuth position. The vector T is passed through 
labelled C0MM0N/NLDAT2. 


The order of the degrees- of- freedom in this routine 
is: 6 fixed system d.o.f. (x, y, z, e x , 0y, 0 Z ) 

and non-linear bifilars d.o.f. (NBIF-location 1763). 
Thus, the maximum number of d.o.f. is 18. 

For azimuth positions up to 5 degrees, the mass matrix 
and the force vector are printed out. 

CALL BIFILR (NREV) 

NREV = Revolution number- location 1762 divided by 
360. 


None 


•| 


None 


| 


RESTRICTIONS: 


Number of input harmonics of rotor hub forces is 
limited to 10. 



BLELPA 


j NAME: 

j 

j 

1 PURPOSE: 


METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To calculate the blade elastic pitch aerodynamic 
contributions. 


The aerodynamic damping and stiffness matrix elements 
are calculated using the appropriate expressions for 
blade pitch in Reference 1. All the terms are 
multiplied by the torsional mode shape, PH(r), 
which is a function of blade radius, except the 
blade pitch terms which are multiplied by the square 
of the mode shape. 


CALL BLELPA 


None 


None 


RESTRICTIONS: 


None 




NAME: 

BLELPD 

PURPOSE: 

To calculate the blade elastic pitch dynamic 
contributions. 

METHOD: 

j 

The dynamic mass, damping and stiffness matrix 
elements are calculated using the equation for the 
blade pitch in Reference 1. All the terms are 
multiplied by the torsional mode shape, PH(r), which 
is a function of blade radius, except the blade pitch 
terms which are multiplied by the square of the 
mode shape. 

USAGE: 

CALL BLELPD 

SUBROUTINES 

CALLED: 

j 

None J 

j 

ERROR RETURNS: 

i 

j 

None j 

RESTRICTIONS: 

| 

None j 



■i 







BLIN4 


NAME: 

| 

i PURPOSE: 
I METHOD: 


To provide a bi variant table lookup with linear inter- 
polation for the airfoil data. 


The Mach number entries are searched to find K such 
that 

m k-i < M « m k 

The angles of attack in the Mj<_i and Mk tables are 
searched to find I and J such that 

A K-1,I-1 < A < A K-1,I 
and 

A K,J-1 < A < A K,J 

Then the coefficient at M^i and A is given by 

C K-1 = ^ A " A K-1,I-1^*^ C K-1,I ” C C-1,I-1^ 

^ A K-1,I ' A K-1,I-1^ + C K- 1,1-1 
and the coefficient at M K and A is given by 

C K = (A - A Kj j. 1 )*(C KiJ - C KjJ . 1 )/(A KjJ - A KjJ _ 1 ) + 
C M-1 

Finally, the coefficient at M and A is obtained from 

C = (M - M k _ 1 )*(C k - C k _ 1 )/(M k - M^) + C |< _ 1 

DC = Ck - C K _i 
DM M|< - M k _i 


: 

i 

i 

; 

| 

j 

i 

■ 

■i 


j 

] 

j 


:i 

.! 

i 




) 



USAGE: 


CALL BLIN4 (T, M, K, X, Y, Z, D, L) 


j 

1 

;j 


T 


M,K 


X 


A two-dimensional array containing the coef- 
ficient, angle of attack pairs. The first point 
in each column represents the number of pairs, 
the second point in each column represents the 
Mach number, and the remaining points in the 
column are the pairs for that Mach number. 

I 

The dimensions of T 

M is the maximum number of Mach numbers in the 
table. j 

K is the maximum column length. j 

Angle of attack. 


Y = Mach number. 


Z = Returned coefficient. 

D = Derivative with respect to M. 

L = Error switch, 

1 No error 

2 Mach number not spanned 

3 Angle of attack not spanned 


i SUBROUTINES None 

; CALLED: 

; ERROR RETURNS: See above 


i RESTRICTIONS: 1. Number of Mach numbers specified must be equal or 
| greater than 2 and cannot be greater than 12. 

2. Number of coefficient/angle of attack pairs must 
i be equal or greater than 5 and cannot be greater 
i than 35. 



NAME: 


CMPRSS 


PURPOSE: 

METHOD: 

\ 

\ USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


To compress a square matrix. 


i 

I 

•I 


The returned matrix of dimension L x L is produced by 
taking the first L rows and columns of the input matrix 
of dimension N x N. 


CALL CMPRSS (A, N, B, L) 
A = Input matrix. 


N = Dimensions of A. 


B = Output matrix, which may have the same location 
as A. 

L = Dimensions of B. 


None 


None 


1 


If L is equal or greater than N, no matrix compression 
is done. 


i 


j 

i 









CMPUTE 


NAME: 

| 

i PURPOSE: 

I 

| 

] METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To control calculations of the generalized forces 
and of the forced response for the linear bi filar 
analysis. 


First, the generalized forces are calculated in GENFOR. | 
Then, the forced response of the linear system of 
equations is obtained in FORCER. The results of the 
analysis are printed out in a call to OUTPUT. 


CALL CMPUTE 


GENFOR, FORCER, OUTPUT 


None 


RESTRICTIONS: 


None 



COMBIN 


NAME: 


To combine the non-linear bifilar mass matrix and 
the force vector developed for the non-linear bifilar 
analysis with the corresponding matrix and vector from 
the linear bifilar analysis. 


The mass matrix and force vector calculated in BIFEXP 
are added to the results of the linear analysis from 
SYSCTL. The fixed system d.o.f. are affected. The 
combined mass matrix and force vector include all 
the degrees-of- freedom (up to the maximum of 72). 

If forces are input to one or two additional aircraft 
stations (as specified in locations 1767 and 1940 
through 2179), their corresponding contributions 
to the final force vector are included by pre-multi ply- 
ing the input forces by the appropriate aircraft 
station mode shapes (see locations 550-749). 

The final combined mass matrix and force vector are 
stored in labelled C0MM0N/NLDAT2. 


USAGE: 

CALL COMBIN (NBML, NRHS, NF, NREV) 


NBML 

= Fixed system modes plus non-linear bifilars j 

(maximum is 28). j 


NRHS 

= Total number of d.o.f. not including non- 
linear bifilars (maximum is 60). | 


NF 

= Number of fixed system modes (maximum is 16). 


NREV 

= Revolution number - loc 1762 divided by 360. 

SUBROUTINES 

CALLED: 

None 

| 

ERROR RETURNS: 

None 

I 

'i 

RESTRICTIONS: 

1 . 

; 

Total number of degrees- of- freedom is 72 (60 
from linear analysis and 12 for the maximum 
number of non-linear bifilars). 


2. 

| 

Harmonic force inputs are possible for only 2 
aircraft stations. 


79 


PURPOSE: 

METHOD: 



I NAME: 

I 

i 

] PURPOSE: 

| 

I METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


CONVER 


To test on the convergence of the time history 
solution for two successive revolutions. 


This routine checks the difference in the displace- 
ments of the first two non-linear bifilar pendulums 
obtained for two successive revolutions. If both 
differences are within .002 radian (corresponding 
to .1146 degree), then the convergence criterion 
is satisfied; IER (see below) is set to 1 and the 
program returns to NLBIF to calculate one more 
revolution after which the harmonic analysis is 
performed and printed out. If one of the differences 
is greater than .002 radian, the time history analysis 
proceeds to the next revolution. 

The revolution number, the two bifilar displacement 
differences and the rotor hub x, y, e z , X, y motions 
are all listed out for each revolution. 


CALL CONVER (I, IER) 

I = Rotor revolution number - location 1762 divided 
by 360. 

IER = Convergence criterion - it is met if equal to 1. 


None 


None 


RESTRICTIONS: 


None 



NAME: 

D I SCON 

! 

j 

PURPOSE: 

1 

! 

To adjust a blade stepwise function so that it is 
defined only over the blade. j 

METHOD: 

| 

The input segments are extended or truncated so that 
the first segment starts at the offset and the final 
segment finishes at the blade radius. Function values 
are unaltered but may be discarded if the segment 
they refer to is completely outside the blade. 

USAGE: 

CALL D I SCON (KSTAR, DELTAR, NSTAR , E, R, RBAR, KBAR, 
N1BAR) 


KSTAR = 

: 

Array of input stepwise function values. 


DELTAR = 

Array of segment lengths over which KSTAR 
is defined. j 


NSTAR = 

Number of entries in KSTAR and DELTAR. 


E 

! 

Blade offset. j 


R 

Blade radius. j 


RBAR 

Array of adjusted segment lengths. 


KBAR 

Array of stepwise function values defined 
over RBAR. J 


N1BAR = 

Number of entries in RBAR and KBAR. 

SUBROUTINES 

CALLED: 

None 

| 

$ 

ERROR RETURNS: 

None 


j RESTRICTIONS: 

None 

i 

] 

.! 

| 




! 
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NAME: 


DISINT 


PURPOSE: 

METHOD: 


USAGE: 


To perform the integration of the product of a step 
function and a radial function over the blade length. 


The stepwise function is first adjusted by DISCON so 
that it is defined over a segment distribution which 
starts at the offset and finishes at the blade radius. 
The integration of the radial function is then performed 
separately over each of these segments and multiplied 
by the value of the stepwise function for that segment. 
The final answer is the sum of these results. 


i.e. / N f(r)S(r)dr = S. / 1 f(r)dr+S 9 r 2 f(r)dr+... 

E 1 E L R x 

R R 

S k / f(r)dr+ S N / N f(r)dr 

R K-1 R N-1 

where 

J.L 

S K = the value of the stepwise function over the K tn 
segment. 

f(r) = the radial function 


I 

\ 


l 

I 

i 


i 

! 


CALL DISINT (COMP, R, N, FPPP, RPPP, NM2, E, CR, SUM) 

COMP = Array of radial function values. 

R = Array of radii at which COMP is defined 

N = Number of points in COMP and R. 

FPPP = Array of stepwise function values. 

RPPP = Array of segment lengths over which FPPP 
is defined. 



USAGE: 

| 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 



NM2 = Number of entries in FPPP and RPPP. 

| 

E = Blade offset. 

CR - Blade radius. 

SUM = Integral of the product of the two functions. 

D I SCON, INTEG I 

J 


None 

None 


i 

j 


! 
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NAME: 

PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


DM AT EX 


To expand the dynamic matrices to include the blade 
elastic torsional mode. 


The original dynamic mass (DMM), damping (DMD) and 
stiffness (DMS) matrices are increased by 3 rows and 
columns and redefined as DMMN, DMDN and DMSN respectively 
to accomodate one collective and two cyclic modes 
associated with the blade elastic torsional mode. 


CALL DMATEX 
None 


None 


RESTRICTIONS: 


The matrices are limited to 30 x 30. 



NAME: 

PURPOSE: 

METHOD: 

| USAGE: 

! SUBROUTINES 
j CALLED: 

j 

I ERROR RETURNS: 

j 

| RESTRICTIONS: 

'I 

\ 


DMDMAT 

To calculate the dynamic damping matrix. 

The matrix is calculated from the expressions given 
in Reference 1, using the dynamic integrals calculated 
in DYNINT. 

CALL DMDMAT 

None 

None 

None 






NAME: 


! PURPOSE: 

j 

| METHOD: 


USAGE: 


| SUBROUTINES 
j CALLED: 


I ERROR RETURNS: 



DMMMAT 

To calculate the dynamic mass matrix. 

The matrix is calculated from the expressions given 
in Reference 1, using the dynamic integrals calculated 
in DYNINT. 

CALL DMMMAT 

None 

None 


RESTRICTIONS: 


None 







NAME: 


DMSMAT 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


To calculate the dynamic stiffness matrix. 

The matrix is calculated from the expressions given 
in Reference 1, using the dynamic integrals calculated 
in DYNINT. 

CALL DMSMAT 

None 

None 

None 





NAME: 

PURPOSE: 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


DYNINT 


To set the initial values of the dynamic integrals 
to zero and to control the calculation of the integrals 
needed to construct the dynamic mass, damping, and 
stiffness matrices. 


The dynamic integrals are obtained from the integration 
of the product of a radial function and a stepwise 
function over the blade length. All integrals have 
the blade mass or the edgewise, flatwise, and torsional 
mass moments of inertia as their stepwise function. 

If the radial function is independent of the blade 
bending modes, then the integrals are calculated in 
DYNIN1. Integrals with singly and doubly subscripted 
radial functions are handled in DYNIN2. 


CALL DYNINT 
DYNIN1, DYNIN2 


None 


RESTRICTIONS: 


None 




I NAME: DYNIN1 


PURPOSE: To calculate the dynamic integrals whose setpwise 

function is either the blade mass or the edgewise, 
flatwise, or torsional mass moment of inertia and 
whose radial functions contain no blade bending mode 
dependent quantities. 


METHOD: These integrals are formed in DISINT from the product 

of a radial function and a stepwise function in- 
tegrated over the blade length. The radial function 
may itself be a product of radial functions. 

The integrals which are a function of blade mass are 
referred to as unsubscripted BMK. The integrals 
dependent on the blade mass moments of inertia are 
referred to as unsubscripted BIK. 


USAGE: CALL DYNIN1 


; SUBROUTINES DISINT 

CALLED: 

) 

I ERROR RETURNS: None 


RESTRICTIONS: None 


I 

i 




NAME: 


DYNIN2 


PURPOSE: 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To calculate the dynamic integrals whose stepwise 
function is either the blade mass or the edgewise, 
flatwise, or torsional mass moment of inertia and 
whose radial function contains one or two blade bending 
mode dependent quantities. 

These integrals are formed in DISINT from the product 
of a radial function and a stepwise function. The 
integrals which are a function of blade mass are re- 
ferred to as BMJ and BMI for singly and doubly sub- 
scripted values respectively. The integrals dependent 
on the blade mass moments of inertia are referred to as j 
BIJ and BII for singly and doubly subscripted values 
respectively. 


CALL DYNIN2 
DISINT 


None 


RESTRICTIONS: 


None 




NAME: 


DYNLST 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To print out the dynamic integrals. 

If the print option is set to 3 or 4, the integrals 
are printed. Otherwise, control is returned to DYNMAT. 

Only those integrals calculated are printed. 

CALL DYNLST 
None 

None 


RESTRICTIONS: 


None 



j NAME: 

| PURPOSE: 

I METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 



DYNMAT 



I 


To control the calculation of the dynamic mass, damping 
and stiffness matrices. 


The dynamic integrals needed to calculate the elements 
of the matrices are evaluated in DYNINT, stored in 
labelled COMMON/ DYNOUT, and printed out by DYNLST. 

DYNMAT then makes calls to DMMMAT, DMDMAT and DMSMAT 
to calculate the mass, damping and stiffness matrices 
respectively. Then, it expands the matrices in DMATEX 
to include the blade torsional elastic mode terms, which 
are calculated in BLELPD. 

In addition, several blade parameters are calculated 
and printed out according to the degrees-of- freedom 
being used. 

CALL DYNMAT 


DYNINT, DYNLST, DMMMAT, DMDMAT, DMSMAT, DMATEX, BLELPD 


None 


RESTRICTIONS: 


None 



NAME: 

PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


EIGER 

J 

To compress and link the rotor matrices to the bifilar 
analysis. I 


The degrees-of- freedom not utilized are eliminated from 
the rotor dynamic and aerodynamic matrices. The 
dynamic and aerodynamic damping and stiffness matrices 
are added together. Then, the 30 x 30 matrices are 
compressed to. K X K in CMPRSS. The final compressed 
matrices (3) are stored in labelled COMMON/INEIG for 
coupling with the bifilar analysis. The matrices can 
be printed out and/or punched out in cards if desired. 


CALL EIGER 
CMPRSS 


None 


I 

| 

1 


RESTRICTIONS: 


None 



PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To calculate the forces and moments needed to create 
the elastic matrix associated with the inboard half 
of a segment. 


The forces and moments are calculated using the ex- 
pressions derived in Reference 2. 


CALL ELI (I) 

I = Blade segment number. 
None 


None 


RESTRICTIONS: 


None 



NAME: 


ELO 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


To calculate the forces and moments needed to create 
the elastic matrix associated with the outboard half 
of a segment. 


The forces and moments are calculated using the ex- 
pressions derived in Reference 2, 


CALL ELO(I) 

I = Blade segment number. 


None 


None 


None 
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NAME: 


EXTEND 


PURPOSE: 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 



To find the values of a function at a set of points, 
given the function values at a different set of points. 


The value of the function at a new point is obtained 
by linear interpolation between the two closest old 
points. 


CALL EXTEND (F, R, M, RSTAR, N) 

F = On input, the old function values. 

On output, the new function values. 

R = The set of points at which F is defined on 
input. 

M = Number of points in F and R. 

RSTAR = The set of points at which F is defined on 
output. 

N = The number of points in RSTAR and F on output. 
None 

None 

N must be less than or equal to 400. 
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NAME: 

PURPOSE: 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


E159X 


To calculate the blade elastic torsional frequency. 


The torsional frequency is found from repeated calls 
to ROOTX which calculates the torsional mode shape. 
The frequency trials start at zero frequency and 
proceed in steps of 10 rad/sec up to a maximum of 
1990 rad/sec. After each trial, a check is made 
on the sign of the root mode shape. If a change 
in sign is found, then the frequency has been found. 
Three more iterations are performed to zero- in the 
frequency value (difference in 2 successive values 
of the root mode shape is within .0001). 


CALL E159X 


ROOTX 


If after 200 trials no sign change in the root mode 
shape has been found, then an error message is 
printed out as follows: 

'OUT OF RANGE' 

followed by the last frequency trial value and the 
root mode shape value. 


1. Frequency upper limit is 1990 rad/sec. 

2. Only the first elastic torsional frequency is 
found. 

3. Three trials allowed to zero-in the final 
frequency value. 



NAME: 

I 

j PURPOSE: 

! 

i 

METHOD: 


FILL 


To find the value/length of a second momentof area 
function over each of the blade segments, given the 
function/ length values over some other segment dis- 
tribution. 

I 

•• | 



The redistributed function G m is defined below for new 
segment RR m : 

= V r fl + Vr r k + VfVi 

G m F k-1 F k "k-l 

Where: 

= Old function values 
= Old blade segments 

= Radial positions of function from hinge (E) 

= New function values 
= New blade segments 


a 

r!1 


m 



USAGE: 



CALL FILL (RR, NS EG, F, RF, NF, E) 

RR « An array containing the blade segment lengths. 

NSEG = The number of segments in RR. 

F = On input, the function values over RF. 

On output, the function values over RR. 

RF = An array containing the segment lengths over 
which F is defined on input. 

NF = The number of segments in RF. 

E = The offset of the blade from the center of 

rotation. 

SUBROUTINES None 

CALLED: 

ERROR RETURNS: None 

RESTRICTIONS: NSEG must be less than or equal to 25. 
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NAME: 


FIXABS 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


; 


To calculate the fixed system absorber matrices 
in the bifilar analysis. 


This routine initializes the matrices to zero. Then, 
it defines needed parameters from the input vector, 

V, and proceeds to evaluate the mass damping and stiff- 
ness elements according to the equations from Reference 
3. 

Printout of final matrices is controlled by location 
■ 10 . 


CALL FIXABS 


None 


Number of fixed system absorber cannot be greater 
than 5. 



NAME: 


FIXSYS 


1 PURPOSE: 

! 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To calculate the fixed system modes matrices 
in the bifilar analysis. 


The fixed system modal mass, damping and stiffness 
elements are calculated according to the expressions 
in Reference 3. 

A printout of the matrices can be obtained if location 
10 is set to 1, 


CALL FIXSYS 


None 


None 


RESTRICTIONS: 


None 



NAME: 


FOLL 


PURPOSE: 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


To convert a blade function from value/length to an 
equivalent radial function. 


The input function is first converted into a cumulative 
function with values at the end points of the input 
segment distribution. The value then assigned to a 
radial point is derived by finding the cumulative 
values at the end points of the segment that the 
radial point represents and subtracting them. 


CALL FOLL (RAD, F, NF, RB, NB, E, RR) 

RAD = Segment lengths of input distribution. 

F = Input stepwise function. 

NF = Number of elements in RAD and F. 

RB = The radial points at which values are required. 

NB = The number of elements in RB. 

E = The blade offset. 

RR = The segment lengths to be associated with 
the radial points in RB. 


EXTEND 


None 


None 
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NAME: 

I 

1 

! PURPOSE: 

I 

j METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 



FORCER 


To solve for the forced response of the linear bi filar 
analysis. 


At first, this routine combines the generalized force 
cosine and sine vectors and prints out the final 
vector, FRC, of dimension 2*NF. Then, it inverts the 
hub impedance matrix, T, by a call to LINV2F and forms 
the fixed system generalized coordinates vector, XQ, 
of dimension 2*NF. The generalized coordinates vector 
of all the system degrees- of- freedom, GAMMA, of 
dimension 90, is calculated from a multiplication of 
the transfer matrix, TRANS F, and the XQ vector. The 
final vector GAMMA is used to calculate and print out 
the forced response amplitudes and phase angles of the 
fixed system absorbers and bifilar pendulums. 

The GAMMA vector is printed out if location 1497 is set 
to 1.0. 


CALL FORCER 


LINV2F 

(This is an "IMSL" package routine which must be 
supplied by the Army). 


None 


Refer to subroutine HUB IMP.. 

jj 

j 


1 





PURPOSE: 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


FREQUN 

To calculate the blade natural frequencies. 


The natural frequencies are the eigenvalues of the 
matrix formed by applying suitable boundary conditions 
to the matrix relating forces, etc., at the blade tip 
to the blade root. 

The eigenvalues are found by assuming an initial 
frequency and then increasing it by a predetermined 
increment until a change in sign of the determinant 
occurrs. Newton's method is then used to improve the 
answer until |F. + ^ - F. | < .001. 

The search for further natural frequencies continues 
until either the requested number has been found or 
the upper frequency limit is reached. 

The starting values for the frequency trials is .15ft 
where ft is the rotor speed in rad/sec. 

The frequency scan interval is ,10ft. 


CALL FREQUN 


PRODM, MIND 


1. 'EXCEEDED UPPER FREQUENCY LIMIT' { 

The program will search for the number of freq- 
uencies requested up to a frequency of 75 cycles/ 
rev. 1 

The program continues using the number of natural 
frequencies that have been found. 

jj 

2. 'DID NOT CONVERGE AFTER SIGN CHANGE - 100 TRIALS' 

Having located a change in sign of the determinant, 
the program was unable to improve the eigenvalue 
sufficiently to satisfy the convergency test in 
100 trials. The program continues using the last 
estimate. I 



3 


FAILED TO LOCATE SIGN CHANGE AFTER 500 TRIALS 


RESTRICTIONS: 



An eigenvalue was not detected in the frequency 
range 0.0 to 75 cycles/rev. The program continues 
with a frequency of 75 cycles/rev. 

4. 'AN EVEN NUMBER OF FREQUENCIES HAVE BEEN MISSED' 

Two frequencies differ by less than .0015Q - 
.00001 (n is the rotor speed). The program 
treats them as one frequency and continues. 


1. Frequencies cannot be calculated at zero rotor 
speed. It is suggested that rotor speed be 
lowered gradually to about 50 rad/sec. Computer 
time increases tremendously as rotor speed 
decreases since the scan interval is set at 

. loo. 

2, See notes on "ERROR RETURNS". 



NAME: 


FULL 


PURPOSE: 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


To find the values of an area function over each of 
the blade segments, given the function/length values 
over some other segment distribution. 


Same method as discussed for subroutine "FILL". 


CALL FULL (RR, NSEG, F, RF, NF, E) 

RR = An array containing the blade segment lengths. 
NSEG = The number of segments in RR. 


F = On input, the function values/ length over L. 

On output, the function values/length over RR. 


RF 

NF 

E 


= An array containing the segment lengths over 
which F is defined on input. 

= The number of segments in RF. 

= The offset of the blade from the center of 
rotation. 


j 


None 


None 

NSEG must be less than or equal to 25. 



NAME: 

1 

| PURPOSE: 
j METHOD: 

ii 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


GENFOR 


To calculate the generalized forces for the linear 
bi filar analysis. 


The generalized force vectors for the cosine (FQC) 
and sine component (FQS) are calculated from the 
force input vectors (locations 110-289) provided 
for the main and tail rotor hubs and 2 additional 
aircraft points and their appropriate mode shapes 
(locations 450-749 and 850-949). 

The force vectors of dimension NF (location 9) are 
passed through labelled COMMON/XFRDAT. 


CALL GENFOR 
None 


None 


RESTRICTIONS: 


None 



GMPRDD 


j NAME: 


| PURPOSE: 


METHOD: 


USAGE: 


SUBROUTINES 

CALLED: 


ERROR RETURNS: 


To multiply two general matrices to form a resultant 
general matrix. 


The M X L matrix B is premultiplied by the N X M 
matrix A and the result is stored in the N X L 
matrix R. 

RESULT: 


R = A X B 
NXL NXM MXL 

CALL GMPRDD (A, B, R, N, M, L) 

A = First input matrix. 

B = Second input matrix. 

R = Output matrix. 

N = Number of rows in A. 

M = Number of columns in A and rows in B. 

L = Number of columns in B. 


None 


None 


RESTRICTIONS: 


None 



NAME: 


HARMON 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To perform harmonic analysis of the time history 
solution. 


The solution acceleration vector calculated in INTEQ 
and passed through labelled C0MM0N/NLDAT2 is harmonic- 
ally analyzed once the convergence criterion is met 
in routine CONVER. 

The harmonics are stored in labelled COMMON/HARM. 


CALL HARMON (IC, Nl) 

IC = Rotor revolution number - location 1762 divided 
by 360. 

Nl = Total number of degrees- of- freedom (maximum 
is 72). 


None 


None 


RESTRICTIONS: 


The number of harmonics output is limited to 10. 



I NAME: 

I 

i PURPOSE: 

| 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


HUBIMP 

To develop the transfer and hub impedance matrices 
for the linear bifilar analysis. 


At first, the final mass (XMFC) and stiffness (XKFC) 
are combined to form matrix A (= XKFC -Wp 2 *XMFC), 
where Wp is the forcing frequency. Then, the damping 
matrix (XCFC) is used to form matrix B (= Wp*XCFC). 

If no fixed system absorbers and no linear bifilars 
are present in the system, the routine calculates 
the hub impedance matrix, E, from matrices A and B 
and then returns to MAINSV; otherwise, it proceeds 
to calculate the transfer matrix, TRANSF, after a 
call to LINV2F, and subsequently the hub impedance 
matrix E. (Additional information can be obtained 
in Ref. 3). 

The transfer and hub impedance matrices are passed 
through labelled COMMON/ XFRDAT. 

Throughout the matrix calculations performed, the 
resulting matrices can be printed out using the control 
switches in locations 10 and 15. 


CALL HUBIMP 


LINV2F 

(This is an " IMSL" package routine which must be 
supplied by the government) 


None 


1) Number of degrees-of- freedom of the final system 
matrices (XMFC, XKFC, XCFC) is limited to 60 


as 

fol 1 ows : 

1 

a. 

Fixed system modes 

16 j 

b. 

Fixed system absorbers 

5 ! 

c. 

Linear bifilars (5X3) 

15 

d. 

Rotor Modes 

24 | 

60 = Total | 



2) The transfer matrix (TRANSF) maximum dimensions 
are (90X32) where 90 represents 2 times the 
maximum d.o.f. of the fixed system absorber 
plus the bifilar and rotor modes and 32 is 2 
times the maximum number of fixed system d.o.f. 

3) The hub impedance matrix (E) has maximum dimen- 
sions of (32X32) obtained from 2 times the maximum 
fixed system d.o.f. 


\ 

i 

| 




NAME: 


INCOND 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To initialize quantities to zero in the bifilar analysis. 


The system mass, damping and stiffness matrices 
(XMFC, XCFC, XKFC respectively) of dimensions 60X60 
are set to zero. In addition, the matrices (10X75) 
used to store the harmonic results (CHARM and SHARM) 
are zeroed out. The input location 1498 controls 
the initialization of the bifilar pendulums, the 
rotor hub and the state variables displacements and 
velocities. 


CALL INCOND 


None 


None 


! 


} 

I 

I 


j 


RESTRICTIONS: 


None 



NAME 


INPUTV 



PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To provide a description of the input locations to 
the bifilar analysis and to read the input data. 


Locations 1 through 2200 are listed and discussed - 
all lines of code are of course commented out. Then, 
a call is made to LOADIT to read the input data. 


CALL INPUTV 


LOADIT 


None 


1 

RESTRICTIONS: None 


j 


; 


j 



i 

i 
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NAME: INTEG 


5 PURPOSE: To convert a function defined at unequally spaced 

} argument points, to an equal function defined at a 

set number of equally spaced points in a given interval 
and integrate it using the trapezoidal rule. 


METHOD: The length of the interval is divided by the number 

of divisions required to obtain the argument spacing. 
Function values at these subdivision points are ob- 
tained by linearly interpolating between the two near- 
est existing function values. The redistributed 
function is then passed to QTFG for integration. 


USAGE: 

CALL 

INTEG (F, R, NSEG, RBAR1 , RBAR2, N, XI) 


F 

- 

Array of function values. 


R 

= 

Array of argument values at which F is 
defined. 


NSEG 

= 

Number of blade segments 


RBAR1 

= 

Lower limit of the interval of integration. 


RBAR2 

= 

Upper limit of the interval of integration. 


N 

= 

Number of subdivisions to be used. 


XI 

= 

Integral . 

SUBROUTINES 

CALLED: 

QTFG 



ERROR RETURNS: 

None 




RESTRICTIONS: 


N must be less than or equal to 100. 



NAME: 

3 

! 

J PURPOSE: 
METHOD: 




To calculate the time history solution of the non- 
linear equations of motion. 


The combined mass matrix (XMT) and force vector (FT) 
of dimensions NT developed in COMBIN are used to 
obtain the time history response of the non-linear 
system. The procedure is discussed below. 

Given [xMt] {q* } = { Ft}, 


first, the matrix and vectors are partitioned as 
follows: 


kr A 

! xmt b 
1 


*1 


'fV 

j 

j 

XMT C 

1 

l I 
1 


• • 

q 2 


_ FT 2_ 

! 


where the square matrix, [xMT/\] , the acceleration 
vector part, {qj}, and the force vector part, {FT 3}, 
have dimensions (NT-ND) and include only the degrees - 
of- freedom associated with the fixed system and rotor j 
modes, while the_unity matrix, [i], the acceleration j 
vector part, ifq^} , and the force vector part, ; 

{FT2}, have dimensions ND and include the rest of the 

d.o.f. of the system. i 

i 

i 

Next, the accelerations are solved for as shown I 

below: ! 


i) 

"xmt a 

{-i} + 

’ xmt b 

{q 2 } - { FTj 

2) 

>MTc 

i q ii + 

I 

{q 2 } - { FT 2 

Solve 

for « 

j^q 2 | from equation 2) above. 


3) {^ 2 } = { FT z} - [ XMT c] {^} 

• • 


Then, the solution vector q^ is given by substituting 
3 ) into 1 ) above. 


4) ( 

-XMT fl ] - 

’ xmt b 

M)K} = 


l FT l} - 

' xmt b' 

Ki 1 
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USAGE: 


116 


which can be written as, 

5) [XXT] = {xi} 

The routine forms the matrix [XXT] and the vector 
{XT} (note that the maximum dimensions of (NT-ND) 
are 40). 

The acceleration vector, {q^} , is solved for in 
the IMSL routineLEQT2F. Subsequently, the acceler- 
ation vector, Jq 2 \ , is evaluated from expression 3) 
above. 1 J 

The analysis proceeds next to integrate the combined 
acceleration vector, {q’j , to obtain the velocity 
and displacement vectors using the expressions 
below: 

6 > Wt+At + 

7 > Wt+At = Mt + At { q }t 

where the time increment. At, is defined from 

8) At = Aij) (deg-loc 1761) , seconds 

(rpm-loc 7) *6 

The resulting velocity vector is loaded into the input 
vector, V, in locations 1740-1759 for the non-linear 
bifilars and in locations 1860-1939 for the remaining 
d.o.f. , while the displacement vector is loaded 
similarly into locations 1720-1739 and locations 1780- 
1859. The rotor hub velocities and displacements are 
also calculated by pre-multiplying the resulting 
vectors by the transfer matrix (input in locations 
450-549). The resulting velocity and displacement 
vectors are printed out at every azimuth position 
up to 30 degrees. The final step in INTEQ is to 
increment the azimuthal angle and return to NLBIF. 

The final results are passed through labelled COMMON/ 
INDAT for the input vector V and through labelled 
C0MM0N/NLDAT2 for the solution acceleration vector 


CALL INTEG (NT, ND) 

NT = Total number of d.o.f. (72 maximum). 






1 SUBROUTINES 
CALLED: 


ERROR RETURNS: 


ND = NT minus number of d.o.f. of fixed system 

(16 maximum) and rotor (24 maximum); maximum 
value of ND is 32. 


LEQT2F 

(This is an "IMSL" package routine which must be 
supplied by the Army). 


None 

The maximum order of the solution vector (as obtained 
from LEQT2F ) is 40. 


RESTRICTIONS: 



NAME: 

i PURPOSE: 

j 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


LINBIF 


To calculate the linear inplane bi filar matrices 


This routine defines needed parameters from the input 
vector, V, to start. Then, it initializes the matrices 
to zero and proceeds with the calculations of the 
mass, damping and stiffness elements according to 
the expressions published in Ref. 3. The matrices 
are of order 9X9 and consist of the following degrees- 
of- freedom: 

1. Fixed system longitudinal 

2. Fixed system lateral 

3. Fixed system vertical 

4. Fixed system roll 

5. Fixed system pitch 

6. Fixed system yaw 

7. Bi filar symmetric mode 

8. Bi filar cyclic (sine) mode 

9. Bi filar cyclic (cosine) mode 

A printout of the bi filar matrices can be obtained 
by setting location 1495 to 1.0. 


CALL LINBIF 


None 


None 


1. Number of bifilars in each kind is limited to 10. 

2. Number of different bifilar kinds (inplane and 
vertical) cannot be greater than 5. 



NAME: 


PURPOSE: 


METHOD: 


USAGE: 


SUBROUTINES 

CALLED: 


LOADIT 


To read a data card, and then after checking the 
characters, to store the data in the specified locations.! 


Each card in the input stream is read and printed out 
before any interpretation is attempted. Each charac- 
ter on the card is then checked for validity. Column 
one must be -, +, 0 or blank; column two must be 1 
through 5; columns three through six must be +, 0 - 9, 
or blank; and columns seven through sixty-six can be 
-, +, 1, 0 - 9, E or blank. 

A minus in column one indicates end of input data for 
that case. 

| 

Column two contains the number of values to be read in. j 

Columns three through six contain the input location 
at which to start storing the values. If the address 
is zero or blank, the next location is used. ] 

Columns seven through sixty-six contain the values to 
be stored. A format of 5E12.4 is assumed. 


CALL LOADIT(X,NFILE) 

X = An array into which the data is placed. 

In the rotor analysis, X = INPUT (which is 
equivalent to blank COMMON). 

In the bifilar analysis, X=V (which is 
equivalent to labelled COMMON/ INDAT . 

NFILE = 5. The read unit file for both rotor and 

bifilar analyses. 


None 


Any error results in the card that caused the error 
being ignored. 


ERROR RETURNS: 



RESTRICTIONS: *1) LOADIT is one of four routines which are computer 

dependent. Coding for both IBM and CDC computer 
systems is retained with appropriate Tines 
commented out. 

2) The maximum number of input locations is 8100 for 
the rotor analysis and 2200 for the bifilar 
analysis. 



j NAME: 

| PURPOSE: 

I 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


LVBIF 


To calculate the linear vertical bi filar matrices 


This routine defines needed parameters from the input 
vector, V, to start. Then, it initializes the matrices 
to zero and proceeds with the calculations of the 
mass, damping and stiffness elements according to 
the expressions published in Ref. 3. The matrices 
are of order 9X9 and consist of the following degrees- 
of- freedom: 

1. Fixed system longitudinal 

2. Fixed system lateral 

3. Fixed system vertical 

4. Fixed system roll 

5. Fixed system pitch 

6. Fixed system yaw 

7. Bi filar symmetric mode 

8. Bi filar cyclic (sine) mode 

9. Bi filar cyclic (cosine) mode 

A printout of the bi filar matrices can be obtained 
by setting location 1496 to 1.0. 


CALL LVBIF | 

1 

None i 


None j 

I 

•? 

1. Number of bifilars in each kind is limited to 10. j 

2. Number of different bifilar kinds (inplane and 
vertical) cannot be greater than 5. 

I 



NAME: 


MAINSV 


PURPOSE: To control the principal logic flaw of the bifilar 

analysis portion of the rotor/bifilar coupled program. 


METHOD: MAINSV first calls INPUTV which calls LOADIT to ] 

read the input data for the case at hand. Then, the 
bifilar matrices and other quantities are set to zero 
in INCOND. The contributions of the fixed system 
modes, the rotor (if used), the fixed system absorbers, j 
and the linear inplane and vertical bifilar pendulums j 
to the final system mass, damping and stiffness matrices j 
are added in SYSCTL. The program then proceeds to 
calculate the hub impedance and transfer matrices in 
HUBIMP, and the generalized forces and the forced | 

response in CMPUTE for the linear analysis case. ; 

If non-linear inplane bifilars are to be analyzed, j 
the program bypasses HUBIMP and CMPUTE and instead it j 
activates NLBIF to calculate the time-history response 
of the system. I 


USAGE: CALL MAINSV (I927SW, NCASE) 



I927SW 

= t o 

Include rotor contributions. 

1 



0 

| 

Do not include rotor contributions. 

‘ 

NCASE 

= 1 

j 

Input vector, V(2200), is initialized 
to zero and the bifilar analysis title ; 
(second one) is read in. 



>2 . 

\ 

Input vector is not set to zero and 
the second title card is not read. 

SUBROUTINES 

CALLED: 

INPUTV, 

INCOND, 

SYSCTL, HUBIMP, CMPUTE, NLBIF 

| 

ERROR RETURNS: 

None 


i 

RESTRICTIONS: 

None 


i 



NAME: 


MATE I 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


H Q'e:;; 


To calculate the elastic matrix associated with the 
inboard half of a rotor blade segment. 


The 13 X 13 matrix is calculated using the expressions 
derived in Reference 2 and the forces and moments 
calculated in ELI. 

CALL MATEI(I) | 

I = Blade segment number. 

MIND, GMPRDD | 


None 

None 




) 

i 

l 

! 


| 

j 



| NAME: 

I PURPOSE: 

I 

| 

| METHOD: 

i' 4 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


MATEO 


To calculate the elastic matrix associated with the 
outboard half of a rotor blade segment. 


The 13 X 13 matrix is calculated using the expressions 
derived in Reference 2 and the forces and moments 
calculated in ELO. 


CALL MATEO(I) 

I = Blade segment number. 
MIND, GMPRDD 


None 


RESTRICTIONS: 


None 




ft*®**™*^ 

NAME: 

j PURPOSE: 
METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


MATF I 

:| 

To calculate the mass transfer matrix which relates 
the blade forces across a concentrated mass. 

I 

| 

The 13 X 13 matrix is calculated using the expressions 
derived in Reference 2. 

.{ 

CALL MATF(I,W) j 

I = Blade segment number. 

| 

W = Rotational speed. ] 

j 

" 

None I 


None 


RESTRICTIONS: 


None 



PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To calculate the transformation matrix associated 
with a discontinuity in blade twist. 


The 13 X 13 matrix is calculated using the expression 
derived in Reference 2 and represents the change in 
twist from segment I to segment 1-1. 


CALL MATR(I) 

I = Blade segment number. 


None 


None 


RESTRICTIONS: 


None 



NAME: 

j PURPOSE: 
METHOD: 


USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 


MIND 

1 

: 

! 

To invert a matrix and calculate its determinant. 


The routine is an adaption of the routine MINV from 
the IBM Scientific Subroutine Package, which uses the 
standard Gauss- Jordan reduction to obtain the inverse. 

The value of the determinant is obtained simultaneous- 
ly as the product of the pivot. j 

To avoid possible overflow conditions, the value of 
the determinant is expressed in the form D X 10* , where i 
D is between 1.0 and 10.0, by repeated calls to OVUN. 


CALL MIND (A, N, D, L, M, IE) 

A - The matrix to be inverted; it is destroyed 
during the computation. 

N = The dimensions of A. 

D = Determinant value. 

L = Work vector of length N. 

M = Work vector of length N. 

IE = Power of 10 to be associated with D. 


OVUN ; 

I 

A determinant of zero indicates a singular matrix. 


RESTRICTIONS: 


None 



PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To perform initial calculations on the input data for 
the blade frequency and mode shape calculations. 


The input segment lengths are used to construct an 
array of radii to the segment centers, and the local 
twist distribution is used to construct the twist dis- 
continuities. 


CALL MI SC 


None 


The sum of the segment lengths is checked against the 
blade radius; a difference of more than 10% produces 
the following warning: 

' INCOMPATABILITY BETWEEN RADIUS AND SUM OF SEGMENTS'. 
The program continues with the input values. 



NAME: 


MODES 


PURPOSE: 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 



To control the calculations of the rotor blade fully- 
coupled frequencies and mode shapes. 


The blade is assumed to consist of a number of span- 
wise segments with the inertial loading on each seg- 
ment concentrated at the center. Each segment is 
then divided into two parts: one inboard and one 

outboard of the concentrated mass. The half segments 
are treated as weightless, with the concentrated mass 
being located at the junction between them. The elas- 
tic properties are assumed to be constant within the 
inboard and outboard halves of the segment. The built- 
in blade twist is incorporated in the model by per- 
mitting angle changes at the junction between segments. 

Relationships between blade forces, moments, and 
elastic deformations at the tip and at the root of 
the blade can be obtained by considering the changes 
in these variables over the blade segments. Elastic 
matrices are used to give the change across the in- 
board and outboard sections of the segments. Trans- 
formation matrices account for the change due to an 
abrupt change in twist at the segment junctions, and 
mass transfer matrices give the change from a position 
just outboard to just inboard of a concentrated mass 
at the center of a segment. The product of these 
matrices for all segments will relate the variables 
at the tip to those at the root. 

By applying suitable boundary conditions and iterating 
on frequency, the blade natural frequencies and hence 
the mode shapes can be derived. 


CALL MODES 

MISC, PINT, FREQUN, MSHAPE, POUT, ORTHOG 

I 

| 


None 


RESTRICTIONS: 


None 



NAME: 

PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


130 | 



MSHAPE 


To calculate the blade mode shapes and their first 
and second derivatives. 


This routine first determines whether -the mode shape 
will be predominantely flatwise, edgewise, or torsional 
by examining the size of the determinant of the matrix 
relating the tip properties to the root properties at a 
frequency just less than the natural frequency. 

This knowledge is used in setting up a series of 
simultaneous equations whose solutions yield the 
required mode shape on substitution into expressions 
derived in Reference 2. 


CALL MSHAPE 


PRODM, MIND, OVUN, SIMLIN, MATF, ELO, MATEO, MATR, ELI, 
MATE I, GMPRDD 


None 


None 




NAME: 


NLBIF 


PURPOSE: 

METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


To control the non-linear bifilar analysis calculations. 


The non-linear bifilar analysis is performed if the 
control switch in location 18 is set to 1 and the 
number of non-linear bifilars in location 1763 is 
greater than zero. 

At first, the routine calculates the total number of 
degrees-of- freedom and proceeds with the calculations 
if it is equal or less than 72. Then, it calculates 
the time history response by calling, in sequence, 
the routines RHS, BIFILR, BIFEXP, COMBIN, INTEQ for 
each rotor revolution until either the maximum azimuth 
angle (input location 1762) is reached or the conver- 
gence criterion is met, as specified in CONVER. In 
either case, the analysis proceeds to analyze the 
harmonic response of the bifilar pendulums, the rotor 
hub and the aircraft stations in HARMON. The results 
are printed out in the routine OUT. The initial values 
of the bifilar, hub and state variables displacements 
and velocities are then listed to be used as starting 
values for the next case. 


CALL NLBIF 

RHS, BIFILR, BIFEXP, COMBIN, INTEQ, HARMON, CONVER, OUT 


If the number of degrees-of- freedom is greater than 72, 
then the non-linear bifilar analysis is not activated 
and the following message is printed out: 

'TIME HISTORY SOLUTION WAS NOT PERFORMED SINCE TOTAL 
NUMBER OF D.O.F. = , I.E. > 72' . 


1. Number of non-linear inplane bifilars is limited 
to 12. 

2. Total number of d.o.f. is 72 (60 from linear 
analysis plus 12 for non-linear analysis). 



NAME: 


ORTHOG 


[ PURPOSE: 
METHOD: 


USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURN: 


To test the orthogonality of the blade mode shapes. 

The orthogonality relation is shown below. 

2 k m k ♦i.k *i,k 

ORTHO , j) = ■■ 

V z k m k (<l> i,k + *j,k ) 
where m^ = Mass of segment k. 

J. L. 

<f>.j ^ = i— mode shape for segment k. 

4>j s k = j~ mode shape for segment k. 

CALL ORTHOG 

None 

None 


RESTRICTIONS: 


None 



PURPOSE: 

To print out the harmonic results for the non-linear 
bi filar analysis. 

| 

METHOD: 

At first, this routine calculates some parameters ] 

from the input vector V. Then, it uses the harmonic 
analysis results from HARMON to calculate the 
amplitudes and phase angles of the harmonic response 
of the non-linear bifilar pendulums, the rotor head 
and the aircraft stations. Appropriate multiplications 
with the mode shapes input for the fixed system and 
the aircraft stations are performed. 

i 

USAGE: 

CALL OUT (NTOT, Nl) 

NTOT = Total number of system d.o.f. 

excluding the number of non-linear bifilars. 

Nl = Total number of system d.o.f. 

SUBROUTINES 

CALLED: 

None 

; 

| 

ERROR RETURNS: 

1 

None 

RESTRICTIONS: 

1. Maximum number of output harmonics is 10. 

| 

2. Maximum number of aircraft stations whose 


response is harmonically analyzed is 4. 



I NAME: 

! 

I PURPOSE: 

j 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 


OUTPUT 

To calculate and print the aircraft and rotor hub 
forced response. 


This routine uses the results from FORCER (specifically 
vector XQ) to calculate the forced response at specific 
aircraft stations (up to 4) where mode shapes are 
loaded in locations 1000 through 1399 and at the rotor 
head using the mode shape in locations 450 through 
549. 


CALL OUTPUT 


None 


None 


The response of the aircraft can be evaluated at a 
maximum of 4 stations. 



PURPOSE: 


METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To express A in the form B x 10*, where 
1.0<B<10.0. 


If A is greater than 10.0, it is continually divided 
by 10.0 until it is less than 10.0. If A is less 
than 1.0, it is continually multiplied by 10.0 until 
it is greater than 1.0. 


CALL OVUN (A, B, I) 

A = Input number to be transformed. 

B = Part of B which is greater than or equal to 1.0 
and less than or equal to 10.0. 

I =* The power of 10 which is associated with B to 
form A. 


None 


None 


RESTRICTIONS: 


None 



NAME: 

i PURPOSE: 

METHOD: 
Column 1 

IT 


w = 


PFMULT 


To generate a point field matrix, relating the steady 
deflections of the center of a blade segment to those 
of the previous segment center. 


The point matrix is defined as W(5X5): 


2 

L/A 

2/A-l 


3 

-L*J/2 

-J 


-2B/(A*J) 1/A 


0 

0 


0 

0 

Where: 


5 Row 

C*J*L/2 + L*y*B ~T 
A*J*L*(l/A+0.5)/3 C*A* J* ( B+2 )+2*y*B 2 

- C* ( B+6 ) - 2*y*B/ ( A* J ) 3 
-L*D 4 

1 5 


•L 1 2 3 4 5 6 7 8 9 10 11 *J/6 


L*(l/A+2)/3 

1 

0 


1) D = Thrust or drag derivative/unit length 

between segment centers. 

2) El = Flatwise or edgewise segment stiffness. 

3) L = Segment length. 

4) M = Segment mass. 

5) R = Segment radial position. 

6) r = Blade lag or coning angle. 

7) ft = Rotor speed. 

8) A = [l-L 2 *ft 2 /(2*EI)] * £M.R i 


9) B = 1/A-l 

10) C = L 2 *D/12 

11) J = 


L/(EI*A) 




For inplane deflections, a point mass matrix is 
added to the fourth row whose elements then become 

W (4,1) = W(4, I)-M*Q 2 * Jw ( 1 , 1 ) +R*y*W ( 5,1 )j 

Conditions at the next center are obtained by pre- 
multiplying W by the matrix representing conditions 
at the previous center. 


USAGE: 


CALL PFMULT (PFIN, PFOUT, LS, AS, EIS, DS, M0MG2S , 
GAMOS, RS, ISWTCH) 


PFIN = 5x5 matrix corresponding to previous center. 
PFOUT = 


LS 

AS 


5x5 matrix corresponding to next outboard 
center. 

Distance between centers. 

N 


l-L* 2 ft 2 /(2*EI)j 


* Z M i R i 


Flatwise or edgewise second moment of area 
between centers. 

Thrust or drag derivative/unit length 
between centers. 


EIS 
DS 

M0MG2S = M*Q^ (centrifugal force). 


GAMOS = y - calculated lag or coning angle. 

RS = Distance of next center from the blade root. 

ISWTCH = 0 for inplane deflections. 

= 1 for out-of- plane deflections. 


I 


SUBROUTINES 

GMPRDD 

CALLED: 


ERROR RETURNS: 

None 

RESTRICTIONS: 

None 
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PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


Given 3 points on the curve y=f(x), choose the best 
2 to use for an estimate of the solution of y=f(x). 


The 3 points are rearranged such that 

*1 “ *1 > *2 " x 2 - y 3 ' x 3 

if y 2 - x 2 — then P oints 1 anc * 2 in the 
c c rearranged order are used. 

If y 2 - x 2 > 0, then points 2 and 3 in the 
L c rearranged order are used. 

CALL PICK (TT, CT) 

TT = Array of argument values. 

CT = Array of function values. 

None 

None 


RESTRICTIONS: 


None 



NAME: 


PINT 


PURPOSE: To print the input data used in the blade natural 

frequency and mode shape calculations. 


METHOD: If the print option is set to 4, the input is printed. 

Otherwise, control is returned to MODES with no 
printing. 


USAGE: CALL PINT 

SUBROUTINES None j 

CALLED: j 


ERROR RETURNS: None 


RESTRICTIONS: None 




NAME: 

POUT 

1 

PURPOSE: 

To print the output of the blade mode subsegment of the 
program. 

METHOD: 

This routine 
program, but 

is not used by the rotor stability 
has been included for completeness. 

| 

USAGE: 

CALL POUT 

I 

j 

SUBROUTINES 

CALLED: 

None 

1 

I 

i 

1 

j 

ERROR RETURNS: 

None 

j 

1 

RESTRICTIONS: 

None 

j 


, 
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NAME: 


PRELIM 


PURPOSE: 

METHOD: 



To prepare the data for use by the rotor aeroelastic 
analysis. 


The two main tasks of PRELIM are to standardize the 
data tables so that all the data tables are defined 
over the same segment lengths or at the same radial 
points, and to control the calculation of other quan- 
tities needed throughout the remainder of the program. 

Transfer of data between the program segments is 
achieved through 3 COMMON blocks. Blank COMMON is 
used primarily to store input data. Quantities cal- 
culated by PRELIM from the input are stored in labelled 
COMMON/DYNINP. Printout control is maintained by 
switches set in labelled COMMON /PRNTSW. All other 
labelled COMMON blocks are used for the transfer of 
data within the PRELIM segment. j 

i 

The following is a more detailed breakdown of the 
tasks performed. j 

1. The input data are read into blank COMMON via 

subroutine LOADIT and immediately stored on a 
temporary work file (Unit 11). Prior to each j 

subsequent case in the run, this temporary 

file is read into blank COMMON; thus, only those 
variables which differ from the previous case 
need be input. 

I 

It should be noted that PRELIM is the third 
computer dependent routine due to different 
requirements for reading input data. Coding 
for both IBM and CDC systems is retained in ! 

the program with the appropriate lines commented 
out. I 

i 

2. The coupling matrix terms needed for the bifilar 
analysis portion are initialized to zero to 
start. For subsequent cases, the coupling 
matrices are retained and used again if the 
rotor is not changed or recalculated if a new 
rotor is employed. 
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3. Specific input locations are set for the coupling 
with the bifilar analysis portion. 

4. Control switches are set and, where applicable, 
inputs are converted from generally accepted 
engineering units to standard units for program 
calculations. 

I I 

5. The input blade segments, which form the basis of 
all radial distributions used by the program, 
are adjusted to include a 0.1 segment to accomo- 
date pitch horn effects. They are then used 

to set up a radius vector representing the 
distance of the mid- points of each segment 
from the center of rotation. The offset is 
automatically included if it exists. 

6. All input tables are extended or redefined over 
the distributions obtained in 3 above. 

7. Pitch horn effects are added to the blade center 
of gravity and blade weight tables. 

| 

8. The thrust for the input pitch angle is cal- 
culated for hover or vertical flight. 

9. An input vector used by MODES to calculate mode 
shapes and frequencies is set up. In multiple 
cases, if this input has not changed from the 
previous case, then no call is made to MODES 
and the program uses the values obtained in 
the previous case. 

10. Pitch-lag and pitch- flap coupling effects are j 

calculated using the blade bending modal compo- j 
nents at the pitch horn. j 

11. Blade lag damper coupling terms are calculated 
from the input blade damper geometry. 

12. The aerodynamic derivatives and steady deflections 
are calculated in hover. 

13. Control system inputs are checked to eliminate 
contradictions. 

1 

14. Blade torsional properties are calculated I 

including cross-beam blade characteristics. j 



15. Selected portions of the input and calculated 
data are printed out by subroutine PROUT. 



USAGE: 

CALL PRELIM (NCASE, NMODE) 


NCASE = 0 Input vector, INPUT (8100), is initialized 
to zero and the rotor analysis title 
(first one) is read in. 


> 1 Input vector is not set to zero, original 
input vector is read from Unit 11 and 
the first title card is not read in, j 


NMODE = Switch for calculating mode shapes and 

frequencies. 

SUBROUTINES 

CALLED: 

] 

LOADIT, SORTAB, EXTEND, FILL, SECAER, QTFG, PICK, 

MODES, REMOVE, STDEFL, FOLL, FULL, E159X, PROUT j 

j 

ERROR RETURNS: 

None 

RESTRICTIONS: 

None 


! 
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NAME: 


PRODM 


I 

> 

i 

| PURPOSE: 


METHOD: 


USAGE: 


SUBROUTINES 

CALLED: 


ERROR RETURNS: 


RESTRICTIONS: 



To calculate the matrix relating blade forces, moments 
and elastic deformations at the tip to the blade root 
and to apply the boundary conditions. 


Since the elastic matrices and twist transformation 
matrices are independent of the frequency, they are 
calculated once for each segment and stored on a work 
file. For the given frequency, the mass transfer 
matrices are calculated and combined with the matrices 
on the work file to produce the tip-to-root relation 
matrix. Boundary conditions are applied, and the re- 
sulting 8X8 matrix is returned. 

Some data are stored in a temporary work file (Unit 8). 


CALL PRODM (A, NEQ) 

A = The returned 8 X 8 matrix 
NEQ = The dimensions of A. 


ELO, MATEO, MATR, ELI, MATE I, GMPRDD, MATF 


None 

None 



NAME: 

| 

i PURPOSE: 
METHOD: 


USAGE: 


SUBROUTINES 

CALLED: 


ERROR RETURNS: 


PROUT j 

To print selected input and calculated quantities. 

All significant input and calculated quantities are 
printed out depending on the setting of the print 

switch. j 

1 

! 

CALL PROUT (NBLD, NSEG, IE, IF, GJ, T, BETOD, GAMOD) 

j 

NBLD = Number of blades. 

NSEG = Number of blade segments. 

IE = Blade edgewise second moment of area. 

i 

IF = Blade flatwise second moment of area. j 

GJ = Blade torsional stiffness. 

T = Blade thrust. j 

BETOD = Blade steady coning angle. j 

j 

GAMOD = Blade steady lag angle. j 

;! 

SKIPLN j 

.] 

None 


| 


RESTRICTIONS: 


None 



NAME: 

QTFG 

PURPOSE: 

To compute the vector of integral values for a given 
table of argument and function values by the 
trapezoidal rule. 

METHOD: 

This routine is a copy of QTFG, and is obtained from 
IBM System/360 Scientific Subroutine Package. 

USAGE: 

CALL QTFG (X, Y, Z, NDIM) 

X = The input vector of argument values. 

Y = The input vector of function values. 

Z = The resulting vector of integral values. 

NDIM = The dimensions of X, Y, and Z. 

SUBROUTINES 

CALLED: 

None 

ERROR RETURNS: 

None 


RESTRICTIONS: 


None 



NAME: 


REMOVE 


PURPOSE: 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To delete a blade flatwise or edgewise mode from an 
array of mode shapes and frequencies. 


The characteristics of each mode are examined until 
the required mode is found. Then, the mode shape and 
its frequency are deleted and the remaining modes 
are compressed in the array. 


CALL REMOVE (SHAPE, MODE, TYPE, N) 

SHAPE = An array containing the type of each mode. 

MODE = Number of modes to be examined. 

TYPE = The type of mode to be removed. 

N = The occurrence of type which is removed, 

i.e. , the 2nd flatwise. 

None 


A message is printed if the requested mode could not 
be deleted. 


RESTRICTIONS: 


None 



NAME: 

PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


RHS 


To shift the damping and stiffness contributions to the 
right-hand-side of the non-linear bifilar equations of 
motion. f 


The stiffness and damping matrices obtained at the 
completion of the calling sequence in SYSCTL are post- 
multiplied respectively by the initial values of the 
state variables displacements and velocities. The 
resulting vector, FRHS, is shifted to the r.h.s. by 
a change in sign and transferred through labelled 
C0MM0N/NLDAT2. It is dimensioned NRHS (see below). 


CALL RHS (NRHS) 

NRHS = Total number of degrees-of-freedom not 
including number of non-linear bifilars 
(maximum values is 60). 

j 

None 


None i 

J 

! 

| 


RESTRICTIONS: 


None 



NAME: 


ROOTX 


PURPOSE: 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 




To calculate the blade elastic torsional mode shape. 


The blade torsional mass moment of inertia and stiff- 
ness distributions are used to calculate the torsional 
mode shape using a given frequency trial and the 
input value of rotor speed. 

The frequency and mode shape trials are printed out 
if location 119 is 6. 


CALL ROOTX (W2, XHI, 02, YTI, YTK) 

2 

W2 = Square of frequency trial, (rad/sec) 

XHI = Calculated blade torsional mode shape. 


non-dimensional 

02 = Square of rotor 

YTI = Blade torsional 

in- lb- sec 2 

YTK = Blade torsional 

None 


2 

speed, (rad/sec) 
mass moment of inertia, 

stiffness, lb-in 2 /in 






j 


None 


None 


I 


i 

; 







NAME': * ' SECAER 


j PURPOSE: To calculate blade section coefficients and derivatives. 


METHOD: Lift, drag, and pitching moment coefficients are input 

against angle of attack for various Mach numbers. For 
a given angle of attack and Mach number, the required 
coefficients and their derivatives with respect to 
Mach number and angle of attack are obtained from 
the corresponding table by linear interpolation. 

The routine assumes airfoil data to be non-symmetric 
if the first input angle of attack value is negative. 

For angles of attack greater than 30°, a Mach number 
of 0.0001 is assumed. 


USAGE: CALL SECAER (ALPHA, AMACH, I, M, FF, DFF, DM) 

ALPHA = Angle of attack (radians). 

AMACH = Mach number. 

I = Not used by this program. 

RFM = 1 for lift coefficient. 

= 2 for drag coefficient. 

= 3 for pitching moment coefficient. 

FF = Returned coefficient. 

DFF = Derivative with respect to angle of attack. 

DM = Derivative with respect to Mach number. 

SUBROUTINES BLIN4 

CALLED: 


ERROR RETURNS: Any error termination from BLIN4 stops execution with 

the message 'TROUBLE IN BLIN4' , followed by 6 numbers 
which represent angle of attack in radians, angle of 
attack in degrees, input Mach number, Mach number used, 
the error switch L from BLIN4 and the switch M 
respectively. See BLIN4 for description of errors. 


RESTRICTIONS: See locations 1850-4548 of rotor stability input 

■ • - — descri ption and BL1N4. - - ~ — 
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SHAKIT 


NAME: 

PURPOSE: 

METHOD: 


USAGE: 

SUBPROGRAMS 

CALLED: 

ERROR RETURNS: 


To control the principal logic flow of the rotor/bifilar 
coupled program. 


MAIN first calls subroutine PRELIM, which controls 
the input, conversion and adjustment of the blade 
data. Then, the rotor blade dynamic and aerodynamic 
matrices are obtained from calls to DYNMAT and AERMAT 
respectively. The dynamic and aerodynamic matrices are 
cleaned out, added together, compressed and stored for 
coupling with the bifilar analysis portion in EIGER. 

The bifilar analysis is subsequently executed by call- 
ing MAINSV. 

After PRELIM is called, the rotor analysis calculations j 
are bypassed if input location 110 is zero or -1. 

It should be added that SHAKIT is the fourth computer 
dependent routine. For CDC use, the first line of 
code i s 

"PROGRAM SHAKIT (INPUT, OUTPUT, TAPE1 , TAPE2, TAPE3)" . j 

For IBM use, this card is not needed; thus, the first 
line contains blank COMMON input data. j 


Program SHAKIT is never referenced in a CALL state- 
ment. 


PRELIM, DYNMAT, AERMAT, EIGER, MAINSV 


None 


RESTRICTIONS: 


None 



NAME: 

! PURPOSE: 

| 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 



SIMLIN 



To solve a system of linear simultaneous equations. 


The solution is obtained by elimination using largest 
pivotal division. The determinant of the coefficient 
matrix is produced during the process. 


CALL SIMLIN (A, B, X, N, NX, C) 

A = Coefficient matrix (destroyed during compu- 
tations) of dimensions XN by XN. 

B = Vector of right-hand value, destroyed during 
the computations. On return, B(l) contains 
the determinant of the coefficient matrix. 

X = Solution vector. 

N = The number of rows in A. 

NX = Dimensions of matrix A. 

C = Work vector. 


None 

! 

j 


t 

None 

■j 

I 


RESTRICTIONS: 


None 



NAME: 


SKIPLN 


PURPOSE: 

METHOD: 

USAGE: 

SUBROUTINES 

CALLED: 

ERROR RETURNS: 


To skip a given number of lines during printing. 


In an effort to prevent the printout of output tables 
going over a page boundary, SKIPLN is used to center 
tables within a page or part of a page. 


CALL SKIPLN(J) 

J = Number of lines to be skipped. 


None 


None 


RESTRICTIONS: 


None 



j NAME: 

| 

| PURPOSE: 

| 

METHOD: 

USAGE: 


SUBROUTINES 

CALLED: 

ERROR RETURNS: 
RESTRICTIONS: 
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SORTAB 


To unscramble an array of alternating X and Y points 
into an array of X points and an array of Y points. 


Even subscripted values in the input array are moved 
to another array. Both arrays are then compressed. 


CALL SORTAB (G, M, KSTAR, NSTAR) 

G = On input, G contains the alternating X and 

Y values. 

On output, G contains only X values. 

M = The number of points in G on input. 

KSTAR = The unscrambled Y points. 

NSTAR = The number of points in G and KSTAR on 

output. 


None 


None 


None 








NAME 


STDEFL 


PURPOSE: 

METHOD: 

USAGE: 


To calculate the inplane and out-of-plane blade steady 
deflections and slopes. 


Blade deflections are calculated by successively 
calculating the deflection of each segment relative 
to the previous one. This is done by constructing a 
point matrix from the mass, stiffness, and geometric 
properties of the segment to relate the displacement 
of the center of a segment relative to the previous 
center. The process is started by assuming suitable 
conditions at the root of the blade. A correction is 
made to these deflections if the pitch bearing does 
not follow the blade root slope. 


CALL STDEFL (QEO, XEO, E, RB, MB, NB, RR, D, IE, IF, 
TH, OMEGA, EB, GAMOT, KGAMA, GAMO, KASE, IEF, ROTDEF) 


QEO 


Returned inplane or out-of-plane steady 
deflections. 

XEO 

= 

Returned slope for inplane or out-of-plane 
deflections. 

E 

= 

Blade offset. 

RB 

= 

Array of distances of the blade segment 
centers from the center of rotation. 

MB 

= 

Mass of each segment. 

NB 

= 

Number of radii in RB. 

RR 

= 

Blade segment lengths. 

D 

= 

Drag derivative for inplane deflections or 
thrust derivative for out-of-plane deflect- 
ions. 

IE 

= 

Blade edgewise second moment of area. 

IF 

= 

Blade flatwise second moment of area. 

TH 

= 

Blade twist. 

OMEGA 

= 

Rotational speed. 



EB 

GAMOT = 

KGAMA = 

GAMO = 

KASE = 
IEF 

ROTDEF = 

SUBROUTINES PFMULT 

CALLED: 

ERROR RETURNS: None 



Blade Young's modulus. ] 

Calculated lag angle for inplane deflections 
or calculated coning angle for out-of-plane 
deflections, j 

Blade lag hinge spring constant for inplane 
deflections or blade flapping hinge 
spring constant for out-of-plane deflections. 

Blade prelag angle for inplane deflections 
or blade precone angle for out-of-plane 
deflections. 

Blade pitch input control (derived from 
input location 115). 

0 for inplane deflections. 

1 for out- of- plane deflections. 

Rotor definition (input location 114). 


RESTRICTIONS: 


None 



NAME: 


SYSCTL 


PURPOSE: To include the contributions of the fixed system 

modes, rotor (if used), the fixed system absorbers 
and the linear inplane and vertical bifilars to the 
final system mass, damping and stiffness matrices 
in the bifilar analysis. 


METHOD: At first, SYSCTL calls FIXSYS to obtain the fixed 

system modes matrices. Then, if rotor contributions 
are desired, it restructures the rotor matrices 
and adds the rotor elements to the fixed system 
matrices using ADDOFR. Subsequently, it includes 
contributions due to fixed system absorbers from 
FIXABS, linear inplane bifilars from LINBIF and 
finally linear vertical bifilars from LVBIF. In 
all cases, the matrices are increased systematically 
by repeated calls to ADDOFR. 

Printout of the matrices is governed by the switches 
in locations 17 and 1490 through 1496. j 

j 

j 

USAGE: CALL SYSCTL i 


SUBROUTINES FIXSYS, ADDOFR, FIXABS, LINBIF, LVBIF 

CALLED: 


ERROR RETURNS: None 


RESTRICTIONS: None 
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FIGURE 2. Output Format - Rotor Blade Input Data 
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FIGURE 3. Output Format - Rotor Blade Characteristics 
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FIGURE 3. Continued. 
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-EIGURE..3 — Concluded. 
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FIGURE 5. Concluded 
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FIGURE 6. Concluded. 
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FIGURE 7. Output Format - Bifilar Nonlinear Analysis Matrices 
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FIGURE 7. Continued 
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FIGURE 8. Concluded. 
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Figure 9, Bifilar Analysis Segmentation Structure. 
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Fig ure 15. Flo w C hart for th e Bifilar Analysis. _ j 
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16c. NLBIF Flow Chart. 



Figure 16. Flow Charts for Subroutines SYSCTL, CMPUTE, NLBIF. 
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Figure 17. COMMON Blocks 
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Figure 17. Concluded. 
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...RHEUaJL. PHEU»H _PHEC2»II_ PHF(2,I).. PHEC3»II — PHFC 3*1) — PHEC4,I) : „PHF(4,I) 





306.45000 -0.00254 0.01581 0.00968 0.02471 
314.15000 -0.00255 0.01584 0.00968 0.02485 
320.00000 -0.00255 0.01585 0.00968 0.02486 



5 

i 



225 






MODEt 2 ) a .209446 LB SEC-SQ/IN 





FIMA l — STIFFNESS... MATRIX 
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_ FINAL ... ...DAMPING MATRIX. 






.FINAL J1ASS MATRIX 
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.0 .0 .0 .0 -.960889D-01 .6 .757781' .0 “ -2" 45428 

18203.9 .0 -17.7265 -172.099 .0 2230.99 -6.14894 
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TOTAL NO. OF DEGREES-OF- FREEDOM IS 26 





N*M*R*R*i l,..+W*M ).#*?„ ...... .1M.656.. 





GENERALIZED FORCES . ORDER IS 
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VERTICAL BIFILAR PENDU LUM( S ) 



£35034D-„0l....-.213538Dr01_, .454660 ^137.744. 

51G358D-02 .326855D-02 .338881D-01 -30.2823 

122553D-02 .391058D-02 .144276 116.499 




SECOND A.C. STATION DISPLACEMENT IN 6 











TOTAL NO. OF DEGREES-OF-FREEDOM IS 25 






—VERTICAL BI FILAR PARAMETERS 






JMLJL.-.S. .29 »NE=„.SI»KRQTQR =L.l 2* NFABS=_.l. KINDS OF... .INPL- .BIF.S_ 0 » KINDS. OF VERT... BIF,.= .-.-l,NUMB. OF N.L- INPU~BIF= 
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THE CONVERGENCE FACTOR TO G = 30.2 
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... VIBRATION LEVELS^AT .... A/C — LOCATIONS. 
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D.l Inplane Bifilar Equations of Motion 

The inplane bifilar rotor head absorber is a single degree- of- freedom 
system. It can be represented as a one degree- of- freedom pendulum as shown 
in Figure D.l-1. The rotating coordinate system shown here is located on 
the main rotor hub. The equivalent pendulum arm, r, and the distance from 
the center of rotation, R, are related by equation (1) (see Figure D.l-2) , 
often referred to as tuning equation. The location of the c.g. relative to 
the inertia frame system is given by equation (2) shown in Figure D.l-2. 
Substituting equation (2) into the Lagrange's equations, the non-linear 
equations of motion are derived, equations (3)-(6) , and shown in Figure 
D.l-3. 

To obtain the linear set of equations of motion, equations (7)- (10) 
(see Figure D.l-4), small angle assumptions were made, e.g. siny^ ~ yk, 
and cosyk ~ 1. Also, the second order terms were neglected. For identical 
bifilars, further simplification of these equations can be made by trans- 
ferring the rotating bifilar coordinate, y^, into the fixed system coordi- 
nates. The equation for the coordinate transformation used in shown in 
Figure D.l-5, equation (11). Thus the transformed equations of motion are 
derived, equations (12) - (17) shown in Figure D.l-5. 
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FIGURE 0.1-3: NON-LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR 


281 



X Equation 


(%+w T )y + + %u)* / 



FIGURE 0.1-3: NON-LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR (CONTINUED) 
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FIGURE D.l-3: NON-LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR (CONTINUED) 
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FIGURE D .1-3: NON-LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR (CONCLUDED) 
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FIGURE D . 1-4: LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR 
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FIGURE D .1-4: LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR (CONTINUED) 
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FIGURE D.l-4: LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR (CONTINUED) 
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FIGURE D .1-4: LINEAR EQUATIONS OF MOTION FOR INPLANE BIFILAR (CONCLUDED) 
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FIGURE D .1-5: INPLANE BIFILAR EQUATIONS IN FIXED SYSTEM COORDINATES 
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FIGURE D. 1-5 INPLANE BIFILAR EQUATIONS IN FIXED SYSTEM COORDINATES (CONCLUDED) 



D.2 Vertical Bi filar Equations of Motion 

The vertical bifilar rotor head absorber is a single degree of freedom 
system. The math model is presented in Figure 0 , 2 - 1 . The equivalent 
pendulum arm, r, and the distance from the center of rotation, R, are 
related by the tuning equation (18), shown in Figure D.2-2, The vector 
location of the dynamic mass (c.g.) is given by equation (19), shown in 
Figure D.2-2. Substituting equation (19) into the Lagrange's equation, 
also assuming small motions and neglecting second order terms, the linear 
set of equations of motion with six hub degrees- of- freedom are derived. 

They are shown in Figure D.2-3, equations (20) - (26 ) . 

For identical vertical bifilars, further simplification of these 
equations can be made by transferring the rotating coordinate, lnt0 
the fixed system coordinates. The equation for the coordinate transfor- 
mation used is shown in Figure D.2-4, equation (27). Thus the transformed 
equations of motions are derived, equations (28)-(36) , shown in Figure 
D.2-4. 
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FIGURE D. 2-3 : LINEAR EQUATIONS OF MOTION FOR VERTICAL BIFILAR 
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FIGURE D. 2-3: LINEAR EQUATIONS OF MOTION FOR VERTICAL BIFILAR (CONTINUED) 
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FIGURE D.2-3: LINEAR EQUATIONS OF MOTION FOR VERTICAL BIFILAR (CONTINUED) 
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FIGURE D..2-3: LINEAR EQUATIONS OF MOTION FOR VERTICAL BIFILAR (CONCLUDED) 
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D.3 Rotor Equations With Hover Aerodynamics 

The rotor equations of motion have been derived with five hub degrees 
of freedom (X, Y, Z, e x , 6y), as well as with four coupled flap-lag flex- 
ible blade bending, two torsional, a rigid flapping and a rigid lead-lag 
degrees of freedom. Definition of the coordinate system is shown in Figure 
D.3-1. A detailed derivation of the rotor equations of motion are pre- 
sented in Reference 1. For completeness sake, these equations are included 
herewith. 

The equations of motion are shown on pp 305 to 331. The airframe 
generalized coordinates have been utilized for the hub flexibilities. The 
first equation is for the fixed system, and the other four are rotor 
equations, i.e. bending, torsion, rigid flapping and rigid lead-lag 
equations (37) -(41) respectively. These equations consist of the acceler- 
ation, velocity and displacement coefficients which account for the left- 
hand-side of the Lagrange's equations of motion. The right- hand- side of 
the Lagrange's equations of motion, the generalized forces, are shown by 
equations (42)- (46) . In every equation, all generalized coordinates and 
physical properties associated with the rotor system have a subscript n, 
the blade number. The definition of all parameters used in all rotor 
equations are described in the list of symbols immediately following the 
equations (pp 332 to 338). 





Airframe Mode Equations 


+ *6Y,A + r6 0 ]) + ♦l.A.l - ♦«X,l Ib * + r U> 
+ ♦z.J^Z.i + +6X,i C(a 2 - rY o )sin * * (e + r * Vo ' Vo’ 003 * 3 


! + '* 9 Y,i l:(e + r * Vo ‘ Vo 131 ”* ‘ (a 2 ‘ r V 0 )o03 *] ) + * 0 YJ t(b 2 ■* r 6 o ) *X,i 

+ * z>1 [(e + r + o 2 y 0 - b 2 B 0 )sin» - (a 2 - rf o )cos*]) + + eXiJ (-( b 2 + r 8 0 ! *Y,i 

+ <J> Z i [(e + r + a 2 Y Q - b^jcosiji + (a £ - rY Q )sini| 0 ) 

+ + 6 Y,j U 0 Y > i C( '" a 2 ( Vo - Vo’ - <r + e)(a 2 - rY o ))3in2 ' 1 ' + ((r + 3)2 
+ 2 (r + e)(a 2 Y 0 - b^) )sin 2 i|/ + a 2 (a 2 - 2 rY 0 )eos 2 ij> + b 2 (b 2 + 2 r 3 Q )] 

+ ♦qX i C%(-a 2 (a 2 - 2rr 0 ) + (r + e) 2 + 2(r + e)(a 2 Y Q - b^) )sin2i|i 

+ (-(r + e)(a 2 - rY Q ) - a 2 (a 2 Y Q - b^) )cos 2 *]) + ♦ eX ,j ( +e Y , i Ps( “ a 2 (a 2 " 2 nr o } 
+ (r + e ) 2 + 2 (r + e)(a 2 Y Q - b^) )sin 2 tj> + (-(r + e)(a 2 - ry Q ) 


" %^- a oY rt ” b S ))cos2i|<] + * 0 v 4[(a (a y - b B ) + (r + e)(a„ - rYj)sin2<J> 


+ ((r + e) 2 + 2(r + e)(a 2 Y Q - b 2 B Q ))cos 2 ^ + & 2 (a 2 - 2rY 0 )sin 2 if» 

♦!> 2 < b 2 * 2re 0 )])>+ X x ar { t 9Y>J 4' E0 (- * 9Y ^os Vin2* ' *ex,i co3 V° 32 ' l ' ) 

+ +ex,j' l 'EO ( - 4>9Y,i cos V 032 '»' + * +0Y,j ( *er,i (co328 o 3ln2 * 

+ (y o cos 2 0 o - J se o sin20 o )sin2ij») + $ QX i (^cos 2 0 o sin24i 

+ (y o cos 2 0 q - * 5& o sin20 o )cos2i |i)) + $ex ^♦0Y,i^ cos2 % sin2,i ' + (y q cos 2 0 o 

- %3 sin20 )cos2t|>) + <}>„„ . (cos 2 0 cos 2 i|> - (y cos 2 0 - %B sin20 )sin2i|>))} 

O O va> 1 o o o o o 

+ 1 Y dr { ♦eYJ ( *0Y,i (q ’EO CO39 o - i'fo 31 ^) 81 ' 12 * * ♦ex.l^'EO 003 ^ 

- 4 'ro 3in, o )co32 * ) ♦ ♦ex,J ( %,i tq 'E0 oo3l, o ' q, FC 3ine 0 )co32 * 1 


305 



+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 

! 

+ 


^0X,i^'" <1, EO cos0 o + f poSi^lsina*) + ♦ 0 YjJ (4 0Y4 (cos 2 V» - Y 0 sin2*) 

^ 0 X i^“ V*? 82 * “ + $ QX *(- Y 0 cos 2 i|> - %sin2)l>) 

+ 9 X,i (Sln2 ' 1 ' + V in2,| ’ >)) + I Z dr { ^0Y ,j^ , F0^9Y,i S8n8 o 3an21 ^ 

*9x,i 3ine o 0082W + *ex > J ,l 'po ( *eY,i sine o 0032 ' , ' - +ex,i sin0 o sln2 1' ) 

*0Y,J^0Y i( sin20 o sin2 + + ^Y o s in 20 o + %3 o sln20 o )sin2f) 

4>0 2ji ( i SSin 2 0o s in24> + (Y o sin 2 0 o + *$B o sin20 o )cos2i|O ) 

^0X,j^0Y,i^ in2e o sin2 * + (Y o sin 2 0 o + ^ o sin29 o )cos2^) 

4> fiY . (sin 2 0 cos 2 i|» - (y sin 2 0 + ^ sin20.)sin2ij/) )>l q. | 

Qa 9 X O O 0 0 0 J Il 

[i5^ 9X,i ^ 0X » jI FA + + ^Z,i^Z,j^S^ + j M A,J q j j 

| E i [^ R ’ emdr ^x > j ((b 2 Y o + a 2 e o )sin * “ * 2 cos1>) + <fr Yj<3 (-(b 2 Y 0 
a 2 B o )cosiJ< - b 2 sini|») + <l> Z J a 2 + *0Y # J^' ft 2 + b P Y o + a 2^ r + e ^ sin,{ ' 

(a 2 U 2 - rY 0 ) + b 2 (b 2 + r$ 0 ))eoB*] + * 0 x,jf ^ a 2 (a 2 ” ry o ) 
b 2 (b 2 + r& Q ) )simj> + ((a 2 + b 2 )r 0 + a 2^ r + ■ e )) c08 *3J 

I X dr{ ^0Y , ( j q, EO COS0 o sin ^ + ^0X, j q, EO COS0 o COS ^ } + I Y dr{ ♦©Y, j ( ’ q 'EO COS0 o 

< ^ , FO S ^ n0 o )sin ^ + ♦0X,j { “ q, EO COS9 o + q, FQ sin0 O ) C0S ^ + ^0Y ,J (Y o sin ^ 
cosip) + <j> ex ^(Y 0 cos^ + sin^)> + I z dr { - <j> 0y jq^sii^sintji 

^X.^’po 3111 ^ 008 * 1 ] Vt| + j^ 1 [ / o'* mdrU X,J (b 2 + rB „ )8in * - *Y,j (b 2 

r6 o ,0 ° 3 * ♦ * 2jJ (r + a^ - b^) ♦ ^[(eU^ - b^) + r(r-+ a + Sa^,) 

^ ' 

b 2 )sl “* * l * 2 ( V, - Vo 1 + rta 2 - r Y 0 ))=°s*] + * 9 X ,j[ fa 2 (a 2 Y 0 ' b 2 B 0 > 

X. 

r(a 2 ” + (e ^ a 2 Y o " b 2 8 0^ + r ^ r + e + 2a 2 Y 0^ + b p COS '<'J } 

I X dr{ - *0Y, t J q, EO COs9 o CO8,, ' + ^0X , J q * E0 COS 0 o S ln ^ + ♦ 9YfJ Ccos 2 0 o3 iru| ( 2 


306 





+ (y cos 2 0 - 3in20 )cost|>) + (j» fl * (cos 2 0 costJ> - (y cos 2 0 

0 0 0 0 v7A >t J ° u u 

- ^ o sin20 o )0in«)} + I Y dr{ *eY,^ q, EO cos0 o " ^FO 311 ^ 008 ’*' 

+ *ex,j ( - i'eo cos6 o * i’ro 3ln0 o lsin ' 1 ' + ♦ex.jV 1 " 1 ’ ' ♦eY.jV 03 ' 1 ' 1 

+ I z drt V, J <1, PO sin9 o <:os ' 1 ' - *ex >3 l’ro sin V ln 'f' + ♦0Y,j (sin2e o sln * 

+ (y sin 2 0 + sin20 )cos<J>) + 4> 0 (sin 2 0 cos^ - (y sin 2 0 

0 0 0 O UAjj O v v 

+ ^ o sin20 o )sint)}] B j + 1 " r ^ 0 ) sin ^ “ t* 2 Y 0 

+ r)cos^} + $ Y j(“Ca 2 Y 0 + r)sinij/ + (. a 2 - ry o )cost|>) + 4 2 > j* 2 & 0 

* *eY,j C( ’ b 2 <a 2 ' n o l * V 6 0 )sln ' 1, + (_(a 2 + r2)6 0 ‘ b ! ( V« + riJcose] 

♦ + ex> jC((a| ♦ r2)S 0 + b 2 (» 2 Y 0 + r))sin* ■*■ (-b^ - n o ) + a^cos*]} 

+ I x 4r<- *0Y,j1'EO Sltl6 o COS * + *9X,j ll 'EO sln9 o sln ' 11 + *6Y ,J fein26 0 s in* 

+ (%Y Q sin20 o - 6 o sin 2 0 o ) cosiit J jP^in20 Q cosijA - (%Y 0 sin29 0 - 

+ I Y ai ' { *9Y,j ( ‘ l 'E0 3in9 o + 1 ' fo CO30 o ) c<>3 ' 1 ' + +ex,j (_ l'E0 3i,l9 o 

- 1 ' po COS 9 o ^ 3 inl ^ ^ * I z 4rt " *0Y,J <1 'fO c ° 30 o co3 ' 11 + ‘ t, 9X,J q, FO COBe o 3in ' ( ' 

+ 4> 0 y jf” J 5Sin20 o sin^ - Osy Q sin20 o + 6 o cos 2 0 o )cosiJ>] + $ 0X j[- ? 5Bin20 o cost|i ' 

+ (%r o sin20 o + 8 o c°s 2 0 o >sinip3}jY| , ‘'|^23 j^ R-ejndr{ ^X j^E,i^ Y o Sin ^ ” coa 1>) 

i *™«L 


+ v^ i sin^ + ♦ Fji 0 o siru|>) + +y»,)^E i^“ Sin * “ Y q cosi<») - v i ± cosf 
" *F,i 0 o cos * ) + *Z,j*F,i + ♦0Y,jC ( *F,l ( Vo + r + e ) + *e iVo* 81 ®* 

* h 

* ( ~ *F,i (a 2 | rr o ) - *E,i (b 2 + r6 o» M3 *} + ♦»X,jt { *F t i<»2 * rr o ) * *S,lK f 
+ r8 „ ),3ln * + ( *E,iVo + ♦p.i^j^o * r + e >> cos «>] + ♦ 0Y>J v i>1 b 2 Sin < , 

+ ^XJ^iV 08 ^ + I x < ^^ 0 YjCcos0 o 8in^ + (y o cos0 o - 0 o sin0 o 
” ^ ' EO ^ cos ^ * p ,i + ^exjCcosV 08 ' 1 ' " (y o CO80 o " 6 o 8in0 o “ <1 ’EO )sin,< '^ , F,i } 3 


307 



+ I Y dr{U eYj q V os * - ♦ex,j q, P 0 sAa * )4, E,i > 

+ I z dr{* 0Y ^[- sine o sinij; -(Y 0 sin © 0 + 8 o cos 0 o + cl' f 0 )cosi|/]<J> 

+ ♦qx, 4^ sin0 o cos* + (Y o sln0 o + & o cos0 o + q’po) 3 ^^ ' E , i } ] <1 T, i | 

+ j« VWfa- *9Y.J I L ) + ♦z,A ]X a| + | y K { - Vj* 

+ ♦z.jV J + j^V,j I I. X I.j ^l^jf 

- Ce + r + a 2 y o - b^sin*) + 4> 0Y>1 ((e + r + a 2 Y 0 - Vo )cos * 

+ (a 2 - ry o )sini/>)] + O<fr 0y E(' r + e ) 2 + 2 ( r + e ^ a 2 Y o ” ^2^0^ 

- a (a - ry o )]sin2ij/ + [-(r + e)Ca 2 - ry Q ) - " b 2 B Q ) ]2cos2i|») 


FA " *0Y,jV 


R " e mdr{ 2 fi 4 . z ^[()> 0 Xji ((a 2 - ry o )cosi|i 


+ * av .([a o (a o y n - b 2 0 Q ) + Cr + e ) ( a 2 - ry o )]2sin2* + [(r + e) 2 


0X,i 2 2 0 

- a z Ca 2 - 2ry o ) + 2(r + eHa^ - b 2 8 o )]cos2ii> - a^ - 2ry o ) - (r + e ) 2 

- 2Cr + e)(a 2 y o - b^))] ♦ ^ex,J^0Y^ (ta 2 ( Vo " Vo } + (r + e)(a 2 

- ry o )]2sin2^ + [- a 2 ( a 2 " 2ry Q ) + ( r + e ) 2 + 2( ' r + e ^ a 2 Y o “ b 2 B o )]cos2ij> 

+ a (a - 2ry ) + (r + e) 2 + 2(r + e)(a y - b B )) + <{> 6X ^1" ( r + e ^ 2 
2 2 Q z. v * 


- 2(r + e)(a y 0 - * z Z Q ) + a, (a, " ry 0 )]sin2* + [Cr + ®K a 2 “ r Y 0 ) 


2 2 


+ a 2 C a 2 Y 0 - ■b 2 B o )]2cos2t(>)]> + I x dr{2ftq’ EO U 0y <(» 9Y ji cos0 o cos2t|> 

+ $ 0Xji cos0 o sin2Tj>] + < ^ 0 x >( j^0Y,i COS0 o sin21, ' + *0X,i COS0 o COS2 ^ 

+ fiU 9YjJ [^ Y>i (co S 2 0 o sin2^ + (2y o cos 2 0 o - 6 o sin20 o )cos2*) 

p 

+ <)> 0X i (cos 2 0 o cos2i|. - (2 y o cos 2 0 o - B o sin20 o )sin2i(;) - 4> 0X>i 3 ia 2e o J 

«• Cfj|Y ji ( cos2 ®o cos2 ^ - ( 2 Y o cos 2 0 o - B o sin20 o )sin2^) 

2 

+ ^ 0X t (- coa 2 6 o sin2^ - (2Y Q cos 9 0 - 0 o sin20 o )cos2 1 j,) + <j> 0Y >i sin2 0 o 

+ I Y ,Jr(2a( *8Y,J C *9Y,i <q ’EO COse o * iV^V 0082 * + ♦ex.l 1 * 1 'eO°° s9 o 
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+ 4 , FO sin0 o )sin2*] + <t>ex f j£*eY,i ( ~ <1 'eO COS0 o 

j + q FO 8in0 o )ain2,j ' + 4> 0Xji (- q , EO cos0 o + q’ FO sine o )cos2^]) 

! 

+ 2y q cos 2 ^ - sin2tp) + <j> QX i (- eos2i// + 2Y Q sin2i|;)] - - 

+ ^0x,j^6Y,i^“ cos2 ^ + 2Y 0 sin2^) + <|tg X i (2Y 0 cos2i|i + sin2<j>)])} 

* ^^W^sy j [ *6Y,i sln9 o cos2 ' 1 ' ” *0x,i sin V os2 ' l ' ;l * * 9 x,jC 

1 “ ^9Y , i 8: ^ n0 o 8 ^ n2 ^ ~ (j> 0x ^ i sin8 o cos2ij/]) + fl(<|> 0Y ,j[^0Y sin2e 0 sin2 ^ + (2Y o sin 2 0 o 

i 

+ 8 sin29 )cos2\j>) + <f> AY . (sin 2 0 cos2^ - (2y sin 2 0_ + 8 rt sin20 )sin2i(i) 

v O vA f X 0 0 0 O 0 

■} 

- ^ x ,i cos2e o ^ + *0XJ ^0Y,i^ sill20 o CO8 ^ “ (2 Y o sin20 o + e o sin20 o )sin2*) 

+ . (-sin 2 0 sin2i|> - (2y sin 2 0„ + 0 sin20^)cos2ij;) 

ua 9 1 O 0 O O O 


& 


mdr{28 (<J> X j (0> 2 Y q + a 2 8 0 )cosi|i + b 2 sin^) 


* +Y ,j«Vo * V' o )sin * - V 08 ^ + *M„J ( V b 2 

+ r6 0 )sin* + » 2 (a 2 S 0 + Vo 1008 * 1 + +0X,J ( ‘ b 2 (b 2 Y o * * 2 6 0 >3la * 

+ ^2(^2 + r0 o )cosi|>))} + I x dr{fi(4> eY ^[cos20 o sin4» + Cy o cos 20 q 
- 8 0 sin2© 0 )cosi|>] + $ 0X j[cos29 o cosij> - (y o cos20 o - 0 o sin26 o )sinij> 3)> 
+ I y dr{Q((|» 0Y j ( y q cosi|> + simf>) + <j» ex ^(cos^ - Y 0 sin4»))> 

+ I ,dr{8 (<j» 0Y [-cos20 o sini|» - (y o cos20 o - $ o sin20 o )cosi|/J 

z 

* * ex jC - cos20 o cost|; + (y o cos20 o - S o sin20 o )sin<j»]) } J $ g § T 

* tin 8 6 mdr{2ft(^ x ^ (b 2 + rB 0 )cosij> + ♦y^Cbj + r8 0 )sin* 

+ ^0Y.,^ b 2^2 + 2^6 o^ COS1,, " *0X,j b 2^ b 2 * 2r ^g ) s i.mp ) } 




+ I x dr{fl(^eY ^cos20 o costj> - $ 0X jCos20 o sini|») > + I Y dr£8(<j> 0Y jCos^i 5 
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- <P eX 'j 8 im(>)} + I^drffiC - <f> 0Y ^cos 20 o cosij> + <j> 0x jCos 20 Q sin^) }J g | 

+ | £ [ ^ em< ir{2fi($ x j[ - (a^ - rY 0 )cosi// + (a^Y Q + r)sinif/] 

+ ♦v ~ ( a Y + r)cosi|» - (a - ry )sin^] + <p .[(r 2 B n + b (a y 
1 >3 2 ° 2 0 0* »J 0 £ 2 0 

+ r) )sin«j> - (b (a - ryJ + a rB )eosi|»3 + <J> , CC^B + b (a y •*• r))cosiJ> 

2 2 u 2 ° vAjj o 2 2° 

+ ^ 2^*2 * + » + I x dr{2ftq' EO C^ 0y ^sine^iruj* 

+ 4 > qx jSin 0 Q costJ>) + n(<J> 0Y j[sin 20 o cosij/ - (y o sin 20 o - B 0 )sinii >3 
+ 4» 0X j[ - sin20 o sinij) - CY Q sin20 o - B 0 )cosif>])} 

+ I Y dr{2n(4» 0YJ ( - 4 , E oSin0 o - q' po cos0 o )sin* + * 0XJ ( - q'^aii^ 

- q' FO cos 0 o )cos^) + ft(-$ 0Y jB Q sin^ - <j) 0x j 3 0 cosi|>)} 

+ I z dr{2nq , FO (<j» 0x ^jCos0 o sirnj) + <j» 0x ^cos8 o cosi/>) + ^(4> 0Y j( - sin20 o cosij» 

+ (Y o sin20 Q + B 0 )sini|)) + * 0Xi j(sin20 o sin* + (Y o sin20 Q + 3 0 )cos*)>] Y j 





i=l i=] 


jT^ £ ~ e mdr{2n($ x j^e,!^© 003 ^ + siruji) + v^cosij; 

+ ♦y i 3 Q cos^ + ♦ £~ ♦ E>i CcoB# - Y 0 sini^) + v^sin* 

+ ^,i 0 o sin ^ + %„>I>E,i Cb 2 + rB o )sin * + b 2 ( ^F,i 8 o 

+ h t i y o )c ° 3 ^ + ♦ex.j 1 - b 2 C *F,i e o + *E,i Y o )0ln * + ^,1^2 + r0 o )cos ^ 

+ ^ 0 Y,J V j,i b 2 C0S ^ ’ Vj V l,i b 2 Sln,,,)} + I X dr{flU 0 Y > j' [(sin 0 o COS ’ 1 ' 

+ Cq , E( ^in 2 e 0 - Y o sin 0 Q + 3 o cos 0 o )sin*) ♦ 

+ (co80 o co6f + (q* E0 - Y o cos0 q + 3 o sin0 o )sini|))* , Fji ] 

+ * 9Xjj r( - sin0 o sini^ + (q , B£) sin20 o -Y o sin0 Q + 3 o cos0 o )cos^)<^^ i 

+ C - cos 0 o simj> + (q'^Q - y o cos 0 o + 3 Q 8 ln 0 o )cosij/)<j)^ ^ 3 )} ■% 
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♦ I Y to «(* 9 I,J £*i »i^ sin 0 Q COS<() C 2 q pqCOs 0 q q gQsin 20 Q ^ 

~ 8 o c08 ®o * V in 6 o )8in <' ) + *'fJ co 86 o cOB * + , t - q’ FQ sin 20 o + q'^cos^ 

- 8 o sine o " Y 0 co se 0 ) s in^)] + *0X,j^ , E,i Csin0 o 8in,p + C " Sq'^co^ 

- q’ EO sin20 o - B o cos0 o + Y o sin0 o )cos*} + *' p>i (- cosO^in* + (q* EO cos20 o 

- q* FO sin20 o - B o sin0 Q - Y o cos0 o )cosij»)]) } + I z dr{a(* QY ^[(- sin« 0 co*# 

+ ^- 2<l, p0 cos2e o + Y o sirl6 o + 8 0 cosQ 0 ) sin 'l')^ , E i + ( - cose o cos* 

+ (a , TO 8in20 o + Y O COS0 O + B o sin0 o )sin^)+ , F + <j> QX j[(.sin0 o simj» 

+ (2q. , po cos 2 0 o +Y o sin0 Q + 3 o cos0 o )cosi|;)^ , e ^ + Ccosd^iwfi + (q' FO sin20 o 

- 1 , eo co 820 o +Y o cos 0 o + & o sin0 Q ) .cos*) *^ #i D 


H HA r B - 

EEU mdr{ii *Z t ^ex t i ( - Ca 2 - r Y 0 )siri ^ - (e + r + a Yq 
n=l i=r 

- b^B^cos*) + * 0Y ^( _ ( e + r + a^Y 0 - b^& o )sin* 

+ U 2 - rT o )eos * )J1 ] 5 i j + j "a.A.jSj I 

+ | 2 |^ R "' e md r {Q 2 ^ j [ - Cb 2 Y 0 + a^g^sin* + b^cos*] 

+ K i^K Y n + *A)eos* + b sin*] + * fiY £(.a (a y + r + e - 2b 0 ) 

1 »«J 2 ° 2 0 2 oi,J 2 2 0 2 0 

- b 2 Y )sini|) + ( - a (a - rY_) + b (b + r3 rt )).c°8*] 

2 0 2 2 0 2 2 0 

* ♦«W r( ‘* ( S ‘ *». ) - b 2 (b 2 + rS o )),in * + o» t <y. + r ♦ * - 2b 2 6 0 > 

- l>|r 0 )eos*3)}+ I x dr<ilV E< J* eY>J = 08 e 0 sln* + ♦ 0x>J cos8 0 C08*) # . 

+ n 2 £* 0Y j(cos20 o cos<j; - (y q cos20 o - 3 o sin20 o )sin*) + * QX ( - eos20 o sin* 

- ( Y oco820 q - 0 o sin2e o )cos^)]> + I Y &r{8 2 £* 0Y ^ (- q , E() cos0 o + q’^sinB^sin* 

+ 0 9X,j ( " *' eQ CO * B 0 * q, F(f in0 o )cO8 ^ } + I Z dr{ft2<1, F 0 (" *0Y, t 3 Sin0 o sin ' , ' ^ 
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~ ij> 9x ^ain© 0 cos^) + ft 2 [ (P 0y ^( “ cos29 0 cosip + Cy o cos20 o - 6 o sin20 o )sin#> 
+ * 0X , j Cco 2 0 s in^ Cy o cos28 o 8 o sin20 o )cosiJ/)3} J I %* 


1 ST 


mdr{ft 2 [ - $ . (b + rB.Jsimp + <p v (b + rB )cosip 

A*J 2 O X S J 2 O 


+ *0Y,j^ r + e ) ( u 2 y 0 + r ~ b 2 6 0^ + r ^ a 2 Y 0 " b 2 V " b 2^ 2rS 0 + b 2 ^ s ^ n< ^ 
+ E - a 2 Ca 2 Y Q + r - b 2 B Q ) - b,^ + 2rB Q ) + r^Jcosf) 

+ *0X,J ([a 2 (a 2 Y o + r - b 2 e o } + b 2 C2r6 o + V ' r2 YoJ sin ^ 

+ E(.r + e)(a 2 Y Q + r - b^) + Ha^ - b 2 B 0 ) - b 2 (b 2 + 2rB Q ) Jcosi|i) J} 

+ X x dr{fl 2 q’^t- ♦ eY>4 co80 o co8 <p + % f j c os0 o sinf) + ft 2 [^ 0Y>J (sin 2 0 o sint|, 

+ Y„cos 2 0 cosip) + <p flY .(sin 2 0„ cosip + y„cos 2 0 sinip)]} 
o o u A 9 j o o o 

+ I Y dr{fl 2 [<|) 0y>J (q* EoC os0 o - q’ FO sin0 o )cosi|» ■+ <P 0XiJ (- ^ e0 ^sQ q 
+ q* TO sin0 o )sinij( +<J> eYj j( - Y 0 cosip - sin*) + <P 0X ^ (Y 0 sinip - cosip)]} 

+ I z dr{a 2 q , po (4> 0Y ^sin0 o cosi(» - <p Qx ^sin0 o sinip) + ft 2 ! <p 0y ^ (cos 2 0 o sinip 


+ y sin 2 0 cos ip) + <j> av .(cos 2 0 cosip - y sin 2 0 sini 

0 0 v7A,j 0 0 0 



B 


lir 


mdr{ ft 2 [ . ( (a - rY Jsirup + (a Y + r)costp) 

AjJ 2 0 2 ° 


+ ♦y E( a 2 Y o + r ) sin * - ^ a 2 ~ r Y 0 )cosip) + <p 0 Y ^[Cr(a 2 8 Q - b 2 Y Q ) 

- b 2 ( a 2 Y 0 + r) + a 2 ([ r + e]B Q + b 2 ) - r 2 B Q )sia* + (b^ - nr Q ) 

+ b 2 (a 2 Y ° + r) " a 2 (a 2 8 o ' r6 o ) + r 2 B 0 )cos«] + 4» eXjJ E(-b 2 (a 2 - ry Q ) 

- b 2 (a 2 Y o + *% + a 2 (a 2 8 0 ’ r6 0 ) ’ r %^ aia * * ^ a 2 6 o“ b 2 Y 0 ) ' b 2 ( Vo * T) 

+ a^( [r + e]$ Q + b ) - r 2 B Q ) cosip]]} $ 
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+ I x dr{a 2 q. , E( ^4> 0Y ^sin0 o cos^ * $ QX sin0 o .8ijnJ»') + 0 2 [ <j> 0Y ^(8 o cos 2 0 o cosi|> 

- ^in20 siruj>) + 4>. v .( - 8 cos 2 © sin^ - ^sin20 cosii>) 1} 

o OA,J o o o 

+ I Y dr{t! 2 t* 9Y>J (- <l' E o sln<S o - q' ro =o S 8 o )cos* + *ex ) j( q, E0 8in6 0 
+ 4 'io co ' e o ) ' ln * • *0Y,j B o cosl1 ' * *ex,j 6 o slnl| ' ;,) 

+ I z d Ha 2 q , FO ( | J) 0Yj jCos0 o cosi|) - <j> ex ^cos0 o sini|O + fl 2 [* 0y ^C^in20 o sin^ 

+ 8 sin 2 0 cosip) + (%sin20 cosi|i - 8 sin 2 8 Sini//)!} I y 

0 O 9X,J 0 O 0 J 

mdr{n2[<j> x ^( - < f , E,i^Y 0 sin 4' - cos tp) - V-^sin^ 

I " ^F,i^o sin ^^ + *y,J ^E,i^ sin ^ + ) + v lji cosij> + <J> y ^cos*) 

+ ♦«,J (( *F,i<Vo + r + * ' a 2 S o> ' *E,iVo )sln * + ( - *F,i (s 2 - 

* *E,i Ci 2 * r6 o ))cOSI <’ ) * +0X,j (U F,i U 2 " n oi - ♦e, 1 (1> 2 + r8 o ),aln * 

+ <* F,i ( Vo + r + « - » 2 9 ,> - *S,iV» IC ° B * ! ' ♦«Y,J v ,,iV la * 

i " ^ 9 x,j v i, i b 2 C o s ^ } ' + I X dr{fl2 ^0Y,^ " sine 0 8in * - ^Y o sin0 o 

- 8 Q cos0 o )cos^)<|i , E)i + 4> 0X>J ( - sin0 o cos^ + (y o sin0 o - & o cos0 o )siniJ>)<J> » E ± 3> 
+ I Y dr{a 2 [^ 0Yjj (^ E ^sin0 o sin4> + (- q' EO sin20 o - e o cos0 Q - q' FO cos20 o 

+ Y o sin0 o )cosi|;) + <j)‘ y (■- cos0 Q sin^ + (- q, F0 Sin29 o " 6 o sin9 o 

! * 

l 

♦ q' EO cos28 0 -Y 0 cos8 0 )cosi(.)) + ♦ eX(J (*' Eil (sin0 o co 8 (. + (q'^si^ 

* 

i + 8^cos0 ■+ q. , „ r .cos20 i , - y sln.0 )sinf) + .(- cos0 cos<|» + (- q , _ n cos20_ 

OOrU ooo o Q 

+ B o sin0 o ♦ q* FQ sin20 o + Y o cos0 o )sinij/) )]> + I z <lr{n2{> QY ^((cos0 o sini{> 

+ (- < 1 , E0 COs29 o + Y o COS0 o + B o si&e 0 ) ftoe * ) * , p, i + * , B,l <l, .FO CO#28 o COS * ) 0 



+ <j> QX ((coBe o cos** ^'k 0' os29 0 “ Y 0 cose 0 ” 8 o sin0 o )sinij»)<> , F ^ i 
^ ♦ , Eji q. , KO c os20oSlii^)]}jq l i, J * 0 1 


( 31 ) 
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/ Blade Pitch Equations 

sr ^mdrf^ i^ t 2 Y o + a 2 B o^ ain ^ “ b 2 C0S ^ + “ ^2 Y o 

x p 

+ a,8 )cos<p - b„sin^ } + a,«J> . + <p v .[((a \ + bf)y + a,(r + e)}simf> 


2 y z ,i 0Y,i 

+ ( - a 2 (a 2 - ry Q ) - b 2 (b 2 + rB Q ) )cosi|/j| + (* 2 (» 2 - ry^ 

+ t> 2 (t> 2 + rB Q ) )simp + ((a^ + b^y^ + a 2 (r + e))cosij>]> 

* *x qr ^ q *Ec/^9Y,i cos9 o 3 * n ^ * *ex, 1 cose 0 oos,, ' )) + I Y drU 9Y,i C - q 'EO c03e o 

* q 'ro sin9 o )sln * + *ex,i ( - q ’B0°° 88 o + q 'Ki 3in9 o )c03 ' 1 ' * *eY,i C V in ' 1 ' 

- cos ill) + ‘frex^V 208 ’*' + sin< ^ } + I z to{q!' , ro ( - <{> 0Yji sine o sini|/ 

- *8X > ^ 3in6 0 003 ' l ’ ), ] *6 q i I 

+ j [ ^ R - e marf ( a | ♦ *§)> + I Y *r{l>] +| ®T I + j M 1 L 1 * 6 PR 9 T j 
+ | ^ R-e mdr{(a| + b^y^ + a 2 r} + I^drtq'ggCose^ + I y dr{ - q , BO cos0 o 


+ q, FO sin0 o + V 


+ I z dr< * q , TO sine o }] <f > 0 3 j + j - M 1 L 2 tan6 3 ^e PR ^ | 

+ | [^ R - e mdr{b 2 r> + ^drlq^s^} + I^dxl-q.' ^sin^ - q'^cosS^ 

+ I Z dr{<1 F 0 COse o } ] j + \ “ M l L 2 tana ^0PR Y j 

+ j £ [i R “ emdr{a 2*F,i + V E ,i } + I X dr{q, EO <), F,i } + r Y dr{ - q Wii } 

1 I ( NE 

+ ^drlq'po*^^) J <P e q T>i ^ + jJ^ M iM*FPR,I 

+ [L 2 tana 1 /(R-e}J(j) ET ^ i ))<|) epR q Tji J 

4 


NA 


,?v 
.. M: 


LlA p 

£ |^ R " e mdr{2fi[4. 0Y i (a 2 (a 2 - ry^sin^ + a 2 (a 2 y o + r + e - to 2 B 0 )cosiJi) 
i 3 l n ** 9 

+ *0X,i^ " a 2^ a 2 Y o + r + e “ ¥o )si ^ + a 2^ a 2 - r y 0 ^ C08 ^l} 10 


314 




+ I x dr{2nq , E0 [(}> QY jCOsO^os^ - < *'e X}i cos9 0 sin, f'] + n[> 0Y *(- cos20 o sinif> 


+ 


- (y o cos20 o - $ o sin20 o )cos^) + <f» 0X ^(- cos20 q cos!J; + (y q cos20 o 

- 0 o sin20 o )siru|»)]} + I Y dr{ft[«J) 0Y (- q» E() cos0 o + q'^sinO^cosil/ 

+ ^0X,i ( ^ , EO° OS0 o ' * V in9 o )2sin ’ i ' + ^0Y,i (Y o COS ^ + Sini,) + *0X,t ( “ y 0 Sin * 

+ cos4<) ]} + I z dr{2fiq , FO [-(|» eYji sin0 o costi» + <|) 0Xji sin0 o sin^3 

+ ^[^qv . (cos20 siru/> + (y cos 29 - 0 sin20 ) cos up) + <j> Qv . . (cos20 cosip 

i * y 0Y,i o r Vo o o o r/ r 0X,i o 

- (V oa26 o - e o sin2e o )sln * ) J ) ]»A I + ! + C l L 2*9PS )4 T [ 

+ ! I / R ” e m&r{2ft[- b (b + rB ) ]} + I v dr{- ficos29 } + I v dr{- n} 

| |_0 22 0 X 0 I 


+ I z dr{ficos20 o }] d> 0 6 | + | - C i L2tan6 3 <(> 0pR 8 j + j [^ R ~ e mdr{2n[a 2 (b 2 + r0 Q )]} 

+ I x dr{- nsin20 Q } + I z dr{nsin29 o }j <j> 0 Y ( + 1 ’ C^tana^gppY j 

iNE 

Elf " e “<ir{2fl(- ♦ Pfl b 2 6 0 + ♦ Efci * 2 6 0 ) - 2^ ij .b 2 } + I x dr{fi[- ♦ ’ S . 1 sin0 <> 
i=l 


- ^’ F ^ i cos0 o ]} + I Y dr{n[()) , E ^ i sin0 o - fp^ 0080 ^ + ^drtfitfVg^sir^ 

+ +'F,i COS9 o ;|! ] °1 L 2 U FPH,1 

+ | + j J^ R ‘ e m<3r{n 2 [e(a 2 Y 0 - t> 2 S Q ) 

- b (b + rB ) + a 2 ]} + I Y dr{- fi 2 cos20 I + I_dr{n 2 cos20 } <p 2 0 [ 

22 0 2 X O L OjyT) 

+ j H K *Q 2 ) 0 T | + | K i L 2 i,,2 0PR e T | 

+ J [^ R " e mdr{fl 2 [a 2 (a 2 Y 0 + r + e) - b^b^ - Pa^)]} + I x dr{n 2 [q» EO coa0 <) 

- Y o cos20 o }} + I Y dr{fi 2 {- q' EC) cos0 o + q , FO 8 in0 () )} + I z <lr{a 2 [- q'^sin^ 

+ Y o cos 20 q ]}] <j> 0 B j + j - K^^tanS^QppB | + | [^ R_e mdr{Q 2 [- (a 2 + b 2 )6 Q 

+ ebj} + I x dr{- ft 2 B o cos20 o ) + I z dr{ft 2 B o cos20 o }j <p 0 Y | >n 
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( a + «Y )*■♦«• 0> 

2 o' E,x 2 


, (HE 

+ - K^tan^+^T + S [/ H ' e mdr{Sl 2 0 Fj . 

^ f / i^l 

+ (r + e)B )]}+I Y dr{- Q 2 $ 4 >' ,cos 9 } + I dr{fi 2 B U’ . cos 9 + <J>* sin 8 )} 
o a ohj,x o i o Jci 9 x o r s x o 

+ I z dr«! 2 [- *' E>i <i’ F0 =°s 2 e o + ( q ' E0 co S 29 o - B 0 sin9 o )*' F>1 ])] *„ ^ j 
( NE 

+ k5 ( “ W*m.x + j *0 (3&) 


316 



Blade Rigid-Body Flapping Equations 


icT a 
6=1 



R-e 


m dr{^x,i (^2 + re o^ sin * “ *Y,i ^ b 2 + r & 0 ) co3 ' 1 ' + ^Z,i^ r + & 2 Y o 

" b 2 V + ^Y.i^^Vo " b 2 B o ) + r(r + e + 2 Vo ) + b 2 )sin ' 1 ' + a 2 (a 2 Y o 

- \Z 0 ) - r(a 2 - ry o ))co S1 |»] + * 0x>i [(a 2 U 2 Y o - b^) + r( &2 - rY 0 ))sin# 

+ (e(a v - b 6 ) + r(r + e + 2a y ) + b 2 )cosi|<]} + I dr{q’ ( 

2°2° 2^ 2 a i!iU 

- <^0y iC ose o cos* + <f> 0x i cos6 0 3irn(() + <j> 0y ^cos^sinip + (y o cos 2 8 o 

- sin20 )cos<|>) + . (cos 2 0 cos 4» - (y cos 2 9 - sin20 )sinif>)} 

0 0 t 7 A,l O 0 O O O 

+ ^^eY.i^W 03 ^ q ’TO sin9 o - V 003 * 1 

♦ 9 ex,l ( ' l’EO e0 " 9 o + 'J’po 313 ^ + V 3in * ) + ^z 9l '^ , p0^9Y»i S ^ n9 o co3 ^ 

- <|>q X ^inS^in^) + <fr 0y .^(sin^sinijj + CY o sin 2 0 Q + %& o sin20 o )cosi|>) 

+ 4 > 0X i ( sin 2 0 o cO 8 ^ - ( Y o sin 2 0 Q + ^ B o sin 20 o ) sini |>) }Jq\ | 

+ J ^ R " e mdr{((a2 + h Z) y ^ + a^)} + I x <3r{q' E() cos0 o } + I Y dr{(- l , EO cos0 o 


+ q* po sin0 ft + y o )> + I 7 dr{- <l’ F n sin0 o > 


FO' 


>] vjrj + j - 


| ^ R ” em<ir ^ r2 + b 2 + 2a 2 ry Q } + I x dr{cos 2 0 Q } + l z dr{sin 2 0 o 


+ | M^L^tan 2 ^^ > + 


ir 


! mdr{- b 2 ( a 2 , “ nr 0 )> + I x drOssin20 o > 


+ I z dr{- J5sin26 o > 


l]r| + | MjL^tanSjtaMjT | 





+ r) + *I,iVo + V l,iV + I X dr{ 9 ’p,i COs 9 o > + r z dr{ - * , E,i 3 in 9 0 , J«T,i 

( NE i 

+ + IV “ a i /( R - e )] 9 ET , l » ; 4 ,i j 

+|^^ R " e mdr{2H[* eY>i ((r(r + e + 2*^) * «C> 2 r 0 - >> 2 S 0 > + 'b|)cosi(> 



- b 2 e 0 ) ” b ^)sini|» + (* 2 (a 2 Y 0 - b^) - r 2 Y Q )cos^)]} 

+ I x dr{ 2 Q(t , E ^HQY i co 8 e o sin 4 » + $ QX ^osB^os^) ii 

+ n ^ 9 Y,i^ COS ^ “ 2 ^ o cos 2 0 o - ^ 8 o sirx 2 Q o )airnj/) + 4 > Qx y (- sin i|i - 2 (Y 0 cos 2 8 t 

- ^ 8 o sin 28 o )cosi(»)]} + I Y dr{J 2 [«j> 9y ^(- q , EQ cos 0 o + q» FO sin 0 o ) 2 sinif> 

+ (<> ex,i ( " q, EO cos9 o + *V te0 « )2cos * 

+ < * > 8 Y,i^- 2 Y o sin * " c ° 8 ^ + ^ex.i^o 008 * + sin ^n} + 

" ^eY,i sin0 o sin '*' " lJ, 0X,i sine o cos ^ + ^^0Y,i^ COS ^ ” (2Y o sin 2 0 o 
+ 0 o 8in26 o )sin^) + <fr QX i (- sin* - (2Y Q sin 2 0 o + B o sin20 Q ) cost/; ) ] }J j 


I [!-' 


mdr{2flb (b + rB )} + I v dr{8cos20 } + I„<lr-£n> + I dr{ - Jlcos20 }|<j> 
2 2 0 X o Y Z o * 


+ j - C^tanS^pg 8 t J t j C^tan^S 3 i J + j [/ R -'mdr{2[)(b 2 (a 2 Y o * r) 


]+0 5 lj 


♦ r 2 e 0 )} ♦ ^{nfaq-g^ine,, +’e 0 - r 0 stn 2 e o )} + i Y 


- 2q' EQ sin0 


~ 2 ^ , F O CO80 o " 8 o )} + h^^FO 003 ^ + V in20 o + 0 o )} ] Y 


C L 2 tan6 tana 
12 3 


,R-e 


mdr{2fl[^ Eji (b 2 + r6 0 )J> 


> ; l # 

+ I x dr{n[*' E>i (q' EO sin20 o + B o cos0 o - Y o sin0 o ) + ♦ , Ffi (l , E0 + & o sin0 o 
- Y 0 cos 8 o )]> + I Y dr{Q[*' E4 (- q’ EO sin20 o - 2q' FO cos 2 0 o + Y o sin0 o - B o cos0 o ) 

+ *’F,i (<1, EO c0829 o ‘ *'F0 Bin2 *o " V O80 o - e o 8in0 o^ } 

+ I 2 dr{fl ^ , E,i (2<1, FO CO820 o * V ln0 o + 6 o CO80 o ) + ♦ , F,i (<1, FO sin20 o * V 0 **#. 


+ 8 sin8,. 
0 o 

l -L 7 Xf 







*«$r{fl 2 [* 0Yii (- b 2 (b 2 + 2r0 o ))cosi|» + $ 0Xji (b 2 (b 2 + 2rB o ))sin*]} 

+ j [{ mar{Q 2 [a 2 (a 2 Y 0 + r + e) - b 2 (b 2 Y 0 - Sa^)]} + I x <lr{n 2 [q' EO cosi§ 

- r n COB29 n * 8 n sin 20 J> + ‘I'eo 0080 " " < 1 , EO 8 in 0 O 3 } + Igdrf^tq'pQSine 


r 

0 1 

14 
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! + Y cos 28 - B sin28 

* 0 0 0 0 


1>]*9%{ 


-- V^SWt 


+ | [f R '* Bdr{(i2 C e<a 2 ’ , 0 ' b 2 6 0 ) + r(r * e) - 1,r6 0 b 2 + 2a 2 rT 0 - b P } 

+ I^dr{H 2 sin 2 9 0 ) + I^dr{- R 2 } + l2dr{fl 2 cos 2 9 o }Jr { + j ^l b 2 b8n2 ^3^j 

* W 


• V 


mdr{ft 2 [b (a - ry ) + a 6 (2r + e)]} + I Y dr{ 


- igfl 2 sin28 } + I_<ir{%ft 2 sin28 }| Y > + I K L 2 tan6 tana y 


H 


12 


+ *E,iV° * *F,i ( V« + r + * + 2b 2 8 o> - v l,i b 2 ]> 

i=l 

+ I x dr{- i sin8 o > + I y dr{fl 2 [f , E ^ i sin8 o + $ ’p^cose^} 

♦ I z dr<- n 2 [+’ F>i cos9 0 j}]q T>1 J+ j^E ^l b 2 ban8 3^FPR,i 


+ [L 2 tamx 1 /(R-e)> EI , il ))(J Til ! , Q 


; 15 
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Blade Rigid-Body Lagging Equations 


R-e 


i C-Ca 2 - rY 0 )sint|> - (a^ + r)cos*) + ^ Y ,i^“^ a 2 Y o + r * sin * 


@ 

+ (a 2 - ry o )cosi|>) + $z,i a 2 e o + * eY [eB Q - b 2 K a 2 ” rY 0 ))sin* 

+ (- b 2 ( Vo + r) “ (a 2 + r2)e o ) cos ' , ' ) + *ex,i C( V Vo + r) + (a 2 

+ r 2 )B 0 )sin* + ([eB 0 - b 2 3(a 2 - ry o ) ) cos* ) } + I x <3r{q’ E0 (- ♦ 9Yji 8ine 0 cos * 
+ * ex i sin8 0 sin*) + * 9Y ^ i (^in20 o sin* + ( 3 sr o sin20 o - 0 o ain 2 0 o )cos*) 

+ * 0X i (?ssin20 o cos* - (%Y o sin20 Q - B o sin 2 0 o )sin*) } + I Y ^{ < t»e Y>i C < l , E o sin9 Q 
+ q» FO cos0 o )cos* + * eXji C- q’ EO aiix0 o - q» FQ cos0 o )sin*> + I z dr{q* F0 ( 

- * QY .cos0 o cos* + * 0Xji cos0 o sin*) + * 0Y>1 (- ^sin20 Q sin* - fey o sin20 o 
+ 0 o cos 2 0 q )cos*) + * 0X i (- i^in2e o cos* + (%y o sin20 o + 0 o cos 2 0 o )sin*) }J 
+ | [^ R ’ e mdr{b 2 r> + + I y dr{- q' E() sin0 o - <l , FO cos0 o > 


£ R ” e mdr{- b 2 (a 2 


W OS V] Vt| + | - M^tam^^ |+ j[ 

- rY 0 )} + I x dr{^in20 o } + I z <ir{- ^in20 Q }] ej + j M i L 2 tan5 3 tan °l ij 
+ j r^ R_e mdr{a 2 + r 2 } + I x dr{sin 2 0 Q } + l z dr{cos 2 0 o >] Y j 

' inn 


+ j M^L^tan 2 ©^ j + 

+ I Z 4r{ *'l,i CO,9 o > ] *T,il + 


/ t- 

I SET 


lmdr{ ^E,i r “ a 2 V l,i } + I X dr{ 4 ,, F,i Sine o > 

i NE 

J J^ M i L 2 tana i ^ ^FPR ,i 


* It J t«» l /(R-«).3* 8l , >1 ))5 T>i | +| i ^[f ? ' e "^{2“C*er,t (a 2 B o sln * 

+ (r + e)a 2 & 0 cos*) + * 0X>i (-(r + eja^sin* + a 2 S o cos*)]} 

k. 

+ I x dr{af- >i 6 o cos2 0 o sin ^ " ♦QX,i 6 o COS29 o COS *^ + I Y dr ^ e o^0Y,i sin ' 1 ' 
+ ^ 9X>i cos*)> + I z dr{ft6 o (* 0Y>i cos20 o sin* + ^0x,i cos29 o cos ' , ' )} ] *i j 
+ | [^ R " e mdr{- 2«a 2 (b 2 + r0 Q )} + I x dr{nsin2e o > + I z <ir{- asin28 Q }] * Q § T 



+ - C 1 L 2 tMO ^0®T 


| 2n(r 2 0 o + t> 2 (a 2 Y 0 + r))} + I x dr{Q(- 2q' E() sin0 


0 + Y 0 sin2V 


- 6 0 )} + I y dr{ne 0 } +I z dr{Q(- 2q» p() cos0 o - Y o sin20 Q - 8 0 )}] 0 J 
+ { C_L 2 tan6 | + j (2^1^ + C^tan 2 ^)^ | 


12 3 

r R-e 


| i I^[i h '' ,nart - 2a(a 2 *E,i + r6 o*P,i> - 2nrv ,,i } * ^^'E.i^l'EO 81112 ^ 

+ 6 0 sin» 0 ) - ♦ , F ,i B o cos0 oJ } + I Y dr{n[^» E ^(- 2q* EO sin 2 0 o - q’ FO sin20 Q 

“ ®o sin8 o ) - + * , F,i (q, EO sin26 o + 2q’ FO cos20 o + ^o c03Q o ] ^ } 

+ I z <ar{^ , E>i (<l , FO 8in20 o + 0 o sin0 o ) + ♦* F#1 (- 2q' FO cos 2 0 o 

-» i ( HE : 

" e o co80 o ):i} J «T,1 ; + £ ( W“V*FPR,i + CL 2 tano i /(R ’ e):i4, ET,i )) ^,il 
J £[£ ” e nKir{n 2 C'> 0Y>i ((- b 2 ( a 2 Y 0 + r) - r 2 0 o )sin<|; +Cb 2 (a 2 - ry Q ) 

+ a 2 r0 o )cosi{») + ♦ 8Xfi ((- b 2 (a 2 - ry Q ) - a 2 r8 0 )sini{i + (- b 2 (a 2 Y 0 + r) 

- r 2 $ o )cosi|»)]}J ^ | * | ^ R “ e mdr{n 2 [- (b 2 + a 2 )8 Q + eb 2 ]} + I x <ir{ 


- « 2 8 o cos 20 o > + I z dr{n 2 8 o cos20 o }] * Q 0 T | 


+ < - K L 2 tano 4 n __0_ 
j 12 1 T 0PR T 


+ J [^ R_e 3hdr{n 2 [b 2 (a 2 - ry Q ) + a^Ur + e)]} + 1^* r{- 3sfi 2 sin20 o } 


+ < K L tanfi tana 0 > + 
12 3 1 l 


| [£ R " 6mdr {ft 2 [e(r - a 2 Y Q ) 


• Y l + K y Y 


+ I z dr{isft 2 sin20 o }J 8 j 

" b 2 r6 o^!_ + * z ^ fi 2 8 o sin 2 0 o }]y J+ J K^tan^ 

e ~ V' ~ ) ~ ” 

lM ^ Vr{fl2[ “ Wo - C) - a 2 0 o^ F ,J }q T,i{ 

+ | ^ ( W““l U FPR.i + C L 2 tan “l /(B - e)] *ET,i )),5 T,i | » 0 M 


■W55»t 
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Blade Bending Equations 


j ^ to ^X,i^E^ Y o sin ^ “ cos ^ + v 1 > j sin '*’ + ♦ F ,j e 0 sin ^ 

i =1 

+ ♦y.J+E.j'- V° 3 * - sin,,) - v ul cos * ' VA 303 * 1 * ♦z.i^F.j 
+ *eY,i [U E,J b 2 Y 0 + +F,j ( Vo + r + e) > sin * + <- *F >J (a 2 - 

- * EjJ (b 2 + r8 0 ))cos«.] + +ex,i rt *F,j (a 2 ‘ ^ * +E,J (b 2 + r6 o ))sin * 

+ + *F,J ( Vo + r + e >5= os *J + *0Y,lVl„) 3ln * * ^ex.l'Vl ,J 003 ** 

+ I x dr{<j> , F ^[<t' 0Yji (cos0 o sinij» + (y o cos0 q - 6 o sin0 o - <l’ E0 )cosi|>) 

+ ^ 0Xji (cos0 o cosi^ - (Y O COS0 O - 6 o sin0 o - q'^isin^)]} 

* V’^Ej^.i^FD' 0 '* - +0X,i' l, FO 3ln *3) + I z ar»' EjJ [* M ,i<- si°9 0 sin* 

- (Y o sin0 o + 0 Q COS0 O + q' F 0 )cosi|/) + <(» QXji (- sin0 o cos^ + (Y Q sin0 o + B 0 cos9 0 

+ i' w )<ta*)]l] \ | ♦ ) [^••«artb 2 * IjJ + a 2 + F ,j } * V^Wr,,) 1 

+ I Y dr{ - ^'fO^ *E,j ^ * ^Z^ r ^ , F0^ , E,J Vt| + j ■ M 1 L 2 * 8 PR U FPE ,3 

+ 0 T |+ j + r) + ♦k.jVo 

+ Vl.J 1 + *F , j COs8 o^ + I Z ar{ - *E,j Sin6 o 1 ] 6 |* | M l L 2 ta “S U FPR,J 
+ [V““ 1 /(R ’ e>] *ET,j ) 8 j + | [i R " endr{ +E,j r ’ Yl,j’ + V rl+ *F.J* la# o > 

+ I Z drt +'E > j COse o , ] Y | + ) M l L 2 tanC, l ( 4 , FPR,J + [ L 2 tan “i /lR ’ e) ^ET,3 ) T | 

* j ^ [i S ‘ e mdr{+ 2 EjJ + 0 2 f>j } + I x drU’ F ,j + , r ,i > + V r{ *’B,i +, E,i ) J «t,tp 

x ~f- NE } 

+ J ( ^FPR,j + C L 2 tan 0 t i / ( R_e ) ]+ET , 4 ) (M 1 ^ *FPR ,i ^ 

4 „ . i NA 

+ [L^an^/CH-e)]^^)} J +^E [^“W{2R[<fr ey - ry^in^ 

+ Vj ( Vo + r + e ’ Vo^ 08 ** + ^ 9X ,i ( - *F,J { Vo + r + e - Vo )sin * 

48 

* ^F,J <a 2 ' rY o )cos * ) ] ) + I X dr{n I*'F,J ( +«,l (oos8 o eos * - ( V os *o - I'eO 


; ;3? 
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i ** 3 0 sin© 0 )3im(() + 4>ex»i^“ cos0 o sini|» - (y o cos0 q - q' E0 - 3 o sin0 Q )cos!j() ) 

> + sin0 o cos * + ( V in9 o “ 8 o CO30 o - q, EO sin20 o )sin,, ' ) 

I 

■’ + ^9x,i^ sin0 O s ^ n ^ + CY o sin0 o " e o cos0 o “ q' EO sin20 o )cos^))]} 

! + J X ‘ir{«C* , Ef< jC0 eY>i (sin0 o cos* + Cq ? EO sin20 o - 2q' F() sin 2 0 o 
' + 6 0 COs9 0 - Y o sin0 o )sin^) + ♦ 0?>1 C- Sixain* + (q' E() sin29 o - 2q' F() sin 2 0 o 
+ 6 q cos0 o - Y o sin0 o )cosi(»)) + ^ 0Y>i (” cos9 o cosi|> + (- q» E() cos20 o 

| 

+ c l'po sin20 o + 9 o sln0 o + Y o cos0 Q )sin(|)) + , *’Qx,i^ COS0 o sin '*' + (- q’ EC) cos20 o 
+ ^F0 Sln29 o + 6 o Sln0 o * Y o cos0 o )cost|())]} + I z dr{ft[«(* » p j(<J> QY ^ cose^osrp 
"( q, FO sin29 o “ q, EO COs29 o + Y o COS0 o + ® o s in0 o )sini|>) + <f> QXji (- eoa6 0 ai»f 
; - (q *F0 sin29 o ■ q * eq cos20 o + ycos0 q + S o sin0 Q )cost|») ) + <j>’ E j((- sin0 o cosi|> 
+ (2q , FO sin 2 0 o + Y o sin0 Q + B o cos0 o )sini{i)<j> 0Yji + <j> 0x ^ i (sin0 o sini|> 

+ ^ 2q Vo sln29 o * Y 0 Slne o + 8 0 cos9 0 )° os t'))]>] li | + j [i R ” e ' nar ^ a8 ^F,j 1> 2 8 o 

* ♦l,jVo ) * “VjV + 1 X arta( * , F,J cos9 o + *'E,J 3in9 o )) 

+ I y dr{n(- 4 »' e jSin0 o + 4> * F jCOs0 Q )} + I z <lr{n( 

- +'EJ 8in9 o * *'f,J cos9 o»] + j - C iV«R ( *FER.J 

+ rL 2 tam i /(R-e)]* ETjJ ) © T | + J [^ R ‘ e mdr{- 2G* E ^(t> 2 + r 3 q )} 

+ I X dr{n[ * , F,J C V os9 o ” V in0 o " q, EO } + *'E,J C V in0 o “ 8 o COs9 o 

! ~ ( l , EO 8in20 o > 3 > + I Y dr{n[ * , E,J (<1, EO sin29 o + 2q ’ro cos20 o “ Y o 3in9 o 

I 

+ S 0 cos« o ) ♦ f p>j (- 4 , eo cos 20 o + 1 , po 3in20 <j + t 0 »»V 8 0 3in « 0 >3> 
I Z '*«[*’ E ,j ( - 2<i' FO = 03 / 6 0 - V ln9 o * 6 o°° 89 o ) + *’F,j ( - q, F0 3in26 o 19 



+ 4'eo cos26 o ' V OS0 o - V in6 o )]) ] i \* ] C i l 2 t “ S 3 ( *fPR,; 


♦ [l^/cr-.)]*^) ej + J(/' 


-aar{ 2 nCa 2 * s>J + rB„*^ + v^r)} 


+ I X dr{a [*'p > j e o oos6 o * ♦’E,J (e o sln0 o + 2 l’EO sl “ 2e o)3> 

+ I Y dr{0£*' KJ (2 4 ' E0 sin 2 e o + q ' p0 sin26 o -f e o sin9 0 ) i r F>j (- l' E0 sin26 o 

- 2 4’po co * 2e o - S O =O80 O )]) + I z dr{n[0' E>J (- q' ro sin29 0 - 8 0 =ine o ) 


+ ♦r,j (2 » , FO cos2e o + «'EO ,t,,29 o + \* j 

+ e o^P*J*E,i ' + I X dr *°^o 

+ ^ * E0 S 1 n6 ° ^ ^ ' F , j ^ ' E , i " 

+ I Y dr{flC<1, EO sin0 o + ** fO cpS0 o + V^EJ^F,! ‘ ^FJ^E,i. )} 

+ I Z ar{!!( '»'E0 Bin9 o - 2 «Vo cos6 o - 6 0 )( - ♦'f.A*. 1 + '♦'E.^'F.i 11 ]^,! 

i NE 

+ J C i^FPR,j + I L 2 tana 1 /( R “e)J^ ETf j ) ^FPR.i 

+ EV“ c ‘ 1 / (B - e) ^ET,i ) » T ,i | + | 2 { 4 ,A. A.J 4 T,j I 

+ j [^W{B 2 [* p>J (a 2 + %> - ♦E,j ( *o (r + e) + b 2 )]) 

+ I x dr{- n 2 ^» E ^(<i’ E0 8in2e o + e o cose o )} + I y dr{fi 2 |> ' EJ (q , E0 8in26 o 

+ q’ FO cos26 o + 6 O COS0 O ) + *' E0 c ° s2e O + ^F0 S±n2Q o + e 0 sine o )]} 

+ I z dr{G 2 [- r EJ q* FO cos26 0 + <l , ¥0 ^e Q + q’ EO cos26 o 

-» 0 »i»e 0 )])]v T j ♦ j -VzW+fpej * [V“°i /(ft<,3 *ET,j )9 i | 

+ j ♦e,j» 2 y 0 - ♦y,j(Vo + r + e + 2 V,j * VjV» 

+ I x dr{- G 2 4 »' e j sin6 0 J + Iyflr.{O 2 (4' , E ^fli»0 o - <f> 1 p ^cos-6 p ) > 20 
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I 

- VriO^ecV] 6 j - | K 1 L 2 ta B « 3 (* ppR>J ♦ [L^/lR-I)]^) 8 J 

/ — ^ / £ 


+ ) K L tarn (<fc . 

) il l MJ 


' * j [/'Wfi 2 [- ♦e.j'Vo - •> - 'f,j‘A ] >] 

♦ [L 2 tana i /(R-e)]* E1 , i3 )Yj+J £ [{ S_e »to{- BVo V 2.i.J } 

♦ i x ar(fl 2 »i» 2 e 0 ** EJ *’ Etl } ♦ V*{a 2 [*- E>J (- ♦' E ,i Bin2e o ♦ ^' F ,i si,l 29 o> 

* ♦’ F> J (,s+ ’E,i sln26 o - ♦'p,i cos 2 e o );l) + I z 4 r( 82 C- ,s+ ’E,j*'F > i 8in2e o 

* fT *“■ 1 / (R -' ) ^ET,J ) £ <*FPR,i + CL 2 tana i /(E -e ):l * E T,i)' J T,l| 

n ( / 


+ [L tana 
2 


+ = O 


w) 



cLT 

* Jspcu£c l u t - 

C D U p] 


<3H 

« * 5 PeU[C L U p + 

OjPJ 


dM 

* *gpc 2 U 2 C w 

M 



3 ( dTT 

3% 

* 3 5P c1 / u Ic l (U' 

! + Uj) - CpU T U p 

+ °L,a U T U P 


- c D>o u|] + Hp= 1 / t [c l>m i4 - C, 

),M U T U P^ 

3 (dT) 
3U p 

* W/JCjVfy - C^U 2 ♦ u£) 

- C L,a U T 


+ C D,a U T U P^ + 

’“"’V'l.yv'p 

- C D,M U P 

3(dT) 

3a 

* ! “ x!U f C L,a U T 

- S.aV 


Kffl) 

SU T 

= j 5 Pp 1 / u Cc l u t u p + C D (U 2 + u 2 ) 

+ C L,a U l 


+ C D/l°p] * 

^At^.mVp 

+ c B y^ 

»w 

3U p 

- ! spc 1 / u [c l (u : 

2 + HP > + Wp 

- c i,.¥p 


- C D.a U P * *<*V C L,lrf + C 1 

3,M U T U P^ 

3(dH) 

3a 

* *" :U I C L.a U P 

+ C D,a U T^ 


XdM) 

3% 

« 3 5PC 2 [C m 2U t ■ 

" C M,a U p] 



* ’ 5(c2l/ v C M,M U T U 


3(dM) 

3Up 

= 3 5pc 2 [C J1 2U p ■ 

' C M,a U T ] 



+ 


1 sf iI = ’ spc2c M.a U 


2 



I 6 u t 


1 «U p 



a^T ajjr 

HA J : 

U M,i (U P C0E * - «ptv 0 - I'bjIsW) 

♦ex,i (u pt T o ' 'J' E0 ] cos <' + y lo *))«i 
{n(ap 0 + bq’ £0 )}* 6 9 T 

<“<» * >*„> + y-r 0 - «•„)}* 

- (»«<*„ - iV + u p S o h 

r WC* , s . i (- »-«»„♦ <r + «>4 ' E0 > - ^E.I^'eO - ♦f.A 

Tj] - u p[*' E(1 (e o + qy) + +' f> 1 (y 0 - 
HA 

Ji { ^0Y,i [b(Y o “ ^'EO )sin ^ - (b + r6 o )cos ^ + ^ex,i Cb(Y o " ^W 008 * 
(b + rP 0 )sini|/] + ♦ x>i t(Y 0 - " cos ^ + *y,1^ Y o 

q’^cosip - siruj/]}^ 

+ fl>H 0 e T 

+ {b ^0 - 4'so nl 

+ {r * < “J’eo , '» 

NE 

+ t I Hjfi'to * ♦ E ,i , %,i 

:j 

8Up 3Up 

“ ai"~ + aT~ 

- I^+ey.il- V“„ + iW ln *> - ♦ex.ilVo + 

+ <®>C * 0 + iVVi 

- {B(a - ry o ) + u p (8 o - d' ro )>6 

t 
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+ {n(rB o - eq' po )}Y y 

NE 

+ CO a + ry o ) + <i>E,i (6 o + " U P*’f,:L (B 0 + q ’FO )}q T,i 

NA 

+ £ {(-(a - rY Q )cosi{/ + (r + e + t>q' F0 + aY o ^ Sin ^^0Y,i 

+ ((a - rY 0 )sinij> + (r + e + bq' F0 + aY 0 ^ cos '^8X,i 

+ *x,i (6 o + »W sin * " *r,i (S o + ^fo 1002 * 

+ {a)«j> e e T 

+ {ay 0 + r + bq' F0 }B 
- 

HE 

+ <_ + j 

«• ■ H + If - ^Vt 


I K KA 

^ R "* [ { U 8X,i r *9Y,j (< - Vo + V’eO + r + e)dli 

+ a i (6 o + + - ♦y.jCVi® 3 ^ 

+ ^ex,j^" a 1 sin ’^ - a j( Y o ~ q' E 0 )cosi|/)dH - (bjSin^ - (a^^ 

” t > i Y 0 ) c °s^)dT] + «(i 0Y ^[(a l cosi|; - a^Yg - q’ E0 )siw(»)dH (b^osil* 

+ ^ a 1 q, F0 ■ b l Y o )sin ^ )dT 3 } ^e e T 

..^*7 

+ ^ex.i^Vn - r - a i^ , w>) 8in * + - q'^Jcos^dH 25 


0 X,j *• ' l o 


o ^ EO 



+ r 


- ( a j ( + e(6 0 + q'^cos*)®-] + * 0y [((- b^ 

+ V'*> )e0B * + e ^ 0 - 4 , B0 )si»*)m + (a i (0 o f q' F0 )cos^ - e($ Q 

+ <l , ro )ai^)dTj + * XfJ sin<fidT - <|> Y>J cos*dT + ^ z>j ((y o - q’ E0 )dH 

" (» 0 + ^'F0 )dT)}$ * ' ‘ 

* { *6X,J C( V Y 0 " q V sin * + (eS 0 - V COS * )dH - (O^q’po + a iYo 

+ r)sin* - (a ? - ryJcos^dT]-. ♦ ey>J [(b i (Y 0 - q’^Jcos* - (eB Q 

* 

- b i )sim|/)dH + ((- b^'^ - a^ - r)cos<|/ - (a - ry o )sin^)dT] 

' *X,/ ((y o * *'bo , co “* + sini|»)dH - q'^costdl] - * y>j [((y 0 - q’^fsin* 

- cos*)aH - q'ygSimi'aT] + * 2>J e 0 an)Y 


HE 


+ £l *F,i + V Y o * + a i B o*'E,i >sln ' 1 ’ 

- <*F,i (Y o - «») + V'l,i - ee o*'E,i )cMI|,)aH - ((*8,1 + b ,< Y o 

q EO^ F,i + a j^®0 + q, FO^'p,i + b i q, F0 9 'E,i^ sin '*' + ^ 4 E,i Y o 
' r i ‘ 1’F0*F,i + (r + e)q' P 0*'F,i )cosl| ' )dT 3 * *eY,J tl( ' *F,i 

+ V Y o - + WE,i )e “* * ( *F,1 (y o - 1’eo’ + V’s.i 

* e6 o*'E,i> al "*)' 3H + «* ♦'*,! * V Y 0 - »'EO>*’F,i * a j (S o * «W*’f,1 

* Vro*'E,i )c0S ’ 1 ’ * U E,i Y o * r i - «'FO*F,i + (r + e)q' F0 4.' F>i ) S in*)dT] 

’ *X,J C( * , E.l Sini1 ' + ♦' E ,i (Y o - 4 , EO >oos4,)dH + ♦ , F .i ,in * * (♦’ F ,i (Y o 

- *’*>> * 4’ P 0*'E,i ,eos,|,)4T J - ♦x.J 11- +'E,i cos * + ♦'E,i Clf o 

- ♦ (♦' F ,i cos * * <*’ F ,i (Y o - «’e 0> + q’po*'E,i )aln *> 4T ] 

- *Z,J C - Ve,!® + ♦'f,1 (6 o + l’F0 )4T l>'>T,i + C(* 6X ,(a i Y 0 * r ♦ e 26 
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LIST OF SYMBOLS 


= q E0 + (% - EA)cos0 o 
= q EQ - ACcose Q 
= q EQ - CGcos6 0 

" EA/b h 

Distance from blade elastic axis to aerodynamic center- 
positive toward leading edge. 

Distance from blade elastic axis to pushrod: positive toward 

leading edge. 

Radial location of blade pitch horn. 

Length of one arm of tail rotor blade pitch spider beam. 
Length of tail rotor pitch beam arm. 

Blade elemental mass 
Airframe mode generalized mass. 

Generalized mass of blade bending modes. 

Generalized mass of fixed system modes. 

Mass at pushrod. 

Mach number. 

Partial derivative of pitching moment with respect to local 
blade tangential velocity. 

Partial derivative of pitching moment with respect to local 
blade vertical velocity. 

Partial derivative of pitching moment with respect to local 
blade angle of attack. 

Blade number. 

Number of blades. 

Number of blade bending modes. 



Number of fixed system modes. 

Partial derivative of drag with respect to local blade angle 
of attack. 

Blade flatwise, edgewise, torsional mass moment of inertia 
matrix. 

Blade torsional mass moment of inertia. 

Blade edgewise second moment of area. 

Elemental blade flatwise mass moment of inertia. 

Elemental blade torsional mass moment of inertia. 

Elemental blade chordwise mass moment of inertia. 

Total blade torsional mass moment of inertia. 

Blade mass moment of inertia about lag hinge. 

Local blade polar second moment of area. 

Blade torsional stiffness. 

Airframe mode generalized stiffness. 

Blade lag hinge spring rate. 

Blade flapping hinge spring rate. 

Pitch beam stiffness or stiffness at main or tail rotor blade 
pushrod. 

Stiffness of tail rotor pitch actuator for pure moment applied 
at pitch beam end. 

Partial derivative of pitching moment coefficient with respect 
to Mach number. 

Distance from blade elastic axis to center of gravity — posi ti ve 
toward leading edge. 

Elemental lift. 

Elemental drag. 

Elemental pitching moment. 



dT Elemental thrust. 

dH Elemental inplane force, 

e Blade offset. 

EA Distance from blade semi-chord to elastic axis-positive toward 

trailing edge. 

G Blade torsional modulus of elasticity. 

h. Partial derivative of drag with respect to local blade tangential 

velocity. 

h 2 Partial derivative of drag with respect to local blade vertical 

velocity. 

q. Fixed system mode generalized coordinate. 

J 

Qj • Blade bending mode generalized coordinate--bending up and 
,J leading positive. 

q F Q Steady blade flatwise deflection- -up positive. 

Q eo Steady blade inplane deflection--lag positive. 

q QX Hub pitch coordinate. 


q eY 


Hub roll coordinate 


q” x Hub lateral coordinate. 

q Y Hub longitudinal coordinate. 

q" z Hub vertical coordinate. 

q Q Collective mode coordinate, 

q Q Reactionless mode coordinate, 

q s Sine cyclic coordinate, 

q Cosine cyclic coordinate. 

V 
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Blade flatwise bending slope steady, and sine and cosine 
coef f i ci ent components . 

Generalized forces. 

Generalized force on j'th fixed system mode. 

Generalized force on blade pitch. 

Generalized force on blade flapping. 

Generalized force on blade Tagging. 

Generalized force on j'th blade bending mode. 

Generalized aerodynamic force. 

Radius of local blade element from offset. 

Radial location of inner snubber of crossbeam rotor. 

Radial location of outer snubber of crossbeam rotor. 

Rotor radius. 

Radius to servo connections on main rotor swash plate. 

Radius to pushrod connections on main rotor swash plate. 

Partial derivative of thrust with respect to local blade 
tangential velocity. 

Partial derivative of thrust with respect to local blade 
vertical velocity. 

Partial derivative of thrust with respect to local blade angl 
of attack. 

Total local blade inflow velocity. 

Local blade vertical velocity. 

Local blade tangential velocity. 



Speed of sound. 

Potential energy. 

Rotor axial velocity. 

Forward-flight speed. 

Displacement at blade pushrod. 

Local blade angle of attack. 

Pitch-lag coup! ing— lead, pitch-up positive. 

Blade rigid-body flapping generalized coordinate— up positive. 
Steady blade coning— up positive. 

Blade rigid-body lag generalized coordinate— lead positive. 

Steady blade lag— lead positive. 

Pitch-flap coupling— flap up, pitch down positive. 

Fraction of critical structural damping of blade bending modes— 
based on modal frequency. 

Fraction of critical structural damping of blade pitch mode- 
based on rotor speed. 

Fraction of critical rigid-body lag damping— based on uncoupled 
lag frequency. 

Fraction of critical structural damping of fixed system modes— 
based on modal frequency. 

Blade pitch generalized coordinate— leading edge down positive. 
Steady blade pitch angle- -leading edge down positive. 

Blade pitch normal coordinate. 

Geometric blade pitch angle— leading edge up positive. 


(f,, E,i q, E0 + *'F,1 q W dr ) 

= + + , F,1+’F,j )dr) 

Air mass density. 



r i 


r i ,J 



Local blade inflow angle. 


<t> 

^X,Y,Z Fixed system translational mode shapes at hub. 

^eX.eY Fixed system rotational mode shapes at hub. 

cf>p Blade flatwise bending mode shape. 

(j> E Blade inplane bending mode shape. 

<f> 0 Blade torsional mode shape. 

<|> e p R Blade torsional mode shape at pushrod radial station. 

4>FPR Blade flatwise bending mode shape at pushrod radial station. 

Blade inplane bending mode shape at tip radius, 
i P Blade azimuthal angle. 

w T1 Blade asymmetric torsional frequency. 

Wj2 Blade symmetric torsional frequency. 

Wp Hub pitch frequency, 

coy Hub yaw frequency. 

Blade edgewise natural frequency. 

Wq Frequency of blade bending modes. 

Uncoupled rigid-body lag frequency, 
o Frequency of fixed system modes, 
w Flutter frequency. 

Q Rotor speed. 

Subscripts 

i Refers to blade element or mode number. 

J Refers to mode or force number. 
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s 

Refers to mode number, 
t 

n Refers to blade number. 

PR Refers to pushrod. 

X 

Y Refers to hub lateral, longitudinal, and vertical directions. 

Z 

9^ Refers to hub pitch and roll directions. 

9Y 

Superscripts 

c Menas coupled, 

u Means uncoupled. 

Differential Notation 

' Differentiation with respect to radius. 

Differentiation with respect to time. 

Second differential with respect to time. 
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D.4 Fuselage Dynamic Representation 

The fuselage math model is presented by equation (47), shown in 
Figure D.4-1. It consists of up to sixteen rigid and/or elastic mode 
shapes. For this modal representation we need the generalized masses 
(Mg-j), frequencies (w n -j) and damping (.q) for these mode shapes. We 
also need the modal components at any point on the aircraft where forces 
and moments are applied. 

Obtaining the generalized coordinates, q-j , one can then evaluate 
the response at any point on the aircraft using equation (48), see Figure 
D.4-1. 
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FIGURE D.4-1: FUSELAGE MATH MODEL 
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j D.5 Fixed System Absorber 

i I 

| 

j The math model of the fixed system absorber is shown in Figure D.5-1. 
i It is a one degree of freedom spring-mass system. The kinetic and potential j 
j energies of the system are given by equations (49) and (50) respectively, 
shown in Figure D.5-1. Substituting into Lagrange's equation results in 
the equation of motion for the fixed absorber shown in Figure D.5-1 by 
equation (51). j 
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Thus : 

+ [c]{x} + 9dM *o &) 


D, 5-1: FIXED SYSTEM ABSORBER 
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J- 5-1: FIXED SYSTEM ABSORBER (CONCLUDED) j 
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D.6 Assembly of Coupled Equations 

To assemble the coupled equations of motion we need to know the modal 
components, of the same mode shapes used in Section D.4, at any aircraft 
| point where the rotor, fixed and rotating absorbers are attached. 

Here we will demonstrate the procedure used to couple the fixed system 
with that of the rotor. This same technique is used to couple the bifilars 
and fixed absorbers with the fixed system. 

Equation (52) in Figure D.6-1 shows the rotor equation of motion in 
general form. We partitioned the mass, stiffness and damping matrices into 
submatrices associated with the hub (attachment point) and rotor degrees of 
freedom, see equation (53) Figure D.6-1. 

The hub degrees of freedom are related to the fuselage (generalized) 
degrees of freedom by equation (54) shown in Figure D.6-1. Thus equation 
(53) can be transformed into equation (55) whose state vector consists of 
the fuselage and rotor degrees of freedom, as shown in Figure D.6-2. 

With the hub degrees of freedom replaced by the generalized degrees of 
freedom, q's, we can now combine the rotor and fuselage equations of 
motion, as shown by equation (56) in Figure D.6-2. 
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Figure D.6^1: Uncoupled Equations of Motion 




Figure D.,d-2; Procedure to Form Coupled Equations (Continued) 
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/ 


/ ji j 

Figure D.,6-2; Procedure to Form Coupled Equations (Concluded) 
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D.7 Forced Response Solution 


The final math model is represented by equation (57) of Figure D.7-1, 
where the state vector, X , consists of the fuselage, rotor, fixed and 
j rotating absorbers degrees of freedom. 

A solution of the form shown by equation (58) is assumed which after 
substitution in equation (57) results in equation (59), see Figure D.7-1. 
Expressing the state vector, X c and X s in terms of the sine and cosine 
components of the fuselage and remaining degrees of freedom, see equation 

(60) Figure D.7-1, and substituting it into equation (59) gives equation 

(61) from which equation (62) immediately follows after grouping like 
degrees of freedom. The form of equation (62) is preferred over that of 
equation (61) since the size of the inverse matrices involved in the 
solution is greatly reduced. From equations (62) the solutions for q and 
Y are expressed by equations (63) and (64) respectively, see Figure D.7-1. 
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Figure D. 7-1: Forced Response Solution (Continued) 
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4 Figure D. 7 - l ; Forced Response Solution (Concluded) 
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D.8 Time History Solution 


We have shown previously the technique used to couple the equations 
of motion of the fixed system, rotor, fixed absorber, and the linear 
bifilar transfered to fixed system coordinates. The final equation we 
arrived at is of the form shown by equation (65), Figure D.8-1. We can 
rewrite equation (65) in the form shown by equation (66), Figure D.8-2, 
where the right-hand-side of the equation is replaced by a forcing 
vector. The non-linear inplane bifilar equations of motion, shown in 
Figure D.l-3, can be rewritten in the compact form shown in equation (67), 
Figure D.8-2. Using the coupling technique described in section D.6 we 
can also couple equations (66) and (67). The resultant coupled mass 
matrix and force vector are shown in Figure D.8-2 by equations (68) and 
(69) respectively. The state vector is expanded and consists of the 
fuselage, rotor, fixed absorber, linear inplane and/or vertical bifilar 
in the fixed system coordinate, and non-linear inplane bifilar degrees of 
freedom. The final equations of motion of the system can be rewritten 
as shown by equation (70) Figure D.8-2. The reason for partitioning the 
matrix as shown in equation (70) is that the submatrix D has been reduced 
to the identity matrix and the only inversion required for solution is 
that of a matrix whose dimension is much smaller than that of the total 
mass matrix as shown in equation (9). This saves considerable computer 
time and in addition it is more stable. 

Solving equations (71) and (72), shown in Figure D.8-3, we get the 
acceleration vector. The velocity and displacement vector are updated by 
equations (73). This procedure is repeated until the state variable, X , 
have converged within the specified constraints. 
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WHERE: 



Rotor Pe&Rets of Freedom 
X Ffl * Fixed Absorber Degrees of Freedom 
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FIGURE D.8-2: COUPLED EQUATIONS OF MOTION (CONCLUDED) 
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FIGURE D.8-3:' SOLUTION OF COUPLED EQUATIONS OF MOTION 



D.9 List of Symbols 
j Fixed absorber mass 

j Mq^ s mg-j Generalized masses 
My Total bifilar mass 

m Individual bifilar mass 

N Total number of bifilars 

n Bifilar tuning 

q.j Generalized coordinates 

R Distance from center of bifilar tracking hole to center of 

rotation 

r Equivalent pendulum arm 

Vertical motion of fixed system absorber 

Xpn Vertical motion of attachment point of fixed absorber and 

A fuselage 

X Rotating coordinate system 

Xj Inertia coordinate system 

w ^ Fixed absorber tuning frequency 

w. Generalized frequencies 

w D Vertical bifilar tuning frequency 

P 

w Inplane bifilar tuning frequency 

3 Vertical bifilar degree of freedom 

y Inplane bifilar degree of freedom 

^ Fixed absorber damping 

z >Q Vertical bifilar damping 

P 

Inplane bifilar damping 






XL • 

1 

Generalized damping 

1 

\ 

| 

* 

Mode shapes 



Bi filar arm angle of rotation 


ft 

Rotor speed 

1 

i 

9 x 

Hub roll 

3 

1 

| 

e y 

Hub pitch 


e z 

Hub yaw 

| 

c x 

Hub longitudinal modal damping 

! 

! y 

Hub lateral modal damping 

| 

? Z 

Hub vertical modal damping 

: 

1 

c ex 

Hub roll modal damping 


c ey 

Hub pitch modal damping 

:! 

■1 

c ez 

Hub yaw modal damping 

i 

w x 

Hub longitudinal modal frequency 

■{ 

! 

w y 

Hub lateral modal frequency 

1 

w z 

Hub vertical modal frequency 

i 

w 

ex 

Hub roll modal frequency 


w 

ey 

Hub pitch modal frequency 


w 

ez 

Hub yaw modal frequency 

■j 

Subscripts 

i 

H, h 

Hub 


HR 

Hub-Rotor 


i 

Fixed system mode 

4 

k 

bifilar 

i 

j 

R 

Rotor 

1 

i 

1 | 

1 
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Differential Notation 


Matri x 
C 

C HH 

C HR 

C RH 


'RR 


'qR 

: Rq 

'qFA 

TAq 

: qLB 


'LBq 


AP 


R 


K 


K, 


HH 


K, 


HR 


Differentiation with respect to time 
Second differential with respect to time 
Definitions 

Damping matrix 

Hub damping sub-matrix 

Hub/ rotor damping coupled sub-matrix 

Rotor/hub damping coupled sub-matrix 

Rotor damping sub-matrix 

Generalized fuselage damping sub-matrix 

Generalized fuselage/rotor damping coupled sub-matrix 

Rotor/generalized fuselage damping coupled sub-matrix 

Generalized fuselage/fixed absorber damping coupled sub-matrix 

Fixed absorber/generalized fuselage damping coupled sub-matrix 

Generalized fuselage/linear bifilar damping coupled sub-matrix 

Linear bifilar/generalized fuselage damping coupled sub-matrix 

Force vector 

Force vector at any point on the aircraft 

Cosine component of generalized force vector 

Sine component of generalized force vector 

Force vector at hub 

Force vector at hub due to rotor 

Stiffness matrix 

Hub stiffness sub-matrix 

Hub/ rotor stiffness coupled sub-matrix 



I 


'RH 


'RR 


qR 

^Rq 

( qFA 

( FAq 

( qLB 

( LBq 


M 

M 

M 

M 

M 

M c 

M 

f 

M 

M 

M 

M 

M, 


HH 


HR 

RH 


RR 


qR 

Rq 

qFA 

FAq 

qLB 

LBq 


Rotor/hub stiffness coupled sub-matrix 

Rotor stiffness sub-matrix 

Generalized fuselage stiffness sub-matrix 

Generalized fusel age/ rotor stiffness coupled sub-matrix 

Rotor/generalized fuselage stiffness coupled sub-matrix 

Generalized fuselage/fixed absorber stiffness coupled sub-matrix 

Fixed absorber/generalized fuselage stiffness coupled sub-matrix 

Generalized fuselage/linear bi filar stiffness coupled sub-matrix 

Linear bifilar/generalized fuselage stiffness coupled sub-matrix 

Mass matrix 

Hub mass matrix 

Hub/rotor mass coupled sub-matrix 

Rotor/hub mass coupled sub-matrix 

Rotor mass sub-matrix 

Generalized fuselage mass sub-matrix 

Generalized fusel age/ rotor mass coupled sub-matrix 

Rotor/generalized fuselage mass coupled sub-matrix 

Generalized fusel age/f i xed.absorber mass coupled sub-matrix 

Fixed absorber/generalized fuselage mass coupled sub-matrix 

Generalized fuselage/linear bifilar mass coupled sub-matrix 

Linear bifilar/general ized fuselage mass coupled sub-matrix 






